Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,239,609 T→G 5.6% 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 5.6% 835.7 / 12.7 286intergenic (+64/‑335)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base T (185/60);  new base G (4/16);  total (189/97)
Fisher's exact test for biased strand distribution p-value = 1.05e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

TGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATT  >  NC_000913/3239460‑3239758
                                                                                                                                                     |                                                                                                                                                     
tGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCtt                                                                                                                                                       <  1:4569514/150‑1 (MQ=255)
  ggCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTcc                                                                                                                                                     >  1:4344018/1‑150 (MQ=255)
   gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  1:3413502/150‑1 (MQ=255)
    cGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                                                                                                   >  1:2958240/1‑150 (MQ=255)
    cGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                                                                                                   <  2:6600865/150‑1 (MQ=255)
      tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                 >  2:2075943/1‑150 (MQ=255)
      tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                 <  1:1422986/150‑1 (MQ=255)
      tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACAAGAGCACCGTTTTCTCCCTCTCCCATCCAggg                                                                                                                                                 >  1:6202048/1‑150 (MQ=255)
      tcgtctttGGCATTCTGGAGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGAACCCTGTCATTAACATGAGCACCGATTACTCCCTCTCCCTTCCAgcg                                                                                                                                                 >  2:3719156/6‑148 (MQ=255)
       ggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                                                                                                >  1:4817070/1‑150 (MQ=255)
        gtgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCGTCCAGGgag                                                                                                                                               >  2:2444372/1‑150 (MQ=255)
          gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                                                                                                             <  2:82163/150‑1 (MQ=255)
           tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            <  2:9562/150‑1 (MQ=255)
           tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  2:557388/1‑150 (MQ=255)
              ggCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                                         <  1:999618/150‑1 (MQ=255)
               gCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAgggaga                                                                                                                                        <  1:3260499/150‑4 (MQ=255)
                  ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                                                                                                     >  1:6346414/1‑150 (MQ=255)
                     tGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                  <  2:5951914/150‑1 (MQ=255)
                      ggTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAg                                                                                                                                 <  2:3913618/150‑1 (MQ=255)
                       gTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAgg                                                                                                                                <  2:408532/150‑1 (MQ=255)
                        tgatgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCAGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAggg                                                                                                                               <  2:4653581/150‑1 (MQ=255)
                        agatgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGACTAAGCCGTGTACCCTGTCAATAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAggg                                                                                                                               >  2:3750057/2‑150 (MQ=255)
                            gaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAt                                                                                                                           <  2:315481/150‑1 (MQ=255)
                             aCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtt                                                                                                                          >  1:3855179/1‑150 (MQ=255)
                               gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                        >  1:4033302/1‑150 (MQ=255)
                               gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                        >  1:4303607/1‑150 (MQ=255)
                                  tattatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCg                                                                                                                     >  1:6294190/1‑150 (MQ=255)
                                  tattatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCg                                                                                                                     >  1:2155860/1‑150 (MQ=255)
                                  tattatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCg                                                                                                                     >  1:5150567/1‑150 (MQ=255)
                                      atCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAcc                                                                                                                 >  1:696320/1‑150 (MQ=255)
                                         aaCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                                                                                                              >  2:6014347/1‑150 (MQ=255)
                                          aCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATg                                                                                                             >  2:6085984/1‑150 (MQ=255)
                                              cTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTgg                                                                                                         >  1:1915941/1‑150 (MQ=255)
                                               tGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGc                                                                                                        >  2:1105363/1‑150 (MQ=255)
                                                gggTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGcc                                                                                                       >  2:3011054/1‑150 (MQ=255)
                                                  gTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTACTTTTTCGCACCGATGCTGGCGTg                                                                                                     >  2:4597411/1‑150 (MQ=255)
                                                  gTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTg                                                                                                     >  1:849499/1‑150 (MQ=255)
                                                    gAAGTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCACCCTTCCAGGGAGAGGGTCGCGGTGAGGGTAATATTTCGCACCGATGCAGGCCTGtt                                                                                                   >  1:4150428/1‑150 (MQ=255)
                                                     aaTTATCGGCATGATACGCCGCGGGGTGGATATTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACTTGCGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGTGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTc                                                                                                  <  2:6445906/150‑1 (MQ=255)
                                                       ttATCGGAATTATAAGCTGGGGGGTGTTTAATTTTTAATGTGCCTAAGCCGTGTACCCGGCCTTTAACATGAGCACCGTTTTCGCCCTCTCCCTTCCTGGGTGTGGGTGGGGGTGAGGTTAATTTTTCGCACCGATGCTGGCCTGTTccc                                                                                                <  2:4650849/150‑1 (MQ=255)
                                                       ttAACGGCATGATAAGCAGCGGGGTGGATAATATTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTccc                                                                                                >  1:3698973/1‑150 (MQ=255)
                                                          tCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTccccgc                                                                                             >  2:2808672/1‑148 (MQ=255)
                                                            ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAc                                                                                           <  2:1644008/150‑1 (MQ=255)
                                                            ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAc                                                                                           >  1:5253243/1‑150 (MQ=255)
                                                             gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                          >  2:3456715/1‑150 (MQ=255)
                                                             gCATGATAAGCAGCGGGGTGGATAATTTTTAATCAGCCTAAGCCGTATACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                          >  2:6167889/1‑150 (MQ=255)
                                                              cATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAccc                                                                                         >  1:2252907/1‑150 (MQ=255)
                                                                    aaGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAccc                                                                                   >  2:3235790/1‑150 (MQ=255)
                                                                     agcagcGGTGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACTGTTTTCTCCCTCTCCCTTCCTGGGAGTGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCct                                                                                  <  2:4323303/150‑1 (MQ=38)
                                                                      gcagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCGTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCGGTTCCCCACACCCTAACCctc                                                                                 >  2:2125582/1‑150 (MQ=255)
                                                                       cagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCctct                                                                                >  1:4066115/1‑150 (MQ=255)
                                                                        tggggggTGTATATTTTTTATTCTGCGTACGCGGTTTACCTTTTCTTTAAATTTCGCCCGTTTTTCTCCCTCTCCCTGCCTGGGTGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                               <  2:2621439/147‑1 (MQ=12)
                                                                        gggggggTGTATTTTTTTTATTGGGCCTATGCGGTGTACCGTGTCTTTAACTTGCGCCCGGTTTGCGCCCTCGCCGTGCGGGGGTGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                               <  2:5734965/147‑1 (MQ=16)
                                                                         gcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTcc                                                                              >  1:5800708/1‑150 (MQ=255)
                                                                               tgtttattttttATTGTGCCTAAGCCGTGTACCGTGTCTTTAAACTGCGCCCCTTTTGCTCCGACTCCCTCCCGGGGTGTGGGTGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAc                                                                        <  2:3610685/143‑1 (MQ=16)
                                                                               tgtaTAATTTTTATTCTGCCTACGCGGTGTACCCTGTCATAAACTTGAGCACCGTTTGCCCCCTCCCCCTGCCCGGGCGCGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAc                                                                        <  2:3421613/147‑1 (MQ=25)
                                                                                  atatTTTTTAATCTCCCAACGCGTTGCACCCTTACTTTAACATGAGCCCCGTTTTCTCCCTCCCCCTCCCTGGTTGTGGTTGGGGGTGTGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACggg                                                                     <  2:6222337/150‑1 (MQ=19)
                                                                                   tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgggg                                                                    >  1:1237445/1‑150 (MQ=255)
                                                                                   tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTGTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCGGGCCTGTTCCCCTCCCCCTAACCCTCTCCCCCCACgggg                                                                    >  2:2723414/1‑150 (MQ=255)
                                                                                      tttttAATCTGCATAAGCGGTGTACCCTGTATTTAAATTGCCCACCTTTTTCACCCTCTCCCTTCCTGGTTGTGGGTGGGGTTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGa                                                                 <  2:2030172/150‑1 (MQ=21)
                                                                                           acggtGCCTACGCGGTGCACCCTTTATTTAACTGGCGCACCGTTTACCCCCCCCCCCTTCCGGGTTGTGGGGGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGa                                                            <  1:128652/146‑1 (MQ=12)
                                                                                              cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTg                                                         >  2:4393393/1‑150 (MQ=255)
                                                                                              cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTg                                                         >  1:1168947/1‑150 (MQ=255)
                                                                                               tGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGa                                                        >  1:2588617/1‑150 (MQ=255)
                                                                                               tGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGa                                                        >  1:6536171/1‑150 (MQ=255)
                                                                                               tGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCACCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCATCTCCCCAAACGGGGCGAGGGGACTGa                                                        >  1:1006035/1‑150 (MQ=255)
                                                                                                  cTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCg                                                     >  1:1097877/1‑150 (MQ=255)
                                                                                                   tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                    >  1:3619216/1‑150 (MQ=255)
                                                                                                   tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGa                                                    >  1:272669/1‑150 (MQ=255)
                                                                                                   tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCGGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                    >  2:6349468/1‑150 (MQ=255)
                                                                                                    ccgCCGTGTCCCCGGGCTTTACCCGGAGCCCCGTTTTCGCCCGCTCCGTCCCGGGGTGCGGGTGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAg                                                   <  2:5228516/148‑1 (MQ=18)
                                                                                                     aGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGt                                                  >  1:2096525/1‑150 (MQ=255)
                                                                                                         gtgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtt                                              >  1:4126361/1‑150 (MQ=255)
                                                                                                           ttACCCTGTAATTAACATGCGCCCCGTTTCCCCCCCCTCCCTTCCCGGGTGTGGGTGGGGGTGAGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttt                                            <  1:5369835/149‑1 (MQ=21)
                                                                                                             aCCCTGTCATTAACTTGAGCACGTTTTTCTCACTCGCCCTACCAGGGTGTGGGTGGGGGGGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  2:1500939/150‑1 (MQ=21)
                                                                                                             aCCCTGTCATTAACATGTGCACGTTTTACTCCCTCGCCCTCCCCGGGTGTGGGTGGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  2:431150/150‑1 (MQ=21)
                                                                                                               ccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCGTCCAGGGAGTGGGTCGGGGTGAGGGTAATGTTTCGCACCGATGCTGGCCTGTTCCCCTCTCCCTAACCCTCTCCCCCCACGGGGCGAGGGGACGGACCGAGTCCTGTtgtgt                                        >  2:4805447/1‑149 (MQ=255)
                                                                                                                cTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTCATTTTTCGCGCCGATGCTGGCCTGTTCCCCTCACCCTAGCCCGCTCCCCAAACGGGGCGGGGGGACTGACCGAGGCCTTTTGTGAt                                       >  2:4772522/1‑150 (MQ=255)
                                                                                                                  gTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgt                                     >  1:5439607/1‑150 (MQ=255)
                                                                                                                   tCATTAACTTGAGCACCGTTTTCTCCCTCTCCCTGCCCGGGAGAGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgtt                                    <  2:5370277/150‑1 (MQ=255)
                                                                                                                   tCATTAACATGAGCACCTTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgtt                                    <  2:4239865/150‑1 (MQ=255)
                                                                                                                   tCATTAACATGAGCACCGTTTTCTCCCTCTCCCGTCCCGGTTGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgtt                                    <  1:4954799/150‑1 (MQ=255)
                                                                                                                      ttAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGGTGGCCTGTTCCCCTGACCCTAACCCTCTCCCCAAAAGGGGCGAGGGGGCTGACCGAGTCCTTTTTTGAtgttttg                                 >  1:283844/1‑147 (MQ=255)
                                                                                                                        aaCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcat                               >  2:228500/1‑150 (MQ=255)
                                                                                                                          cATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcatca                             >  2:5028658/1‑150 (MQ=255)
                                                                                                                          cATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCTCTCCCCAAACGGGGTGAGGGGAATGACCGAGTCCTTTTTTGATTTTGtcatca                             >  1:1918398/1‑150 (MQ=255)
                                                                                                                          cATGACCACCGTTTACACCCACCCCCTACCCCGGACAGGTTCGGGTTTAGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACCGGGCGAGGGGACTGACCGAGTCCTTTTTTGAAGTTGtcatca                             <  2:892066/150‑1 (MQ=25)
                                                                                                                           aggcgCACGGTTTTCTCCCTCTCCCTCCCTGGGTGTGGGTGGGGGTTTGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAg                            <  1:3737334/146‑1 (MQ=25)
                                                                                                                              cccACCGTTTTCACCCACCCCCTCCCGGGTCGTGGGTCGGGTTGTGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCt                         <  1:2526062/148‑1 (MQ=25)
                                                                                                                               ccccggTTTGCCCCCCCCCCCTCCCCGGGTGTGGGGGGGGGGGCGGGCAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTg                        <  1:17969/145‑1 (MQ=21)
                                                                                                                                cACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTgg                       >  2:1355036/1‑150 (MQ=255)
                                                                                                                                cACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCTCACCCAAACCCACACCCCAAACCGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTgg                       >  1:5052362/1‑150 (MQ=255)
                                                                                                                                  cctttttCTCCCTCTCCCTTCCCGGGTGCGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                     <  1:3919315/147‑1 (MQ=255)
                                                                                                                                    gTTTCCGCCCACCCCCTGCCCGGGAGTGGGTGGGGGTGCGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGc                   <  1:6239818/150‑1 (MQ=25)
                                                                                                                                     ttttCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTTGTTGTCATCAGTCTGGAAGcc                  >  2:733486/1‑150 (MQ=255)
                                                                                                                                          tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACg             >  2:1965369/1‑150 (MQ=255)
                                                                                                                                            cctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           >  1:5988632/1‑150 (MQ=255)
                                                                                                                                            cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           >  1:394855/1‑150 (MQ=255)
                                                                                                                                             ctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCGGGCCTGTTCCCCTCCCCCTAACCCTCTCCCCACACGGGGCGAGGGGACTGACCGAGTCCGTTTTTGCTGTTGGCAGCAGTCTGGCAGCCGCACgggg          >  2:1851295/1‑147 (MQ=255)
                                                                                                                                                    ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt   >  2:1211135/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: