Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,239,609 T→G 7.9% intergenic (+64/‑335) 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,6090TG 7.9% 515.8 / 15.3 172intergenic (+64/‑335)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base T (116/34);  new base G (5/8);  total (121/51)
Fisher's exact test for biased strand distribution p-value = 4.92e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

TGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATT  >  NC_000913/3239460‑3239758
                                                                                                                                                     |                                                                                                                                                     
tGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCtt                                                                                                                                                       >  1:224161/1‑150 (MQ=255)
  ggCGTGGTGTTTGGCAATCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCACTCCCTTcc                                                                                                                                                     >  1:3849317/1‑150 (MQ=255)
   gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTTTCGGCATGATAAGCAGCGGGGTGGATAATTTTTATTCTGGCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    >  1:165126/1‑150 (MQ=255)
   gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  1:22083/150‑1 (MQ=255)
     gtggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                                                                                                                  <  1:1853161/150‑1 (MQ=255)
     gtggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                                                                                                                  <  1:2579088/150‑1 (MQ=255)
     gtggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                                                                                                                  >  1:4000494/1‑150 (MQ=255)
     gtggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                                                                                                                  >  2:240995/1‑150 (MQ=255)
     gtggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                                                                                                                  <  2:3212641/150‑1 (MQ=255)
     gtggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                                                                                                                  >  2:814807/1‑150 (MQ=255)
     gtggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCCCCCTCTCCCTTCCAgg                                                                                                                                                  <  2:2832881/150‑1 (MQ=255)
       ggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                                                                                                >  1:1642580/1‑150 (MQ=255)
       ggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                                                                                                <  2:3475307/150‑1 (MQ=255)
        gtgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                                                                                               <  1:510457/150‑1 (MQ=255)
        gtgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                                                                                               <  2:899189/150‑1 (MQ=255)
          gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGagag                                                                                                                                             <  2:3062233/150‑1 (MQ=255)
           tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  1:1023187/1‑150 (MQ=255)
             tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGt                                                                                                                                          >  2:3148860/1‑150 (MQ=255)
             tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGt                                                                                                                                          >  1:1029100/1‑150 (MQ=255)
               gCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCg                                                                                                                                        >  2:625302/1‑150 (MQ=255)
                        tgatgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAggg                                                                                                                               >  2:1966748/1‑150 (MQ=255)
                           tgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCCTTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTaa                                                                                                                            >  1:718149/1‑150 (MQ=255)
                            gaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAt                                                                                                                           >  2:436561/1‑150 (MQ=255)
                             gcGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCCTTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtt                                                                                                                          >  1:3126254/2‑150 (MQ=255)
                                  tattatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCg                                                                                                                     >  1:2283889/1‑150 (MQ=255)
                                  tattatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCg                                                                                                                     >  2:2799100/1‑150 (MQ=255)
                                  tattatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCGTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCg                                                                                                                     >  2:2961800/1‑150 (MQ=255)
                                          aCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATg                                                                                                             >  2:309682/1‑150 (MQ=255)
                                            gCCTGGGTGAATTATCGGCATGATAAGCTGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACTTGAGCACCGTGTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGTTGAGGGTAATTTTTCGCACCGATGCt                                                                                                           <  2:603188/150‑1 (MQ=255)
                                             ccTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTg                                                                                                          >  1:1918048/1‑150 (MQ=255)
                                              cTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTgg                                                                                                         >  2:1925760/1‑150 (MQ=255)
                                              cTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTgg                                                                                                         >  2:2513218/1‑150 (MQ=255)
                                              cTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTgg                                                                                                         >  2:2513276/1‑150 (MQ=255)
                                                gggTGAATTATCGGCATGATACGCCGGGGGGTGTATAATTTTTAAGCTGCCTACGCGGTGTACCCGGTCATTAACATGCGCCCCGTTTTCCCCCACGCCCTCCCCGGGAGAGGGTGGGGGGGTGGGTAATTTTTCGCACCGATGCTGGcc                                                                                                       <  2:3973466/150‑1 (MQ=255)
                                                gggTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGcc                                                                                                       <  2:1010087/150‑1 (MQ=255)
                                                        tATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTACTTTTTCGCACCGATGCTGGCCTGTTcccc                                                                                               >  2:1377611/1‑150 (MQ=255)
                                                          tCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTc                                                                                             >  2:2264821/1‑150 (MQ=255)
                                                             gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                          >  2:1966315/1‑150 (MQ=255)
                                                              cATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAccc                                                                                         >  2:1885410/1‑150 (MQ=255)
                                                                tGATAAGCAGCGGGGTGGATAATTATTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTa                                                                                       >  2:3147265/1‑150 (MQ=255)
                                                                  aTAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGGTGGCCTGTTCCCCTCACCCTAAc                                                                                     >  1:880501/1‑150 (MQ=255)
                                                                  aTAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAc                                                                                     >  1:3160957/1‑150 (MQ=255)
                                                                  aTAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAc                                                                                     >  2:4012658/1‑150 (MQ=255)
                                                                            gggTGGATAATTTTTAATCTGCCTAAGCCGAGTACCCTGTCAATAACATTAGCACAGTTTTCTACCACTACCATCCAGGGAGATGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTTGCCTGTTCCCCACAACCAAACCCACACCCCa                                                                           >  1:702396/1‑150 (MQ=255)
                                                                                 gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgg                                                                      >  1:1550107/1‑150 (MQ=255)
                                                                                        tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAgg                                                               >  2:3400539/1‑150 (MQ=255)
                                                                                           aaTCTGCCTACGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCGCCCTCCCCCTTCCCGGGTGTGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGa                                                            <  2:4016486/150‑1 (MQ=37)
                                                                                            aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAggggcc                                                           >  1:1882334/1‑148 (MQ=255)
                                                                                             tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACt                                                          >  2:2707104/1‑150 (MQ=255)
                                                                                              cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTg                                                         >  1:1736454/1‑150 (MQ=255)
                                                                                                 ccTACGCGGTGGACCCTGTCATTAACATGAGCCCCGTTTTCGCCCTCCCCCTGCCCGGGTGCGGGGGGGGGGGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGAcc                                                      <  2:1918033/150‑1 (MQ=21)
                                                                                                 ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGAcc                                                      >  2:3223846/1‑150 (MQ=255)
                                                                                                   tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                    >  2:2732163/1‑150 (MQ=255)
                                                                                                    aaGCCGTGTACCCTGTCATTAACATGTGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGTGGGTAATTTTTCGCACCGTTGCTGGCCTGTTCCCCTCAGCCTAACCCTCTCCCCAAACGGGGCGAGGGGAGTGTCCGAg                                                   >  2:1851176/1‑150 (MQ=255)
                                                                                                    aaGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAg                                                   >  2:3896464/1‑150 (MQ=255)
                                                                                                    aaGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCACACCCAAACCCTCACCCCAAACGGGGCGAGGGGCCTGACCGAg                                                   >  2:1850958/1‑150 (MQ=255)
                                                                                                     aGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGt                                                  <  2:556973/150‑1 (MQ=255)
                                                                                                           gtACCCTTTAATTAACATGACCACGTTTTACCCCCACCCCCTCCCGGGGTGTGGGTGGGGGTGTGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttt                                            <  1:436561/150‑1 (MQ=21)
                                                                                                            tACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCttttt                                           >  1:2488895/1‑150 (MQ=255)
                                                                                                             ccccTGTAATTAACATGCGCACCTTTTTCTCCCTCCCCCTTCCGGGGTGCGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  1:4107323/149‑1 (MQ=25)
                                                                                                             aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCTGGGTGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  1:3098732/150‑1 (MQ=255)
                                                                                                             aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  1:3939326/150‑1 (MQ=255)
                                                                                                             aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          >  1:3367403/1‑150 (MQ=255)
                                                                                                               ccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTGTTTTGa                                        >  1:2211198/1‑150 (MQ=255)
                                                                                                                 tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtg                                      >  2:1741024/1‑150 (MQ=255)
                                                                                                                 tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtg                                      >  2:1740552/1‑150 (MQ=255)
                                                                                                                  gTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgt                                     >  1:238603/1‑150 (MQ=255)
                                                                                                                       tAACATGAGCACCGTGTTCGCCCTCTCCCGTCCAGGGAGAGGGTCGGGGTGAGGGTCATTTTTCGCACCGAGGCTGGCCTGTTCCCCTCACCCTGCGCCTCTCCCCCCCCGGGGCGAGGGGCCTGCCCGAGTCCTGTTGTGAGGTTGtcg                                >  2:889348/1‑149 (MQ=255)
                                                                                                                          cATGCGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcatca                             <  2:3425443/150‑1 (MQ=255)
                                                                                                                          cATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcatca                             <  1:3014993/150‑1 (MQ=255)
                                                                                                                          cATGAGCACCGTTTCCTCCCTCTCCCCTCCCGGGAGAGGGTCGGGGGGCGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcatca                             <  1:3363343/150‑1 (MQ=37)
                                                                                                                           aTGAGCACCGTTTTCTCCCTCTCCCTGCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGCCTGGCCGAGTCCTGTTTTGTTGTTGGCATCAg                            >  2:1856514/1‑150 (MQ=25)
                                                                                                                              aGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAACGGGGCGAGGGGCCTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCt                         >  1:2199308/1‑150 (MQ=255)
                                                                                                                              aGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCACACCCTAACCCACACCCCAAACGGGGCGAGGGGACGGACCGAGTCCGTTTTTGATGTTGTCATAAGTCt                         >  1:6316/1‑150 (MQ=255)
                                                                                                                              aGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCTCACCCAAACCCACACCCCAAACCGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCt                         >  1:3462851/1‑150 (MQ=255)
                                                                                                                               gccgcGTTTTCTCCCTCGCCCTGCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTg                        <  2:3740183/146‑1 (MQ=255)
                                                                                                                               gccccGTTGTCGCCCTCGCCCTGCCCGGGTGAGGGTGGGGGGGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTg                        <  2:2316665/147‑1 (MQ=25)
                                                                                                                               gCACCGTTTCCGCCCGCACCCGCCGCGGGTGTGGGTGGGGGGGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTg                        <  2:3853620/150‑1 (MQ=25)
                                                                                                                                 aCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCCAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTTGTTGTCATCAGTCTGGa                      >  1:950239/1‑150 (MQ=255)
                                                                                                                                  ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                     >  1:2861177/1‑150 (MQ=255)
                                                                                                                                  ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                     >  1:2772868/1‑150 (MQ=255)
                                                                                                                                  ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                     >  1:1029481/1‑150 (MQ=255)
                                                                                                                                  ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCTCACCCTAACCCTCACCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTTATTTTGTCATCAGGCTGGaa                     >  1:2039419/1‑150 (MQ=255)
                                                                                                                                  ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCGGTTCCCCTCACCCTAACCCACTCCCCAAACGGGGCGAGGGGGCTGACCGAGTCCTTGTTTGATGTTGTCATCAGTCTGGaa                     >  1:1265771/1‑150 (MQ=255)
                                                                                                                                          tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGACCTGTTTTGAAGTTGTCATCAGTCTGGAAGCCGCACg             >  2:1114483/1‑150 (MQ=255)
                                                                                                                                          tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCGCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACg             >  2:1446859/1‑150 (MQ=255)
                                                                                                                                           ccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCTCCCTTACCCTCTCCCCAAACGGGGCGTGGGGCCTGACCGAGTCCTTTTTTGCTGTTGTCATCAGGGGGGGAGCcgcgcgg            >  2:2456717/1‑149 (MQ=255)
                                                                                                                                           ccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTAGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCGCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGGCATCAGTCTGGAAGCCGCACGt            >  1:2091335/1‑150 (MQ=255)
                                                                                                                                           ccctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCACACCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTTGTTGTCATAAGTCTGGAAGCCGCACGt            >  1:1467123/1‑150 (MQ=255)
                                                                                                                                            cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           >  2:2295139/1‑150 (MQ=255)
                                                                                                                                             ctctccctTCCCGGGAGAGGGTGGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg          <  2:4079598/150‑1 (MQ=255)
                                                                                                                                             ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg          >  2:3261725/1‑150 (MQ=255)
                                                                                                                                             ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg          >  2:3330349/1‑150 (MQ=255)
                                                                                                                                                tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGcg       >  2:1242419/1‑149 (MQ=255)
                                                                                                                                                tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCTCACCCTAACCCTATCCCCAAACGGGGCGAGGGGGCGGACCGAGTCCTTTTTTGATGTTGTCATCAGTCGGGAATCCGCACGATGGCt       >  1:927861/1‑150 (MQ=255)
                                                                                                                                                  cctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCttt     >  2:2905745/1‑150 (MQ=255)
                                                                                                                                                  cctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAGGTTGTCCTCAGTCTGGAAGCCGCACGTTGGCgtt     >  2:3007230/1‑150 (MQ=255)
                                                                                                                                                   ctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTa    >  1:3222616/1‑150 (MQ=255)
                                                                                                                                                   ctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTa    >  1:3885925/1‑150 (MQ=255)
                                                                                                                                                    ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTTGTTGTCATCAGTCTGGAAGCCGCACGTTGGCGTTAt   >  1:46335/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATT  >  NC_000913/3239460‑3239758

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: