Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,239,611 C→A 9.1% intergenic (+66/‑333) alx → / → sstT putative membrane‑bound redox modulator/sodium:serine/threonine symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,239,6110CA 9.1% 565.1 / 40.6 176intergenic (+66/‑333)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base C (100/56);  new base A (11/5);  total (111/65)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

Only 100 of 176 total aligned reads displayed.

GCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAT  >  NC_000913/3239463‑3239757
                                                                                                                                                     |                                                                                                                                                  
gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                   <  1:643968/150‑1 (MQ=255)
   tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                >  1:4113364/1‑150 (MQ=255)
   tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                >  2:1969550/1‑150 (MQ=255)
    ggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                                                                                               <  1:1511569/150‑1 (MQ=255)
    ggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                                                                                               <  1:3602144/150‑1 (MQ=255)
     gtgtTTGGCATTCTGGTGATGACGTTTATTATCATCGCCTGGGTGAATTTTCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTTCCCTGTCATTAACATGAGCAGGGTTTTCTCCCTCTCCCTTCCAGGgag                                                                                                                                              >  2:371624/1‑150 (MQ=255)
     gtgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                                                                                              >  2:1546605/1‑150 (MQ=255)
     gtgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                                                                                              >  2:2261390/1‑150 (MQ=255)
       gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                                                                                                            <  2:1895713/150‑1 (MQ=255)
        tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                           >  2:2113177/1‑150 (MQ=255)
        tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAgg                                                                                                                                           <  1:4124132/150‑1 (MQ=255)
         cgggCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAggg                                                                                                                                          <  2:3436450/148‑1 (MQ=255)
          tGGCATTCTGGTGATGACGTTTATTATCACCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAGTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCACCCTCTCCCTGGAAGGGAGAggga                                                                                                                                         <  1:3774242/150‑2 (MQ=255)
          tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGt                                                                                                                                         >  2:121051/1‑150 (MQ=255)
          tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAggga                                                                                                                                         <  1:830792/150‑2 (MQ=255)
            gCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCg                                                                                                                                       <  2:3444837/150‑1 (MQ=255)
              aTTCTGTTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGTGGGTCggg                                                                                                                                     <  1:3021702/150‑1 (MQ=255)
                  tGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                 <  2:3815899/150‑1 (MQ=255)
                  tGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                 >  2:73235/1‑150 (MQ=255)
                   ggTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGATTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAg                                                                                                                                >  1:770175/1‑150 (MQ=255)
                          aCGTTTTTTATAAACGCCTGGTTTAATTATCGTAATTATAAGCAGCGGGTTTTATTATTTTTAATCTGCATAAGCGGTTTACCCTGTAATTAAATTGAGCACGTTTTTCTCCCTCTCCCTTCCTGGTTGTGGGTCGGGGTGAGGGTAAtt                                                                                                                         <  1:1958815/150‑1 (MQ=255)
                          aCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtt                                                                                                                         <  2:3309684/150‑1 (MQ=255)
                            gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                       >  2:2722012/1‑150 (MQ=255)
                            gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCGTCCAGGGAGAGGGTCGGGGTGAGGGTCAtttt                                                                                                                       >  2:2155944/1‑150 (MQ=255)
                            gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAACTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                       >  2:873051/1‑150 (MQ=255)
                             tttATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttttt                                                                                                                      >  2:1266682/1‑150 (MQ=255)
                              ttattaTCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTCATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGGGTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTGTTc                                                                                                                     >  2:1873218/1‑150 (MQ=255)
                              ttattaTCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTc                                                                                                                     >  1:3196320/1‑150 (MQ=255)
                              ttattaTCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTc                                                                                                                     >  2:2497640/1‑150 (MQ=255)
                                    tCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCg                                                                                                               >  2:1712909/1‑150 (MQ=255)
                                      aaCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCgct                                                                                                             >  2:4019331/1‑148 (MQ=255)
                                      aaCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                                                                                                             >  2:2665688/1‑150 (MQ=255)
                                              gttGAATTATCGGCATGATACGCTGGGGGGTGTATATTTTTTAATCTGCCTAAGCGTGTTACCCTGTCATTAACATGAGCACCGTTTGCGCCCCCCCCCTTCCCGGGAGTGGGTCGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCt                                                                                                     <  2:3748242/148‑1 (MQ=255)
                                               gTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTg                                                                                                    >  2:3068005/1‑150 (MQ=255)
                                                 tatttATCGGCATGATAAGCAGGGGGGTGGATAATTTTTATTCTGCCTACGCCGTTTACCCTGACATTAACATGAGCACCGTTTACTCCCTCTCCCTTCCTGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGtt                                                                                                  <  2:2805931/147‑1 (MQ=255)
                                                   aTTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcc                                                                                                >  1:864250/1‑150 (MQ=255)
                                                         ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAc                                                                                          >  2:909879/1‑150 (MQ=255)
                                                          gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                         >  1:446169/1‑150 (MQ=255)
                                                          gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                         >  2:1424267/1‑150 (MQ=255)
                                                          gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                         >  1:2059250/1‑150 (MQ=255)
                                                          gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                         >  1:892367/1‑150 (MQ=255)
                                                              gATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTaa                                                                                     >  2:3102057/1‑150 (MQ=255)
                                                                tAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAcc                                                                                   >  1:15138/1‑150 (MQ=255)
                                                                     tggggggTGTATAATTTTTAATCTGCATAAGCGGTGTACCCTGTCATTAACATGAGCACCGTTTCCGCCCTCGCCCTCCCTGGGTGTGGGTGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                              <  1:3812118/147‑1 (MQ=37)
                                                                     agcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                              >  1:2435688/1‑150 (MQ=255)
                                                                      gggggtTGTATAATTTTTAATCTGCCTACGCCGTGTACCCTGTCTTTAACATGCGCACGGTTTTCTCCCACTCCCTCCCCGGTTGTGGGTCGGGTTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTcc                                                                             <  2:697736/148‑1 (MQ=25)
                                                                       cGGGGGGGATAATTTTTAATCGGCCTAAGCCGTGTACCCTCTCATTAACAGGCGCACCGTGTGCGCCCGCTCCCTGCCCGGGAGAGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTccc                                                                            <  2:4116353/150‑1 (MQ=25)
                                                                            tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGTAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAc                                                                       <  2:2598929/150‑1 (MQ=255)
                                                                            tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAGGCTGGCCGGTTCCCCTCACCCTAACCCTCTCCCCAAAc                                                                       >  2:1381226/1‑150 (MQ=255)
                                                                               tttttttttATTGTGCCTATGCGTTGTACCCTTTTTTTAATTTGAGCACGGTTTTCTCCCTCCCCGTTCGTGGTTTTGGGTGGGGGTTGGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACggg                                                                    <  2:2868672/146‑1 (MQ=16)
                                                                               atatTTTTTATTCTGCCTAGGCGTTTTACCCTGACTTTAACTTGAGCACCGTTTTCTCCCCCTCCCTGCCCGGGTGTGGGTGGGGGTGAGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACggg                                                                    <  2:3002899/150‑1 (MQ=21)
                                                                               aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACggg                                                                    >  1:4197632/1‑150 (MQ=255)
                                                                                tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgggg                                                                   >  1:2762269/1‑150 (MQ=255)
                                                                                 aaTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGc                                                                  >  1:2236416/1‑150 (MQ=255)
                                                                                    ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGGTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCgtg                                                               >  1:153886/1‑148 (MQ=255)
                                                                                    ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAg                                                               >  1:1178824/1‑150 (MQ=255)
                                                                                    ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAg                                                               >  1:208714/1‑150 (MQ=255)
                                                                                    ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAg                                                               >  1:2286425/1‑150 (MQ=255)
                                                                                    ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAg                                                               >  1:3252694/1‑150 (MQ=255)
                                                                                     tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAgg                                                              >  2:1318700/1‑150 (MQ=255)
                                                                                          tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACt                                                         >  1:105532/1‑150 (MQ=255)
                                                                                            tGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGATCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGa                                                       >  1:342935/1‑150 (MQ=255)
                                                                                              ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACCGGGCGAGGGGACTGAcc                                                     >  2:1618983/1‑150 (MQ=255)
                                                                                                tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                   >  1:1691463/1‑150 (MQ=255)
                                                                                                aacgcggtgtACCCTTAAATTAACTTGAGCACCTTTTTCTCCCCCCCCCTGCCGGGTTGTGGGTGGGGGTGTGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                   <  2:3218722/144‑1 (MQ=21)
                                                                                                      gtgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtt                                             >  1:3906461/1‑150 (MQ=255)
                                                                                                           aCCCTTACATTAAAATGCCCCCCGTTTCCCCCCGCCCCCTTCCGGGGTGTGGGGGGGTGTGCGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                         <  2:1441387/150‑1 (MQ=18)
                                                                                                           aCCCTGTCATTAACCTGAGCCCCGTGTCCTCCCGCGCCGGGCCCGGGAGCGGGTGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                         <  2:977728/150‑1 (MQ=21)
                                                                                                           aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                         <  2:3473044/150‑1 (MQ=255)
                                                                                                             ccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGCCTGACCGAGTCCTGTTTTGa                                       >  2:248358/1‑150 (MQ=255)
                                                                                                             ccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                                       >  1:1589035/1‑150 (MQ=255)
                                                                                                              cTGTCATTAACATGAGCACCGTATTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                      >  1:2496977/1‑150 (MQ=255)
                                                                                                               tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtg                                     >  1:4149982/1‑150 (MQ=255)
                                                                                                                   aTTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgttgt                                 >  2:1406736/1‑150 (MQ=255)
                                                                                                                     tAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtca                               <  2:3166272/150‑1 (MQ=255)
                                                                                                                     tAACATGAGCACCGTTTGCGCCCTCGCCCTTCCCGGTTGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtca                               <  1:564622/150‑1 (MQ=25)
                                                                                                                       aCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACCGGGCCAGGGGACTGACCCAGTCCTTTTTTGATGTTGtcata                             >  1:3687516/1‑149 (MQ=255)
                                                                                                                        cATGAGCACCTTTTTCCCCCCCTCCCTTCCTGGGTGTGGGTGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcatca                            <  1:344171/150‑1 (MQ=25)
                                                                                                                          tGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGt                          >  2:3361310/1‑150 (MQ=255)
                                                                                                                             gCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTCATTTTTCGCACCGATGCTGGCCGGTTCCCCTCACCCTAACCCTCTCCCCAGACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTg                       >  2:3466985/1‑150 (MQ=255)
                                                                                                                              cACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCACACCCTAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTTTCATAAGTCTgg                      >  1:601654/1‑150 (MQ=255)
                                                                                                                              cACCGTTTTCGCCCTCCCCCGTCCCGGGAGTGGGTGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTgg                      <  1:1319064/150‑1 (MQ=25)
                                                                                                                                cctttttCCCCCTCCCCCTCCCCGGGAGAGGGTGGGGTTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                    <  1:690177/147‑1 (MQ=25)
                                                                                                                                 cGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGTTGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCACAAACGGGGGGAGGGGACTGACCGAGTCGTTTTTTGTTGTTGTGATCAGTCTGGAAg                   >  1:2548158/1‑150 (MQ=25)
                                                                                                                                 cGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAg                   >  1:4053801/1‑150 (MQ=255)
                                                                                                                                 cGTTTCCCCCCGCCCCCTGCCCGGGAGTGGGTGGGGGTGGGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAg                   <  1:1266682/150‑1 (MQ=25)
                                                                                                                                  gTTTCCCCCCCCCCCCTACCTGGGTGTGGGTGGGGGTGTGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGc                  <  1:688046/150‑1 (MQ=25)
                                                                                                                                      tctccctctcccTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCACACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCa              <  2:1049425/150‑1 (MQ=255)
                                                                                                                                        tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGCCGGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACg            >  1:2963998/1‑150 (MQ=255)
                                                                                                                                        tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTGGGCATCAGTCTGGAAGCCGCACg            >  2:1095979/1‑150 (MQ=255)
                                                                                                                                          cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCCCCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTGTTGATGTTGTCCTCAGTCTGGAAGCCGCACGtt          >  2:2087345/1‑150 (MQ=25)
                                                                                                                                          cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt          >  2:2087338/1‑150 (MQ=255)
                                                                                                                                          cctctccctGGAAGGGAGAGGGTAGGGGTGAGGGTAATTTTTCGCACCGATGCGGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtg          >  1:3378907/1‑149 (MQ=37)
                                                                                                                                           ctctccctCCGGGGGTGTGGGGGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg         <  1:2001549/150‑1 (MQ=37)
                                                                                                                                            tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        >  2:2027765/1‑150 (MQ=255)
                                                                                                                                            tctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        >  2:541711/1‑150 (MQ=255)
                                                                                                                                            tctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        >  2:1628184/1‑150 (MQ=255)
                                                                                                                                            tctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        >  2:1628266/1‑150 (MQ=255)
                                                                                                                                             ctccctTCCTGGGTGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGc       <  1:1726977/150‑1 (MQ=255)
                                                                                                                                              tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt      >  2:2382900/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAT  >  NC_000913/3239463‑3239757

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: