Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,822,1720CT27.6% 112.4 / 18.5 58D29N (GAC→AAC) osmEosmotically‑inducible lipoprotein OsmE
Reads supporting (aligned to +/- strand):  ref base C (30/12);  new base T (0/16);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 3.80e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAGCCTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCGCTTGTCGTTTCAGTAGCAAC  >  NC_000913/1821924‑1822411
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
cGCAGCCTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCtt                                                                                                                                                                                                                                                  <  1:73696/251‑1 (MQ=255)
  cAGCCTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCtttt                                                                                                                                                                                                                                                >  2:225467/1‑250 (MQ=255)
        gtgGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTAc                                                                                                                                                                                                                                          <  2:578267/249‑1 (MQ=255)
            aTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCAt                                                                                                                                                                                                                                    >  2:641529/1‑251 (MQ=255)
              cAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATa                                                                                                                                                                                                                                   <  2:318457/250‑1 (MQ=255)
                     gTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTAc                                                                                                                                                                                                                           >  1:313087/1‑251 (MQ=255)
                      tATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTAc                                                                                                                                                                                                                           >  1:538165/1‑250 (MQ=255)
                      tATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTAc                                                                                                                                                                                                                           >  2:559387/1‑250 (MQ=255)
                       aTTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACaa                                                                                                                                                                                                                         <  1:219601/251‑1 (MQ=255)
                          cAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATCAGCCGTACAAccc                                                                                                                                                                                                                      >  2:266367/1‑251 (MQ=255)
                                gggTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTCCCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCACCCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACGGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGc                                                                                                                                                                                                                <  1:160922/250‑1 (MQ=255)
                                  gTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTCCCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTCCAACCCGCCAGCAt                                                                                                                                                                                                              <  2:585045/251‑1 (MQ=255)
                                    cTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATgg                                                                                                                                                                                                            >  1:426553/1‑251 (MQ=255)
                                    cTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATgg                                                                                                                                                                                                            >  1:143285/1‑251 (MQ=255)
                                     tGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGt                                                                                                                                                                                                           >  2:208786/1‑251 (MQ=255)
                                         aaCCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTGGACCCAGGAGGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCGCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAAt                                                                                                                                                                                                       <  2:480994/251‑1 (MQ=255)
                                                    atgaCATGTACGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGACTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGc                                                                                                                                                                                           >  2:617884/1‑251 (MQ=255)
                                                          tGTCGGGTATCATTTAACGCGCAAAAGTGGGTTTCTGCTTTACCACCACGTGGACCCGGGATGTAGGTCGGGCAAGAACCGCAGGCTGGGATCATGCGCACTTCAGACGAAGTTTTACCCGCAATTTGGGAACCCGGCCCCCGGCTAATGCCTTTTTGGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTCCGATCATAAGCCGTACACCCCGCCAGCATGTTTAATACCGCCGCTGCACTc                                                                                                                                                                                       <  2:576008/249‑1 (MQ=255)
                                                            tCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGa                                                                                                                                                                                    >  2:432948/1‑250 (MQ=255)
                                                                gTATCATCTAACGCGACAAAGTAGGTTGCTGCTTTCCCAGCACGTTGCCCCAGGATGTAGGTCGGGCAAGTACCGCGACCATGGATCATGCTCACTTCAGAGGAGGGTTTACCCGAAATTGGGCCAACCGGCGCCCGGCTCGTGCCTTTTTTGACGTTTTGCACCACGGGCTGGACAACCTGGTTTTTGGTATGATCATAAGCCGTACAACCCGCCAGCATTGTTAATACCGCCGCTGCACTCAGAATTcc                                                                                                                                                                               <  1:161062/251‑1 (MQ=255)
                                                                 tATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCt                                                                                                                                                                              >  2:345735/1‑251 (MQ=255)
                                                                    catcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGcc                                                                                                                                                                           >  1:573253/1‑251 (MQ=255)
                                                                     atcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTAGCAAGTCCCGCGAGCATGAATCGTGCTCACTGCAGACGAGGGTTTGCCCGCAATTTGTGCAGCCTGCGCCCGGCTCATGCCTTTTTTGCCGTCTTACCCCACGGTCTGTACAAACGGTTCTTTGGTACAATCATAAGCCGTACAACCCGCCAGCAGGGTTAATGCCGCCGCTGCACTCAGAATTCCTGcc                                                                                                                                                                           <  1:4418/250‑1 (MQ=255)
                                                                        aaaCGCGACAAAGTAGTTTTCTTCTTTACCAACACGTTCACCCAGGATGTGGGTCGGGCAAGTCCCGCGAGCATGGTTCATGCTCACTTCAGACGAAGGTTTCCCCGCAATTTGTGCACCCTCCGCCCGGCTCATCCTTTTTTTGACGTCTTTCCCCACAGGCTGTACAACCTGTTTTTTGTTACAATCAAAAGCCGTACAACCCGCCAGCAGGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCatat                                                                                                                                                                       <  2:621048/250‑1 (MQ=255)
                                                                                aaaaGTAGGTTTCGGCTTTACCACCACGTTGCCCCAGGATTAAGGGCTGCAAAGTACCGCGAGCATGGATCATCCTCCCTTCAGACGAAGGTTTCCCCCCAATTTTGCAACCCGGCGCCCGGCTCATGCTTTTTTTGTCGTCTTTCACCACGGGCTGTACAAACGGGTCTTTGGTACGACCATAAGCCTTACAACCCGCCACCAGGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTc                                                                                                                                                               <  1:423544/250‑1 (MQ=255)
                                                                                        gTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCg                                                                                                                                                       >  2:168170/1‑251 (MQ=255)
                                                                                         tttCTGTTTTCCCATCCCGTTGCCCCAGGATGTAGGTCTGGCAAGTACCGCGAGAATGGATCATGCTCACTTCAGACGAGGGTTTCCCCGCAATTTGTGCACCCTGCGCCCGGCTCATGCTTTTTTGGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCg                                                                                                                                                       <  1:345735/250‑1 (MQ=255)
                                                                                              gCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCctct                                                                                                                                                 >  1:490329/1‑251 (MQ=255)
                                                                                              gCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCctct                                                                                                                                                 >  2:139012/1‑251 (MQ=255)
                                                                                                 ttACCATCACGGTGACCCAGGATGTACGTCTGGCAAGTACCGCGAACATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGCCGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTCCAACCCGCCAGCATGGTTACTACCGCCGCTGCACTCAGAATTCCTACCATATTCTTGTTCATATTCCGTCATCTTGt                                                                                                                                              >  2:580284/1‑251 (MQ=255)
                                                                                                           gttgccccaggtggatggTCTGGCAAGTACCCCGAGAATTCTTCATGCTCCCTTCATACGAAGGTTTCCCCGCATTTGTTACAACCTGCCCCCGGCCATGGTCTTTTGTTACGTTTTTCACCACAGGTTGTACAACTTGTTTTTTGGTCCGATAATAAGCCGTACACCCCGCCAGCATTGTTAATCCCGCCGCTGCACTCATAATTCCTGCCATATTCTTGTTCATATTCCGGCCTCGTGTTTATCAGCgt                                                                                                                                    <  2:338016/235‑1 (MQ=255)
                                                                                                           gTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCatat                                                                                                                                                           >  1:363895/1‑228 (MQ=255)
                                                                                                                        gaagttcgggCAGGTACCGCGAGCATGGTTCATGCTCACTTAATCCGAAGTTTTACCCGCATTTTGTTCACCCTGCGCCGGGCTTTTGCTTTTTTTAACGTCTTTCACCACAGGCTGTACAAACTGGTTTTTGGTACGATAATAAGCCGTACACCCCGCCAGCATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCt                                                                                                                       <  1:270647/243‑1 (MQ=255)
                                                                                                                          aGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTgg                                                                                                                     >  2:80409/1‑251 (MQ=255)
                                                                                                                                         cgcgGGCATGGACCATGCTCCCTTCAGACGAAGGTTTACCCGCAATTTGTGCACCCTGCCCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTTTTTGGTACGATAATAAGCCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCt                                                                                                      <  1:309696/251‑1 (MQ=255)
                                                                                                                                            gagaagggatcatgccccctcaaGACGAAGGTTTACCCGCAATTGGTGCACCCTGCGCCCGGCTAATGCCTTTTTTGACGTTTTTCACCACGGGCTGTAAAACCGGGTTTTTGTTACAATCATAAGCCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCtttt                                                                                                   <  1:169913/229‑1 (MQ=255)
                                                                                                                                            gAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCtttt                                                                                                   >  1:550191/1‑251 (MQ=255)
                                                                                                                                                  ggATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGc                                                                                             >  2:638139/1‑251 (MQ=255)
                                                                                                                                                          cccACTTCAGACGAAGGTTTCCCCCAATTTGGTGAACCCGGCGCCCGGCTAATGCCTTTTTTGCCTTCTTTCCCCACGGGCGGTAAAAACGGGTTTTTGGTACGATCATAGGCCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAAc                                                                                     <  1:173541/249‑1 (MQ=255)
                                                                                                                                                           tCACTTCAGACGAAGGTTTCCCCGCAATTTGTGAAACCTGCGCCCGGTTCATGCCTTTTTTGCCTCTTTCCACCACAGGCTGTACAACCTGGTTTTTGTTACGATCAAAAGCCGTACACCCCGCCAGCATGGTTAATCCCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGaa                                                                                      <  1:391839/249‑1 (MQ=255)
                                                                                                                                                                tCAGAGGAGGGTTTCCCCGCAATTTGTGCACCCTGCGCCCGCCCCTTGCCTTTTTGCACGTTTTCCACCACGGCCTGTAAAACCGGTTTTTGGGTCCGATAAAACCCCGAACACCCCCCAAGCATGTTTAATCCCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAAc                                                                                <  2:313087/250‑1 (MQ=255)
                                                                                                                                                                      gAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCg                                                                         >  1:238260/1‑251 (MQ=255)
                                                                                                                                                                        aGGTTTACCCGCAATTTGTGCACCCGGCCCCCGCCTAATGCCTTTTTTGACGTTTTTCCCCACAGCCTGTACAACCTGGTTTTTGGTACGATCAAAAGCCGTCCACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCt                                                                       <  1:559387/251‑1 (MQ=255)
                                                                                                                                                                              ccccGCAATTTTGGCACCCTGCGCCCGGCTAATGCTTTTTTTGCCTTTTTTCCCCACAGGTTGTACAACCTGGTTTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCag                                                                  <  1:208786/249‑1 (MQ=255)
                                                                                                                                                                                 cGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCCTATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAggc                                                              >  2:38172/1‑250 (MQ=255)
                                                                                                                                                                                    aaTTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTc                                                           >  1:535495/1‑251 (MQ=255)
                                                                                                                                                                                       ttGTGCACCCGGCGCCCGGCTCATGCTTTTTTTGACGTCTTCCACCACAGGCTGTACAAACGGTTTTTGGTTACGATAAAAAGCCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCt                                                         <  1:432948/250‑1 (MQ=255)
                                                                                                                                                                                                     gcgcCCGGCTCAGGCCTTTTTTGACGTCTTCCCCCCCAGGCTGTACAAACGGGTTTTGGGTACGATCATAAGCCTTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGtttt                                             <  1:168170/251‑1 (MQ=255)
                                                                                                                                                                                                     gcgcCCGGCTAATGCTTTTTTTGACGTCTTCCACCCCAGGCTGTACAAACTGGTTTTTGTTACGATAATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGtt                                               <  1:122484/249‑1 (MQ=255)
                                                                                                                                                                                                        cccGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAac                                          >  1:108766/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGTCGCTTTTTCAAGCCCGTGAACGAAACGGCACCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCg                      >  2:616374/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ttcccccacAGGCTTAAAAAACGGGTTTTTGGTACAATCAAAACCCTTAAACCCCCCCAGCATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCGCTTg                  <  1:38172/245‑1 (MQ=255)
                                                                                                                                                                                                                                   ttCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCGCTTGTc                >  2:640483/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cctcttccccaatctaaatcTGGTTTTTGTTACTATCAAAAGCCGTACACCCCGCCAGATTGTTAATTACCGCCGCTGCACTCAGAATCCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCGCTTGTCg               <  2:339122/232‑1 (MQ=255)
                                                                                                                                                                                                                                           caGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTg                         >  1:71606/1‑235 (MQ=255)
                                                                                                                                                                                                                                                tGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATACGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCGCTTGTCGTTTCAGTAGc     >  1:413209/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   aaaaaCGGGTCTTTGGTACGATCATAAGCCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCGCTTGTCGTTTCAGTAGCAAc  <  1:139012/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    cAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCGCTTGTCg               >  1:339122/1‑236 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
CGCAGCCTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCGCTTGTCGTTTCAGTAGCAAC  >  NC_000913/1821924‑1822411

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: