Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,239,609 T→G 6.4% 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 6.4% 572.9 / 16.9 191intergenic (+64/‑335)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base T (121/49);  new base G (3/10);  total (124/67)
Fisher's exact test for biased strand distribution p-value = 8.12e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

TGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAT  >  NC_000913/3239460‑3239757
                                                                                                                                                      |                                                                                                                                                    
    cGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                                                                                                  <  1:3945031/150‑1 (MQ=255)
    cGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                                                                                                  <  1:1936079/150‑1 (MQ=255)
     gtggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                                                                                                                 <  1:3269369/150‑1 (MQ=255)
       ggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                                                                                               <  2:1775926/150‑1 (MQ=255)
        gtgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTGGAAGGgag                                                                                                                                              <  1:3141187/150‑1 (MQ=255)
        gtgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTGGAAGGgag                                                                                                                                              <  1:1875969/150‑1 (MQ=255)
         tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                             <  2:520702/150‑1 (MQ=255)
         tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                             <  2:3809233/150‑1 (MQ=255)
         tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                             >  1:1974873/1‑150 (MQ=255)
           tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                           <  1:1118251/150‑1 (MQ=255)
           tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                           >  1:3726733/1‑150 (MQ=255)
           tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAgg                                                                                                                                           <  2:3781738/150‑1 (MQ=255)
            ttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAggg                                                                                                                                          <  1:2753134/150‑1 (MQ=255)
             tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGt                                                                                                                                         <  1:257100/150‑1 (MQ=255)
             tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGt                                                                                                                                         <  2:2101896/150‑1 (MQ=255)
             tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCACCCTCTCCCTGGAAGGGAGAggga                                                                                                                                         <  2:3892425/150‑2 (MQ=255)
              ggCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                                        >  2:2681976/1‑150 (MQ=255)
               gCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCg                                                                                                                                       <  2:779926/150‑1 (MQ=255)
                cATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                                      >  1:1374209/1‑150 (MQ=255)
                cATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                                      <  1:3982543/150‑1 (MQ=255)
                  ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                                                                                                    >  2:4036587/1‑150 (MQ=255)
                   tCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAGCA‑TGAGCACCGTTTTCTCCCTCGCCCTTCCAGGGAGAGGGTCGGGGt                                                                                                                                   <  2:4111193/150‑1 (MQ=255)
                   tCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGt                                                                                                                                   <  2:229948/150‑1 (MQ=255)
                     tGGTGATGCCGTTTATTATCAACGCCTGGGTGAATGATCGGCATGATACGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGa                                                                                                                                 <  2:1468035/150‑1 (MQ=255)
                     tGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                 >  1:1548752/1‑150 (MQ=255)
                     tGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                 <  2:2148339/150‑1 (MQ=255)
                      ggTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAg                                                                                                                                >  2:3536013/1‑150 (MQ=255)
                            gaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATACGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGATTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAt                                                                                                                          >  2:3223042/1‑150 (MQ=255)
                             aCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAACTTTTAATCTGCCTAAGCCGTGTACCCTGTCACTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtt                                                                                                                         >  1:98312/1‑150 (MQ=255)
                               gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                       >  1:2385688/1‑150 (MQ=255)
                               gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                       >  1:1693782/1‑150 (MQ=255)
                               gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                       >  2:3488060/1‑150 (MQ=255)
                               gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                       >  2:3407523/1‑150 (MQ=255)
                               gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                       >  2:2816703/1‑150 (MQ=255)
                                tttATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttttt                                                                                                                      >  1:3835203/1‑150 (MQ=255)
                                tttATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttttt                                                                                                                      >  1:3968465/1‑150 (MQ=255)
                                       tCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCg                                                                                                               >  1:2555215/1‑150 (MQ=255)
                                       tCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCg                                                                                                               >  2:1724074/1‑150 (MQ=255)
                                        cAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGa                                                                                                              >  2:1851328/1‑150 (MQ=255)
                                         aaCGCCTGGGTGAATTATCGGCATGATCAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                                                                                                             >  2:1356585/1‑150 (MQ=255)
                                         aaCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                                                                                                             >  1:4041601/1‑150 (MQ=255)
                                               tGGGTGAATTATCGGTATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGc                                                                                                       >  2:2702143/1‑150 (MQ=255)
                                               tGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATTCTGGc                                                                                                       >  2:3992068/1‑150 (MQ=255)
                                                            ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAc                                                                                          >  2:3092609/1‑150 (MQ=255)
                                                            ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAc                                                                                          >  1:1615002/1‑150 (MQ=255)
                                                             gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                         >  1:228426/1‑150 (MQ=255)
                                                                 gATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTTTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGTGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTaa                                                                                     <  2:1019373/150‑1 (MQ=255)
                                                                 gATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTaa                                                                                     >  2:2000093/1‑150 (MQ=255)
                                                                   tAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAcc                                                                                   >  1:2895255/1‑150 (MQ=255)
                                                                   tAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTGTCTCCCTCCCCCTCCCTGGGATAGGGTCGGGGTTTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAcc                                                                                   <  1:3031629/150‑1 (MQ=39)
                                                                                 gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCACCCCAAACgg                                                                     >  2:1438801/1‑150 (MQ=255)
                                                                                   tattttttAATCTGCCTACGCCGTGTACCCTGAAATTAACA‑TGAGCACCTTTTTCTCCCCCCCCCTACCCGGGCGTGGGTGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgggg                                                                   <  2:3943990/147‑1 (MQ=25)
                                                                                   aaaTTTTTAATCTCCATAAGCCGTTAACCCTGTAATTAACA‑TTACCACCTTTTTCCCCCTCACCCTACCCGGTAGAGGCTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgggg                                                                   <  1:2030272/149‑1 (MQ=21)
                                                                                       ttttAATCTGCCTAAGCCGTGTACCCTGTCTTTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGTGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAg                                                               <  2:2058799/150‑1 (MQ=255)
                                                                                        tttAATCTGCCTAGGCGGTGCACCCTTTAATTAACT‑TGCGCCCCGTTTTCCCCCGCCCCCTGCGCGGGTGTGGGGGGGGTTGTGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAgg                                                              <  1:2439267/150‑1 (MQ=16)
                                                                                        tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAgg                                                              >  1:1189117/1‑150 (MQ=255)
                                                                                         ttAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAggg                                                             >  2:2282382/1‑150 (MQ=255)
                                                                                         ttAATCTCCATAAGCGTTTAACCCTGAC‑TTAACATTGAGCACGGTTTCCCCCCTCCCCCTCCCGGGTTGTGGGTGGGGTTTTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAggg                                                             <  1:3663114/150‑1 (MQ=12)
                                                                                         ttAATCGGCCTACGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCGCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAggg                                                             <  2:2421012/150‑1 (MQ=255)
                                                                                           aaTCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGGGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGCCGAGGGGa                                                           <  2:671206/150‑1 (MQ=255)
                                                                                            aTCCGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGAc                                                          >  2:4006678/1‑150 (MQ=255)
                                                                                             tCTGCCTAAGCCGTGTACCCTGTCATTATCA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACt                                                         >  1:3507548/1‑150 (MQ=255)
                                                                                              cTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTg                                                        >  1:1411323/1‑150 (MQ=255)
                                                                                                  cTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCg                                                    >  1:1990000/1‑150 (MQ=255)
                                                                                                   tAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCACACCCCAAACGGGGCGAGGGGACTGACCGa                                                   >  1:1817218/1‑150 (MQ=255)
                                                                                                    acgCCGTGAACCCTTTCATTAACT‑TGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAg                                                  <  1:499476/148‑1 (MQ=255)
                                                                                                    aaGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAg                                                  >  2:2292540/1‑150 (MQ=255)
                                                                                                          tgtACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCttg                                            >  2:3558619/1‑149 (MQ=255)
                                                                                                             ccccTTTAATTAACA‑TGTGCACGGTTTTCTCCCTCTCCGTTCCCGGGTTTGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                         <  2:2554253/149‑1 (MQ=25)
                                                                                                               ccTTTCATTAACT‑TGTGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                                       <  2:4101139/150‑1 (MQ=25)
                                                                                                               ccTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGGTGGCCTGTTCCCCTCACCCTAACCCACTCCCCGAACGGGGCGATGGGACTGACCGAGTCCTTTTTTGa                                       >  2:308324/1‑150 (MQ=255)
                                                                                                               ccTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTGTTTTgt                                       >  1:1579869/1‑149 (MQ=255)
                                                                                                                cTGTCATTAACA‑TGAGCACCGTTTACTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                      <  2:3731424/150‑1 (MQ=255)
                                                                                                                cTGACATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGTGGGTCGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                      <  2:3457445/150‑1 (MQ=255)
                                                                                                                 tatCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCATGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtg                                     >  1:2673210/3‑150 (MQ=255)
                                                                                                                     aTTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgttgt                                 >  1:3360074/1‑150 (MQ=255)
                                                                                                                      ttAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtc                                <  1:4076485/150‑1 (MQ=255)
                                                                                                                      ttAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAATTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGACTCCTTTTTTGATGTTGtc                                >  2:3234815/1‑150 (MQ=255)
                                                                                                                      ttAACA‑TGAGCACCGTTTTCGCCCACTCCCTTCCCGGGAGAGGGTGGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtc                                <  2:3400795/150‑1 (MQ=255)
                                                                                                                       tAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtca                               >  1:2941131/1‑150 (MQ=255)
                                                                                                                           a‑TGAGCACCGTTTTCTCCCTCCCCCTCCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAg                           <  1:1033436/150‑1 (MQ=255)
                                                                                                                               aGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCt                        >  2:1661651/1‑150 (MQ=255)
                                                                                                                                gCACCGTTTTCTCCCTCTCCCTCCCTGGGTTTGGTTCGGTTTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTg                       <  1:3443767/150‑1 (MQ=37)
                                                                                                                                gCACCGTTTTCTCCCTCTCCCGGCCCGGGAGAGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTg                       <  2:249765/150‑1 (MQ=255)
                                                                                                                                          ctccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCACTCCCCCAACGGGGCGAGGGGACTGACCGCGTCCATTTTTTGTGTTGTCATCAGTCTGGAAGCCGCAc             >  1:711009/1‑150 (MQ=255)
                                                                                                                                           tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACg            >  1:2101743/1‑150 (MQ=255)
                                                                                                                                           tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACg            >  1:2662560/1‑150 (MQ=255)
                                                                                                                                           tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACg            >  1:3605497/1‑150 (MQ=255)
                                                                                                                                            ccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCGAACCCTCTCCCCACACGGGGCGAGGGGACTGACCGAGTCCGTTTTTGATGGTGTCATCAGTCTGGAAGCCGCACGt           >  2:1180635/1‑150 (MQ=255)
                                                                                                                                             cctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt          >  2:366870/1‑150 (MQ=255)
                                                                                                                                             cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTTATGTTGTCATCAGTCAGGAAGCCGCACGtt          >  2:419829/1‑150 (MQ=37)
                                                                                                                                              ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg         >  1:3111199/1‑150 (MQ=255)
                                                                                                                                              ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg         >  1:3115128/1‑150 (MQ=255)
                                                                                                                                              ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCACACCCTAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGAGTACTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg         >  1:919608/1‑150 (MQ=255)
                                                                                                                                                ctccctTCCCGGGAGTGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGc       <  2:347755/150‑1 (MQ=255)
                                                                                                                                                ctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCACCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGc       >  2:1217016/1‑150 (MQ=255)
                                                                                                                                                ctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCGCACCCAAACCCGCTCCCCCAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGGTGGc       >  1:1072649/1‑150 (MQ=255)
                                                                                                                                                 tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt      >  2:3124422/1‑150 (MQ=255)
                                                                                                                                                   cctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCttt    >  2:929153/1‑150 (MQ=255)
                                                                                                                                                     ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACAGATGCTGGCCTGTTCCCCTCACCCTAACCCTCACCCCAAACGGGGCGAGGGGACTGACAGAGTCCTATTTTGATGTTGTCAACAGTATGGAAGCCGAACGTTGGCTTTAt  >  2:1263634/1‑150 (MQ=37)
                                                                                                                                                      |                                                                                                                                                    
TGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACA‑TGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAT  >  NC_000913/3239460‑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: ‑