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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,239,6110CA 9.6% 386.7 / 26.7 125intergenic (+66/‑333)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base C (75/38);  new base A (5/7);  total (80/45)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 126 total aligned reads displayed.

CGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATT  >  NC_000913/3239464‑3239758
                                                                                                                                                    |                                                                                                                                                   
cGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                                                                                                   >  2:1340672/1‑150 (MQ=255)
  tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                 <  2:2680945/150‑1 (MQ=255)
    gtgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                                                                                               >  1:831320/1‑150 (MQ=255)
    gtgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                                                                                               >  2:2623929/1‑150 (MQ=255)
     tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGgaga                                                                                                                                              <  2:2897059/150‑1 (MQ=255)
      gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGagag                                                                                                                                             <  1:1409039/150‑1 (MQ=255)
      gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGagag                                                                                                                                             <  1:2717830/150‑1 (MQ=255)
      gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGagag                                                                                                                                             <  1:2717832/150‑1 (MQ=255)
       tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  2:71179/1‑150 (MQ=255)
       tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  1:1620061/1‑150 (MQ=255)
       tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  2:571272/1‑150 (MQ=255)
       tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAgg                                                                                                                                            <  2:652059/150‑1 (MQ=255)
       tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACAGGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  2:1687471/1‑150 (MQ=255)
        ttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAggg                                                                                                                                           >  2:3033904/1‑150 (MQ=255)
         tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGt                                                                                                                                          >  2:1907502/1‑150 (MQ=255)
          ggCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                                         <  2:1442380/150‑1 (MQ=255)
          ggCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                                         <  1:53178/150‑1 (MQ=255)
          ggCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                                         <  2:1442342/150‑1 (MQ=255)
          ggCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                                         <  1:1131504/150‑1 (MQ=255)
           gCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCg                                                                                                                                        >  2:2105847/1‑150 (MQ=255)
           gCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCg                                                                                                                                        >  2:2403956/1‑150 (MQ=255)
           gCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCg                                                                                                                                        >  2:2803457/1‑150 (MQ=255)
           gCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCg                                                                                                                                        >  2:2835380/1‑150 (MQ=255)
                 tGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                  <  2:75813/150‑1 (MQ=255)
                 gggTTATGACTTTTATTATCAACGCGTGGGTTATTTATGGTCATTATAAGCTGGGGGGTTTATATTTTTTATTCTGCCTATGCCGTGTACCCTGTAATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCTGGGAGTGGGTCGGGGTGa                                                                                                                                  <  1:2565987/149‑1 (MQ=255)
                      atgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTa                                                                                                                             >  2:1489113/1‑150 (MQ=255)
                      atgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTa                                                                                                                             >  2:427039/1‑150 (MQ=255)
                      atgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTa                                                                                                                             >  2:928541/1‑150 (MQ=255)
                      atgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGGTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTa                                                                                                                             >  2:427095/1‑150 (MQ=255)
                            tttATTATAAACGCCTGGTTGAATTATCGGCATGAAAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCGTCCAGGGAGAGGGTCGGGGTGAGGGTAAttttt                                                                                                                       >  2:1612661/1‑150 (MQ=255)
                             ttattaTCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTACGCCGTGTACCCTGTCATTAACCTGAGCACCGTTTTCGCCCTCGCCCTCCCCGGGAGAGGGTGGGGGTGAGGGTAATTTTTc                                                                                                                      <  2:3114789/150‑1 (MQ=255)
                               attatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGc                                                                                                                    >  2:2926096/1‑150 (MQ=255)
                                ttatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCa                                                                                                                   >  2:2789149/1‑150 (MQ=255)
                                ttatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTCTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCa                                                                                                                   >  1:3038871/1‑150 (MQ=255)
                                     aaCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                                                                                                              >  1:2318970/1‑150 (MQ=255)
                                          cTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTgg                                                                                                         >  1:2436789/1‑150 (MQ=255)
                                           gggttGAATTAGCGCCATGATACGCCGGGGGGTGGATAATTTTAATTCTCCATACCCGGTGTACCCTGAATTTAAAATGCGCCCCGTTTTCGCCCTCCCCTTGCCCGGGAGTGGGTCGGGTTGTGGGTAATTTTTCGCACCGATGCTGGc                                                                                                        <  2:3078023/146‑1 (MQ=255)
                                              gTGAATTATCGGCATGATACGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTg                                                                                                     >  1:3181234/1‑150 (MQ=255)
                                                        ggCATGATAAGCAGCGGGGTGGATAATTAGTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCGTCCCGGGAGAGGGTCGGGGTGAGGGTAATTGTTCGCACCGATGCTGGCCGGTTCCCCTccc                                                                                           >  2:848422/1‑148 (MQ=255)
                                                         gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                          <  2:2121870/150‑1 (MQ=255)
                                                         gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                          >  2:720372/1‑150 (MQ=255)
                                                            tGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTa                                                                                       >  2:205214/1‑150 (MQ=255)
                                                              atacgatgcGGGGTGTTTATTTTTTATTCTGCACATGCCGTGTCCCCTGTCATTAACTTGAGCACGGTTTTCTCCCTCTCCGTTCCTGGTTGTGGTTGGGGGTGTGTGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAc                                                                                     <  2:2848366/143‑1 (MQ=255)
                                                                  gctggggggTGGTTTTTTTTTACTCTGCCTACGCCGTGTACCGTGTCTTTAACCTGCGCCCGGTTTTCGCCCTCTCCCTCCCGGGGTGGGGGGGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCctc                                                                                 <  2:2011975/145‑1 (MQ=19)
                                                                  gcagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCctc                                                                                 >  2:2548446/1‑150 (MQ=255)
                                                                   ccggggggTTTATATTTTTTAATCTCAATACGCCTTGACCCCTGTATTAAAATTGCCACCGGTTTTCCCCCTCCCCCTCCCCGGTTGTGGGTGGGGTTGTGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCctct                                                                                <  2:1663895/146‑1 (MQ=255)
                                                                        gggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCa                                                                           >  2:3139186/1‑150 (MQ=255)
                                                                        gggTGGATAATATTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCa                                                                           >  2:475146/1‑150 (MQ=255)
                                                                           tgtatatTTTTTAATGTGCATACGCGTTTTACCCTTTATTTAAATTTTGCACGTTTTTCTCCCTCTCCCTTCATTTGATTGTGTTGTTGTTTGTTTAATTTTTCGAACCGATGCTGGCCTGTTCTCCTCACCCTAACCCTCTCCCCAAAc                                                                        <  2:2838818/147‑1 (MQ=255)
                                                                             gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgg                                                                      >  2:1342758/1‑150 (MQ=255)
                                                                             gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgg                                                                      >  1:1143428/1‑150 (MQ=255)
                                                                              atatTTTTTAATCGCCCTAAGCCGTGTACCCTTAATTTAACATGAGCACGTTTTTCTCCCTCCCCCTCCCCGGGAGTGGGTCGGGGTTTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCGccccgtttggg                                                                     <  1:762157/150‑8 (MQ=21)
                                                                              aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                              >  2:2425622/1‑77 (MQ=255)
                                                                              aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                              >  2:2425637/1‑77 (MQ=255)
                                                                              aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                              <  1:2425622/77‑1 (MQ=255)
                                                                              aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                              <  1:2425637/77‑1 (MQ=255)
                                                                               tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCATGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgggg                                                                    >  1:1913806/1‑150 (MQ=255)
                                                                               tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACGCTAACCCTCTCCCCAAACgggg                                                                    >  1:412950/1‑150 (MQ=255)
                                                                               tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgggg                                                                    >  1:2147617/1‑150 (MQ=255)
                                                                               tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgggg                                                                    >  2:734529/1‑150 (MQ=255)
                                                                               tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCGGTTCCCCTCACCCTAACCCTCTCCCCAAACgggg                                                                    >  2:3080821/1‑150 (MQ=255)
                                                                                attttttACTCGGCCTCCGCGGTGCCCCCTGTCTTTACCTGGCGCCCGGTTGCCGCCCTCGCCCTCCCGGGGTGGGGGGGGGGGTGTGGGCAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGc                                                                   <  2:359704/148‑1 (MQ=12)
                                                                                   ttttAATCTGCCTAAGCCGTGTACCCTGTCATTACCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAg                                                                <  2:3231757/150‑1 (MQ=255)
                                                                                   ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAg                                                                >  2:1948379/1‑150 (MQ=255)
                                                                                    tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGGGAGGGTAATTTTTCGCACCGATGCGGGCCTGTTCCCCGCACCCTAACCCTCTCCCCAAACGGGGCGCgg                                                               >  1:1102808/1‑150 (MQ=255)
                                                                                    tttAATCGGCATAAGCCGTGAACCCTGACATTAACATGAGCACCGTTTTCGCCCCCGCCCTTCCTGGTTGTGGGTGGGGTTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAgg                                                               <  1:701731/150‑1 (MQ=25)
                                                                                        aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGAc                                                           >  1:2914143/1‑150 (MQ=255)
                                                                                             ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCACACGGGGCGAGGGGACTGAcc                                                      >  2:1583538/1‑150 (MQ=255)
                                                                                              cTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCg                                                     >  2:1248112/1‑150 (MQ=255)
                                                                                               tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                    >  1:38195/1‑150 (MQ=255)
                                                                                               tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                    >  1:811428/1‑150 (MQ=255)
                                                                                               tAAGCCGTGTACCCTGTCATTAACATGAGCAACGTTTTCTCCCTCTCCCTTACAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAAACCTCACACCAAACGGGGCGAGGGGACGGACCGa                                                    >  1:38190/1‑150 (MQ=255)
                                                                                                aaGCGGTGTACCCTGTCATTAACATGAGCACCGTTTTCGCCCTCTCCCTTCCTGGGTGTGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAg                                                   <  1:752873/150‑1 (MQ=37)
                                                                                                         aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  2:1597046/150‑1 (MQ=255)
                                                                                                            gtgacATTAACATGAGCACGTTTTTCTCCCTCTCCCTTCCCGGTTGAGGGTGGGGGTGTGTGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                       <  2:2676571/146‑1 (MQ=25)
                                                                                                            cTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                       <  2:795781/150‑1 (MQ=255)
                                                                                                            cTGTCAGTAAGCTGAGCACCGTTTTCTCCCTCGCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                       <  2:1911769/150‑1 (MQ=255)
                                                                                                             tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtg                                      >  2:2263860/1‑150 (MQ=255)
                                                                                                             tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtg                                      >  2:1737693/1‑150 (MQ=255)
                                                                                                             tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCGGTTCCCCTCACCCGAACCCACTCCCCAAACCGGGCGAGGGGACTGCCCGCGGCCTTTTTTGAtg                                      >  1:298160/1‑150 (MQ=255)
                                                                                                                 aTTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCTCACCCTAACCCTCTCCCCAAACCGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgttgt                                  >  1:1233462/1‑150 (MQ=255)
                                                                                                                          aGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCt                         >  2:2651392/1‑150 (MQ=255)
                                                                                                                          aGCACCGTTTTCGCCCTCTCCCTTCCTGGTAGTGGTTGGGGGTGTGCGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCt                         <  1:2123084/150‑1 (MQ=37)
                                                                                                                            cACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCACTCCCCAAACCGGGCGAGGGGACGGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTgg                       >  1:173351/1‑150 (MQ=255)
                                                                                                                              ccGTTTTCTCCCTCTCCCTTCCAGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                     <  1:2936620/150‑1 (MQ=255)
                                                                                                                              ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                     <  2:1356598/150‑1 (MQ=255)
                                                                                                                                    ttCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGc                >  1:2946912/1‑150 (MQ=255)
                                                                                                                                     tctccctctcccTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGCCTGACCGAGTCCTTTTTTGATGTTGTCATCAGGCTGGAAGCCGCa               >  2:1727432/1‑150 (MQ=255)
                                                                                                                                      ctcccccccccgGCCGGGGTGTGGGTCGGGGTTAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCAc              <  2:2565233/150‑1 (MQ=37)
                                                                                                                                      cgccccctccctGCCGGGGTGCGGGTGGGGGGGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCAc              <  2:2561408/148‑1 (MQ=37)
                                                                                                                                       tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACg             >  1:2009998/1‑150 (MQ=255)
                                                                                                                                         cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           >  2:2044087/1‑150 (MQ=255)
                                                                                                                                         cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAACCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           >  1:693332/1‑150 (MQ=255)
                                                                                                                                         cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCTCACCCTAAACCTCTACCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTTTTTGTCATCAGTCTGGAAGCCGCACGtt           >  1:693308/1‑150 (MQ=255)
                                                                                                                                         cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCTCACCCAAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCATTTTTTATTTTGTCATCAGTCTGGAAGCCGCACGtt           >  1:1522955/1‑150 (MQ=255)
                                                                                                                                             tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTGCTTTTTTGTTGTTGTCATCAGTCTGGAATCCGCACGTTgggt       >  2:1461118/1‑148 (MQ=255)
                                                                                                                                             cccctCCCGGGGTGCGGGTGGGGGGGAGGGTCATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt       <  1:2283407/149‑1 (MQ=37)
                                                                                                                                              ccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCtt      >  2:2650171/1‑150 (MQ=255)
                                                                                                                                                  tCCAGGGAGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAtt  <  1:2092582/150‑1 (MQ=255)
                                                                                                                                                  gccAGGGCGAGGGTCGGGGTGAGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAtt  <  1:2092616/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATT  >  NC_000913/3239464‑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: