Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,239,609 T→G 11.2% 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,6090TG11.2% 695.6 / 42.9 269intergenic (+64/‑335)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base T (158/68);  new base G (9/19);  total (167/102)
Fisher's exact test for biased strand distribution p-value = 1.89e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

TGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATT  >  NC_000913/3239460‑3239758
                                                                                                                                                      |                                                                                                                                                     
  ggCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTcc                                                                                                                                                     >  1:6392172/1‑150 (MQ=255)
   gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAATCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCAGTCAATAACATGAGCACCGTTTTCTCCCACTCCCTTCCa                                                                                                                                                    >  1:1239958/1‑150 (MQ=255)
   gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  2:4271405/150‑1 (MQ=255)
    cGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                                                                                                   >  1:205471/1‑150 (MQ=255)
      tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                 <  1:725043/150‑1 (MQ=255)
       ggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                                                                                                <  2:1498576/150‑1 (MQ=255)
         tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                              <  2:4184662/150‑1 (MQ=255)
         tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                              <  2:1857999/150‑1 (MQ=255)
         tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGgaga                                                                                                                                              <  1:1677672/150‑1 (MQ=255)
           tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAgg                                                                                                                                            <  1:5202062/150‑1 (MQ=255)
             tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGt                                                                                                                                          >  2:6205104/1‑150 (MQ=255)
             tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGt                                                                                                                                          >  1:5345561/1‑150 (MQ=255)
              ggCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCACCCTCTCCCTGGAAGGGAGAgggag                                                                                                                                         <  1:3459032/150‑3 (MQ=255)
                aaTTCTGGTTATGACGTTTATTATCAACGCCTGGTTGATTTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCATAAGCCGTGTACCCTGTCATTAACATTAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                                       <  1:4039477/149‑1 (MQ=255)
                 cttCTGGTGATGACGTTTATTATCAACGCCTGGGTGACTTATCGGCAGGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCGGTCATTAACCTGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCggg                                                                                                                                      <  2:5712501/149‑1 (MQ=255)
                  ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                                                                                                     >  2:2192472/1‑150 (MQ=255)
                   tCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGt                                                                                                                                    >  1:856886/1‑150 (MQ=255)
                    cTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                                                                                               >  1:3264866/1‑138 (MQ=255)
                    cTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                                                                                               <  2:3264866/138‑1 (MQ=255)
                     tGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                  <  2:4646109/150‑1 (MQ=255)
                       gTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAgg                                                                                                                                <  2:1459127/150‑1 (MQ=255)
                       gTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAgg                                                                                                                                <  2:6574042/150‑1 (MQ=255)
                         gatgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGt                                                                                                                              >  2:5070669/1‑150 (MQ=255)
                               gcttattaTCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                        >  2:940248/3‑150 (MQ=255)
                               gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                        >  2:1138438/1‑150 (MQ=255)
                               gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                        >  2:2855393/1‑150 (MQ=255)
                               gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                        >  2:2893506/1‑150 (MQ=255)
                                tttATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttttt                                                                                                                       >  2:2291075/1‑150 (MQ=255)
                                tttATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttttt                                                                                                                       >  1:3029857/1‑150 (MQ=255)
                                tttATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttttt                                                                                                                       >  2:541859/1‑150 (MQ=255)
                                 tttttatCAACGCGTGGGTGATTTTTCGGCATTATAAGCTGGGGGGTGTATATTTTTTAATCTGCCTAAGCGGTGTACCCTTTCATTAACATGCGCACGTTTTTCTCCCTCTCCCTTCCCGGTTTTGGGTGGGGGTGAGGGTAATTTTTc                                                                                                                      <  2:4968303/147‑1 (MQ=255)
                                 ttattaACAACGCTGGGGTTATTTATGGGCTTGATAAGCGGGGGGTTGGTTATTTTTTATTCTGCATACGGCGTTTACCCTGTCATTAACATGAGCACGTTTTTCACCCTCCCCCTGCCCGGGTTTGGGTGGGGGTTTGGGTAATTTTTc                                                                                                                      <  2:4845914/150‑1 (MQ=255)
                                  tattatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCg                                                                                                                     >  1:3080311/1‑150 (MQ=255)
                                   attatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGc                                                                                                                    >  1:1622370/1‑150 (MQ=255)
                                    ttatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCa                                                                                                                   >  2:4901244/1‑150 (MQ=255)
                                      atCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAcc                                                                                                                 >  1:2186255/1‑150 (MQ=255)
                                      atCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAcc                                                                                                                 >  2:1080416/1‑150 (MQ=255)
                                            gCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCt                                                                                                           >  2:1595945/1‑150 (MQ=255)
                                             ccTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTg                                                                                                          >  2:1413841/1‑150 (MQ=255)
                                                       ttATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTccc                                                                                                >  1:4193695/1‑150 (MQ=255)
                                                         aTCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCt                                                                                              >  1:286855/1‑150 (MQ=255)
                                                            ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAc                                                                                           >  1:6684319/1‑150 (MQ=255)
                                                             gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                          >  2:1747802/1‑150 (MQ=255)
                                                              cATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAccc                                                                                         >  2:3096925/1‑150 (MQ=255)
                                                                 taTAAGCCGGGGGGTTTATAATTTTTAATCTGCAAAAGCCGTGTACCCTCTCATTAACTTGAGCACCTTTTTCACCCTCCCCCTCCCCGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTaa                                                                                      <  2:6358597/149‑1 (MQ=25)
                                                                   tAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCCCCCTAAcc                                                                                    >  1:4776764/1‑150 (MQ=255)
                                                                    aaGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAccc                                                                                   >  2:5110778/1‑150 (MQ=255)
                                                                       cagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCACCCACTCCCATCAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAAcccaca                                                                                >  1:2476393/1‑147 (MQ=255)
                                                                            gggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCa                                                                           >  1:2563473/1‑150 (MQ=255)
                                                                                ggATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACg                                                                       >  2:3457925/1‑150 (MQ=255)
                                                                                ggATAATTTTTAATCGGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCGCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACg                                                                       <  2:4388250/150‑1 (MQ=255)
                                                                                 tttattttttAATCTGCCTACGCGGTTTACCTTGTAATTAACTTGCGCAGGTTTTTCCCCCTCCCCCTGCCGGGGTGTGGGTGGGGTTGTGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgg                                                                      <  2:254276/145‑1 (MQ=16)
                                                                                 gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGATTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCGCTCCCCAAACgg                                                                      >  2:5885855/1‑150 (MQ=255)
                                                                                   tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgggg                                                                    >  2:560922/1‑150 (MQ=255)
                                                                                      tttttATTGTGCCTACGCGGTGTCCCTTTCAATTAACTGGCGCACTTTTTTCCCCCTCGCCCTCCCGGGTTGTGGGTGGGGGTTTGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGa                                                                 <  2:5253124/150‑1 (MQ=12)
                                                                                      tttttAATCTGCATAAGCCGTGTACCCTTTCATTAACATGAGCACCTTTTTCTCCCTCTCCCTTCCTGGGTGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGa                                                                 <  1:6149053/150‑1 (MQ=37)
                                                                                       ttttAATCTCCCTCAGCCTGGTACCCTTTCATTAACCTGCGCCCGTTTTTCGCCCCCGCCCTGCCGGGGTGTGGGTGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAg                                                                <  1:1007463/150‑1 (MQ=18)
                                                                                        tttATTCTGCCTACGCCTTTTACCCTGTAATAAACTTGAGCACGTTTTTCGCCCACGCCCTTCCTGGGTGTGGGTGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAgg                                                               <  2:29331/150‑1 (MQ=18)
                                                                                         ttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAggg                                                              <  2:6632318/150‑1 (MQ=255)
                                                                                                   tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                    >  2:4525306/1‑150 (MQ=255)
                                                                                                     aGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGt                                                  >  2:6479915/1‑150 (MQ=255)
                                                                                                          tgtACCCTGTCATTAACTGGCGCCCCGTTTCCCCCCTCCCCCTGCCTGGTTGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCttt                                             <  2:809952/150‑1 (MQ=21)
                                                                                                          tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCACACGGGGCGAGGGGACTGACCGAGTCCttt                                             >  1:461309/1‑150 (MQ=255)
                                                                                                          tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCttt                                             >  1:2745003/1‑150 (MQ=255)
                                                                                                             ccccggccATTAACCGGCGCACCGTTTTCGCCCGCCCCCTGCGGGGGAGTGGGGGGGGGTGTGGGCAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  2:4734489/143‑1 (MQ=21)
                                                                                                             ccccTGACTTTAACCTCCGCACGGTTTCCCCCCCCCCCCGCCCGGGGCGTGGGGGGGGGTGTGGGCAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  1:5221661/149‑1 (MQ=16)
                                                                                                             aCCCTGTAATTAACTTGCGCCCGGTTTTCGCCCGCCCCCTGCCGGGTGGTGGGGGGGGGTGTGGGACATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  1:2613729/150‑1 (MQ=21)
                                                                                                             aCCCGGCCATTAACATGAGCACCGGTGTCTCCCTCTCCCGTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTGTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  2:748364/150‑1 (MQ=25)
                                                                                                                cTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                       <  2:4442072/150‑1 (MQ=255)
                                                                                                                 tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCCCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtg                                      >  1:5010081/1‑150 (MQ=255)
                                                                                                                  gTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCATTTTTgtttt                                     >  2:6272866/1‑146 (MQ=255)
                                                                                                                   tCATTACCTTGTGCCCGGTTTTCTCCCTCACGATTCCCGGGTGTGGGTGGGGGTTTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgtt                                    <  1:1061539/150‑1 (MQ=21)
                                                                                                                    cATTAACCTGAGCACCGTTTTCTCCCACCCCCCCCCGGGGTGTGGGTGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgttg                                   <  1:5923604/150‑1 (MQ=25)
                                                                                                                       taaaCATGAGCACGTTTTTCCCCCTCCCCCTTCCCGGGTGTGGGTGGGGGTGTGGGCAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCTAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtc                                 <  1:490773/148‑1 (MQ=21)
                                                                                                                        tAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcc                                >  2:102500/1‑149 (MQ=255)
                                                                                                                          aCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCCCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcatc                              >  2:1188286/1‑150 (MQ=255)
                                                                                                                            aTGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAg                            >  2:2686066/1‑150 (MQ=255)
                                                                                                                                ccACCGTGTCCCCCCCCCCCCTCCCCGGGTGTGGGTGGGGGTGAGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTg                        <  1:2342440/149‑1 (MQ=25)
                                                                                                                                  aCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCACACCCTAACCCACTACCCAAAAGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGa                      >  1:630304/1‑150 (MQ=255)
                                                                                                                                   ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTATCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATAAGTCTGGaa                     >  1:916509/1‑150 (MQ=255)
                                                                                                                                   ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCGAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                     >  1:3821656/1‑150 (MQ=255)
                                                                                                                                   ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCGGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                     >  1:3715449/1‑150 (MQ=255)
                                                                                                                                       tttCTCCCTCTCCACACCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCg                 <  1:4826877/150‑1 (MQ=255)
                                                                                                                                             cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           >  2:6126770/1‑150 (MQ=255)
                                                                                                                                             cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCGGGCCTGTTCCCCTCAACCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           >  2:1966116/1‑150 (MQ=255)
                                                                                                                                              ctctccctTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg          <  1:4454872/150‑1 (MQ=255)
                                                                                                                                              ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCTCACCCAAACCCACACCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTATTGATGTTTTCATAAGTATGGAAGCCGCACGTTg          >  1:4571878/1‑150 (MQ=255)
                                                                                                                                              ctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTAGGGAAGCCGCACGTTg          >  1:2543112/1‑150 (MQ=37)
                                                                                                                                               tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg         >  2:4296352/1‑150 (MQ=255)
                                                                                                                                                 tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                       <  1:333867/118‑1 (MQ=255)
                                                                                                                                                 tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                     <  2:6709712/136‑1 (MQ=255)
                                                                                                                                                 tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt       >  2:3377070/1‑150 (MQ=255)
                                                                                                                                                     ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt   >  2:3563102/1‑150 (MQ=255)
                                                                                                                                                     ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt   >  2:5006458/1‑150 (MQ=255)
                                                                                                                                                     ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt   >  2:1801688/1‑150 (MQ=255)
                                                                                                                                                     ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCACTCCCCCAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt   >  1:1459749/1‑150 (MQ=255)
                                                                                                                                                     ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGGCCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt   >  1:5688067/1‑150 (MQ=255)
                                                                                                                                                     atCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt   <  1:3584411/149‑1 (MQ=255)
                                                                                                                                                     atCCAGGGAGAGGGTCGGGGTGAGGGTAAGTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt   <  2:3581441/149‑1 (MQ=255)
                                                                                                                                                      tCCTGGGAGTGGGTCGGGTTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAtt  <  2:4066649/150‑1 (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: