Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,823,581 G→A 5.2% V63V (GTC→GTT)  recA ← DNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,823,5810GA 5.2% 2326.9 / 68.3 691V63V (GTC→GTT) recADNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor
Reads supporting (aligned to +/- strand):  ref base G (340/315);  new base A (18/18);  total (358/333)
Fisher's exact test for biased strand distribution p-value = 8.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 734 total aligned reads displayed.

GATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTG  >  NC_000913/2823434‑2823729
                                                                                                                                                     |                                                                                                                                                    
gATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATGACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTcga                                                                                                                                                     >  4:4759471/1‑149 (MQ=255)
 atatCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTcgac                                                                                                                                                    >  3:1326888/1‑149 (MQ=255)
    tCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCAACGAt                                                                                                                                                 >  2:2272643/1‑149 (MQ=255)
         gCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGc                                                                                                                                            >  2:977745/1‑149 (MQ=255)
         gCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCAACGATACGGc                                                                                                                                            <  3:2319071/149‑1 (MQ=255)
          cccAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGcc                                                                                                                                           >  1:469380/1‑149 (MQ=255)
           ccAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGccc                                                                                                                                          >  4:3898929/1‑149 (MQ=255)
           ccAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGccc                                                                                                                                          >  1:6619522/1‑149 (MQ=255)
            cAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCa                                                                                                                                         >  3:4113045/1‑149 (MQ=255)
            cAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCa                                                                                                                                         >  1:2914428/1‑149 (MQ=255)
              gTTTACGTGCGTAGATTGGGTCCAGCGCGTGCTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATc                                                                                                                                       >  2:2785851/1‑149 (MQ=255)
                  aCGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCAACGATACGGCCCATCGGCa                                                                                                                                   <  2:4579596/149‑1 (MQ=255)
                             ttGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGcc                                                                                                                        >  1:1084345/1‑149 (MQ=255)
                                 gTCCAGCGCGGGTTCAACATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCaa                                                                                                                    <  3:6903410/149‑1 (MQ=255)
                                     aGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgcgc                                                                                                                <  3:1150621/149‑1 (MQ=255)
                                        gcgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGat                                                                                                             <  3:4194944/149‑1 (MQ=255)
                                           tgtTCAGCATCGATAAACGCACAGGTTTTACCTTCAGGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATc                                                                                                          >  4:6741883/1‑149 (MQ=255)
                                            gtTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgcgc                                                                                                                >  2:6957180/1‑142 (MQ=255)
                                            gtTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATcc                                                                                                         >  3:2653775/1‑149 (MQ=255)
                                             ttCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCa                                                                                                        >  2:5912934/1‑149 (MQ=255)
                                             ttCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCa                                                                                                        <  4:5923608/149‑1 (MQ=255)
                                              tCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAg                                                                                                       <  4:7071642/149‑1 (MQ=255)
                                              tCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAg                                                                                                       <  2:6789179/149‑1 (MQ=255)
                                               cAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGt                                                                                                      >  3:6273814/1‑149 (MQ=255)
                                               cAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGt                                                                                                      <  4:3579097/149‑1 (MQ=255)
                                                 gCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGATGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGa                                                                                                    >  2:518524/1‑149 (MQ=255)
                                                 gCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGa                                                                                                    >  1:1561691/1‑149 (MQ=255)
                                                   aTCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGcccc                                                                                                                      <  4:6104463/129‑1 (MQ=255)
                                                   aTCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGaaa                                                                                                  <  2:1329208/149‑1 (MQ=255)
                                                   aTCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGaaa                                                                                                  >  3:3689582/1‑149 (MQ=255)
                                                   aTCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCAACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGaaa                                                                                                  <  1:430168/149‑1 (MQ=255)
                                                      gATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCg                                                                                               >  3:3240877/1‑149 (MQ=255)
                                                       aTAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGa                                                                                              >  3:6090570/1‑149 (MQ=255)
                                                        tAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGaa                                                                                             <  1:2076382/149‑1 (MQ=255)
                                                         aatcGCTCAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCCTACATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAAc                                                                                            <  2:4165751/146‑1 (MQ=255)
                                                          aaCGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAAcc                                                                                           >  2:5726785/1‑149 (MQ=255)
                                                              cacaGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTa                                                                                       <  1:6872928/149‑1 (MQ=255)
                                                               acaGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTag                                                                                      >  3:6509740/1‑149 (MQ=255)
                                                                   gTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGAt                                                                                  <  3:814122/149‑1 (MQ=255)
                                                                   gTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGAt                                                                                  <  4:3199108/149‑1 (MQ=255)
                                                                    ttttACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATg                                                                                 >  1:3264266/1‑149 (MQ=255)
                                                                      ttACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGt                                                                               <  1:4956501/149‑1 (MQ=255)
                                                                       tACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGtt                                                                              <  1:7172292/149‑1 (MQ=255)
                                                                       tACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGtt                                                                              <  1:236994/149‑1 (MQ=255)
                                                                          ctacACGCTG‑CGCTGC‑GGAGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCACGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTcc                                                                           >  4:557684/4‑149 (MQ=255)
                                                                          cTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGAt                                                                                                                                                 >  2:3732397/1‑79 (MQ=255)
                                                                          cTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTcc                                                                           >  4:6893883/1‑149 (MQ=255)
                                                                          cTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTcc                                                                           <  2:2405103/149‑1 (MQ=255)
                                                                            tCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCcac                                                                         >  1:3948766/1‑149 (MQ=255)
                                                                            tCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCcac                                                                         <  1:6049319/149‑1 (MQ=255)
                                                                            ccACGCTGACGCTGCAGGTGATCGCC‑GCAGCG‑CAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCcac                                                                         >  3:170285/2‑149 (MQ=38)
                                                                             cACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTcc                                                                           >  3:547683/1‑146 (MQ=255)
                                                                             cACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCcaca                                                                        <  1:3835172/149‑1 (MQ=255)
                                                                             cACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCcaca                                                                        >  4:1435944/1‑149 (MQ=255)
                                                                                gctg‑cgctgc‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCCTCGGCAGACCACCTGTCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATcc                                                                     <  2:6164899/149‑1 (MQ=255)
                                                                                gctg‑cgctgc‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATcc                                                                     >  1:6321351/1‑149 (MQ=255)
                                                                                gctg‑cgctgc‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAAACGGTAGAGATGGTTTCCACATcc                                                                     >  1:3595492/1‑149 (MQ=255)
                                                                                       ctgc‑gGCGGTCACCTGCAGCGGCAGGGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGaa                                                               <  2:2116767/149‑1 (MQ=255)
                                                                                        tgc‑gGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAAc                                                              <  2:7193601/149‑1 (MQ=255)
                                                                                            ggcgATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGt                                                           >  3:3421541/1‑149 (MQ=255)
                                                                                               gATCACCTGCAGCGTCAGCGTGGTTTTACCGAAAGATTCCGGTCCGTAGATTTCGACGATAAGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCtt                                                        <  2:484690/149‑1 (MQ=255)
                                                                                                 tCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCa                                                      <  2:1220221/149‑1 (MQ=255)
                                                                                                  cACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCAc                                                     <  4:5132158/149‑1 (MQ=255)
                                                                                                   aCCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCAcc                                                    >  3:5022751/1‑149 (MQ=255)
                                                                                                   aCCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCAcc                                                    <  2:1668024/149‑1 (MQ=255)
                                                                                                     cTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCa                                                  >  1:788962/1‑149 (MQ=255)
                                                                                                      tGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAg                                                 <  3:3904090/149‑1 (MQ=255)
                                                                                                      tGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCACATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAg                                                 <  1:1301331/149‑1 (MQ=255)
                                                                                                      tGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGGCCGTAGATTTCAACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAg                                                 >  2:6564774/1‑149 (MQ=255)
                                                                                                             tcagcgtGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCatg                                          >  4:2155452/1‑149 (MQ=255)
                                                                                                              cagcgtGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCatga                                         <  4:2396111/149‑1 (MQ=255)
                                                                                                              cagcgtGGTTGTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCAGCCCCAGGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCCCatga                                         >  4:4633762/1‑149 (MQ=255)
                                                                                                                gcgtGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCatgatg                                       <  4:3689582/149‑1 (MQ=255)
                                                                                                                  gtGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGa                                     <  4:3326674/149‑1 (MQ=255)
                                                                                                                  gtGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGa                                     <  4:2927912/149‑1 (MQ=255)
                                                                                                                   tGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAg                                    <  4:6424681/149‑1 (MQ=255)
                                                                                                                     gTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCACACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGcc                                  <  3:4682501/149‑1 (MQ=255)
                                                                                                                       tttACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCtt                                <  4:5381412/149‑1 (MQ=255)
                                                                                                                         tACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTa                              <  4:1138722/149‑1 (MQ=255)
                                                                                                                          aCCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTAc                             >  2:2977383/1‑149 (MQ=255)
                                                                                                                            cGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCCGGCGCATGATGGAGCCTTTACCa                           >  4:1606813/1‑149 (MQ=255)
                                                                                                                            cGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCa                           >  3:5895080/1‑149 (MQ=255)
                                                                                                                            cGGAAGATTCCGGTCCGTAGATTTCAACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCa                           >  1:5579465/1‑149 (MQ=255)
                                                                                                                              gAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTa                                                                                       <  4:1960001/87‑1 (MQ=255)
                                                                                                                              gAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTcc                                                                           <  3:922325/99‑1 (MQ=255)
                                                                                                                              gAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCaaa                         >  3:2136657/1‑149 (MQ=255)
                                                                                                                              gAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCaaa                         >  1:907138/1‑149 (MQ=255)
                                                                                                                               aaGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAAt                        <  4:6267466/149‑1 (MQ=255)
                                                                                                                                aGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAAtt                       >  4:4215524/1‑149 (MQ=255)
                                                                                                                                    tCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGttt                   <  2:3350424/149‑1 (MQ=255)
                                                                                                                                    tCCGGTCCGTAGATTTCAACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGttt                   <  3:2954919/149‑1 (MQ=255)
                                                                                                                                     ccGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTtc                  <  3:679146/149‑1 (MQ=255)
                                                                                                                                      cGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAg                                                 <  1:2626130/117‑1 (MQ=255)
                                                                                                                                      cGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTtct                 >  2:3932179/1‑149 (MQ=255)
                                                                                                                                          ccGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAAt             <  1:5912934/149‑1 (MQ=255)
                                                                                                                                              aGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCt                                 >  2:5523062/1‑125 (MQ=255)
                                                                                                                                              aGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTgg         >  4:765859/1‑149 (MQ=255)
                                                                                                                                               gATTTCAACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGACAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGc        >  3:5350026/1‑149 (MQ=255)
                                                                                                                                                     gacgaTACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGtg  >  1:5160099/1‑149 (MQ=255)
                                                                                                                                                     gacgaTACGGCCCAGCAGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAACGCGCACCGGTAGAGATGCTTTCCACATCCATGGAACGGTCTTCAACCAGGCGCATGATCGACCCTTTGCCAAAGTGTTTCTCAATCTGGCCCAGtg  >  4:2572836/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTG‑CGCTGC‑GGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTG  >  NC_000913/2823434‑2823729

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑