Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,967,2820TC26.0% 79.2 / 10.6 51G52G (GGA→GGG) cheYchemotaxis protein CheY
Reads supporting (aligned to +/- strand):  ref base T (24/13);  new base C (0/13);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 3.63e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTT  >  NC_000913/1967036‑1967503
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                              
atAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCAt                                                                                                                                                                                                                            >  1:345051/1‑251 (MQ=255)
          ccTCACATGCCCAGTTTCTCCAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑gcctgc                                                                                                                                                                                                                 >  1:91873/1‑251 (MQ=255)
            tCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAAAACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑gcctgcct                                                                                                                                                                                                               >  1:307919/1‑251 (MQ=255)
                   cccAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAACGTCGCGGCGGTAAATGGCTTCAACACATAGACACTGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CCCCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACtt                                                                                                                                                                                                        >  1:306145/1‑251 (MQ=255)
                    ccAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACtt                                                                                                                                                                                                        >  1:411191/1‑250 (MQ=255)
                     cAGTTTCTAAAAGATTTGTTTGAGTTTTTCCCCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACGGGCCCCCGC‑TTGCGCCGCAGCAATGATGTCCTCTTTCTTCGCTTCTGCAGT‑CACCATAACCACTGGCAATCCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGCCATGTTCCAGTCGGAGATAACAAATCCAGAACC‑GCCGGCCTGCAACTTAt                                                                                                                                                                                                      <  1:361394/251‑1 (MQ=255)
                        tttCTCAAAGATTCTCTTGCGTTTTTACTCCAGCGTCGCGGCGGTCAATGCCTTAACCAAATAGCCACTGGACCCCGCCTTG‑ACCGCCGCAATGTCGTTCTCTTTCTTCGCGTCTGCAGT‑CAACATTCACACTGGCAATGCCCACATCTCGCCCTCCGCACGAATTGTTTTCAGCAATTCCAGGCCCTCCATAATCGGCATGTTCCCGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTCTTg                                                                                                                                                                                                    >  2:404905/1‑250 (MQ=255)
                              aaaGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGC‑TTGGGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAAGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCATCAACTTATTGAGAGCgg                                                                                                                                                                                             >  1:160576/1‑250 (MQ=255)
                                  aTTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGAc                                                                                                                                                                                         >  1:374234/1‑251 (MQ=255)
                                    ttttttgAGTTTTCCCCCCCGCGCCGCGGCGGTAAAGGGCTACCCCACATAGCCCCGGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTCCGCTTCTGCAGT‑CCCCATTAACCCTGGCAATGCCGACATCGCGCCACCCCCACAATTTGTTTTCAGCAATCCCAGGCCATCCATATGGGGCATGTTCCAGTCGGAGATAACAAATCCAAAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACg                                                                                                                                                                                        <  2:377686/246‑1 (MQ=255)
                                              ttttCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAACCACTGGCCCCCGC‑TTGCGCCGCATCAATGATGTTCTCTTTCTTCGCTTCTGCAGA‑CACCATTCACACTGGCACTGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCCATTCCAGGCCATCCATATTGGGCATGTTCCAGTAGGAGATAACACATCCATAACC‑GCCTGCCTGCCACTTATTGCGAGCGTCGACGcccacctccgc                                                                                                                                                                             >  1:32235/1‑242 (MQ=255)
                                                  cctccAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTcc                                                                                                                                                                         >  2:331494/1‑251 (MQ=255)
                                                       aGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCCTCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCCCTTCCTCacc                                                                                                                                                                    >  2:4922/1‑249 (MQ=255)
                                                             gcggcgGTAAAGGGCTTCCCCCAATAGCCCCTGCCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACGGGAAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACattatt                                                                                                                                                              <  1:29163/251‑1 (MQ=255)
                                                                gcgatAAATGCCTTCACCACATAGCCTCTTGCCACAGA‑TTTCCCCGCAGCAAGTATGTTCTTTTTCTTCGCTTCCGCAGT‑CACCATTAACACGGCCAATCCCCACACCGCGCAATCCGCACGAATTGTTTTCAGCAATACCAGGCAATCCACATTGGGCATGTTCAAGTCGGAGTTACCAAATCCAAAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCGCAACATAATTGaa                                                                                                                                                           <  1:88829/247‑1 (MQ=255)
                                                                     aaagggCTTCACCACATAGCCCCTGGCCCCCGC‑TTGCCCCGCAGCAATGATGTCCTCTTTCTCCGCTCCTGCAGT‑CCCCATTAACCCTGGCAATCCCGACATCGCCCCATCCGCACGAATTGTTTTCAGCAATCCCAGCCCATCCATATTGGGCATGTTCCAGTCGGAGAAAACAAACCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCa                                                                                                                                                      <  1:265513/247‑1 (MQ=255)
                                                                       aTGGCTTCACCACATAGCCACTGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGc                                                                                                                                                    >  1:208851/1‑251 (MQ=255)
                                                                             tccccacATAGCCACGGCCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CCCCATTAACACTGGCAATGCCGCCATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCAAGCCCCCCCATATTGGGCATGTTCCATTCGGAGATAACAAACCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATGATTGAATCCCAGCTCtt                                                                                                                                                <  1:373735/246‑1 (MQ=255)
                                                                                    aTAGCCACTGGCCCCCGC‑TTGCGCCGCACCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCACCCGCACGAATTTTGTTCAGCAATTCCAGGCCATCCATATTGGGCATGTCCCAGCGGGAAAAAACAAATCCAAAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGt                                                                                                                                       <  1:181650/251‑1 (MQ=255)
                                                                                     gaGCCCCTGCCCCCCCC‑GTGCGCCGCAGCAATGATTTTCCCTTTCTTCGCTTCTGCATT‑CCCCATTACCACTGGCAATGCCGCCATCGGGCCGGCCGCACAAATGGTTTTCGGAAATTCCAGGCCACCCATATGGGGCAGGTTCCGGCGGGAGAAAACAAACCCAAAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTGTCAGCAGGtt                                                                                                                                      <  2:363931/250‑1 (MQ=255)
                                                                                       cccACTGGCCCCCGC‑TGGCGCGGCAGCAATGATGTTCTCTTTCTTCGCTTCGGCAGT‑CACCATGACCACTGGCAATGCCGACATCGCCCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGAAAACAAACCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTAc                                                                                                                                    <  1:186894/250‑1 (MQ=255)
                                                                                          aCTGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCTGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGc                                                                                                                                  >  2:134705/1‑250 (MQ=255)
                                                                                           cggccccccGC‑TTACGCCCAAGCAATGATTTTCTCTTTCTTCGCTTTCGCAGT‑CCCCATTAACACGGGCAATGCCAACATCGCGCCACCCGCACGAATTTTTTTCAGCAATTCCAGGCCACCCATATGGGGCATGTCCCAGTCGGAGATAACAAACCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATGATTGAATCCCAGCTCTTTCAGCAGGTTACGCAc                                                                                                                                <  2:80666/247‑1 (MQ=255)
                                                                                              gCCCCCGC‑TTGCGCCGCAGCAATGATGATCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTAt                                                                                                                             >  2:256395/1‑251 (MQ=255)
                                                                                                   cGC‑TTGCGCCGCAGCATTGATGTCCCCTTTCTTCGCTTCTACAGT‑CCCCATTAACCCGGGCAACGCCGACATCGCCCCACCCGCACAAATTGTTTGCACCAATGCCAGGCCATCCATATTGGGCATTTTCCAGCGGGAGATAACAAACCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCATGTTACGCACTATGCGTc                                                                                                                        <  1:136702/251‑1 (MQ=255)
                                                                                                              gcagcaATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCAtggtgg                                                                                                               >  1:281264/1‑250 (MQ=255)
                                                                                                                agcaATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACCCGGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGaaa                                                                                                            <  1:120222/251‑1 (MQ=255)
                                                                                                                    atgatgTTCTCTTTCTTCGCCTCTCCTGT‑CACCAATAACACTGGAAATGCCGACATCGCGCCATCCCCAAGAATTGTTTTCAGCAATCCCAGGCCATCCATATTGGGCATGTCCCAGTCGGAGATAACAAACCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTc                                                                                                        <  2:99475/251‑1 (MQ=255)
                                                                                                                       atgTTCTCTTCCTCCGCTTCTGCATT‑CACCATTAACACGGGCAATGCCGCCACCGCGCCATCCGCACGAATTGTTTGCAGCAATCCCAGGCCATCCATATGGGGCATGTTCCAGCCGGAGAAAACAAATCCAAAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGGCGCATGGTGGAAAAGTCATc                                                                                                     <  1:151636/251‑1 (MQ=255)
                                                                                                                           cctctTTCTTCGCTTCTG‑AGTTCACCATTACCCCTGCCAAGCCCGACACCGCCCCATCCGCCCGAATTGTTTTCAACAATTCCAGGCCATCCATATTGGGCAGTTTCCAGTCGGAGATAACAAACCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCcac                                                                                                  <  2:407590/249‑1 (MQ=255)
                                                                                                                                ttcttcGCTTCTGCAGT‑CACCATTAACCCGGTCAATGCGGACATCGCGCCTTCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGGCGGAGATAACAAATCCAGAACC‑GCCTGCCTGCAACGTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTGCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAAc                                                                                               <  1:434683/248‑1 (MQ=255)
                                                                                                                                ttctccgcttctccctt‑CACCATTAACACGGGAAAAGCCGACATCGCCCCATCCGCACGAATTGTTTCCATCAATTCCAGGCCATCCATTTGGGGCATGTTCCAGTCGGAGATACCAAACCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCa                                                                                             <  2:160576/234‑1 (MQ=255)
                                                                                                                                    tcGCTTCTGCATT‑CACCATTAACACTGGCAATGCTGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAACCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCaaa                                                                                           <  1:134705/248‑1 (MQ=255)
                                                                                                                                    tcGCTTCTGCAGT‑CCCCATAAACCCGGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGAAATTCCAGCCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCaaa                                                                                           <  1:434110/248‑1 (MQ=255)
                                                                                                                                    tcGCTTCTGCAGT‑CACCATTAACACTGCCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCTTTTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCaaaa                                                                                          <  1:331494/249‑1 (MQ=255)
                                                                                                                                             gcagccacaattaaaccTGGCAAAGCCGAAATCGCGCCATCCGCCCGAATTGTTTTCAGCAATTCCGGCCCCACCATATGGGGCATGTTCCAGTGGGAGATAACAAATCAATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTc                                                                                <  1:251254/235‑1 (MQ=255)
                                                                                                                                               gT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCttt                                                                             >  2:69592/1‑251 (MQ=255)
                                                                                                                                                           cacTGGCAATGCCGCCACCGCTCCATCCGCACGAATTGTTGGCAGAAATTCCAGGCTATCCATTTTGGGCATGTTCCAGTCGGAGATAACAAATCCATATCC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTTGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCAttt                                                                  <  1:409321/251‑1 (MQ=255)
                                                                                                                                                                    cgccgccATCGCCCCACCCGCACGATTTGTTTCAAACAATCCCAGCCCACCCATATTGGGCATTTTCCAGTCGGAGATAACAAACCAATACCC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCt                                                         <  1:4922/250‑1 (MQ=255)
                                                                                                                                                                        gACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCTCCATTTCACACTCCTgctt                                                     >  2:111279/1‑248 (MQ=255)
                                                                                                                                                                           aTCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAACAATTTAAGTTCTTTATCCGCCCTTTCACACTCCTGATTTa                                                   >  2:299245/1‑250 (MQ=255)
                                                                                                                                                                                gcCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACg                                             >  1:356711/1‑251 (MQ=255)
                                                                                                                                                                                             aTTGTTTCCAGAAATCCCGGCCCACCCATATTGGGCATGTTCCAGTCGGAGAAAACAAACCCAAACCC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCg                                <  1:69592/251‑1 (MQ=255)
                                                                                                                                                                                              tttttttcagaaatctcgggCCATCCATTTTGGGCATTTTCCTGCCGGAAAAAAAAAATGCAAAACC‑TCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGTTTCCTCAACAGTATTGAAGCCCAGCGCTGTCAGCAGGTTACTCACTATGCGTCGCATGGTGGAAAAGTCATCCACACCCAAAAATTTAAGTTCTTTATCCGCCATTTCACACGCCTGATTTAAATACGTATCGCCTGTCCg                                <  1:386792/232‑1 (MQ=255)
                                                                                                                                                                                                     cAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAAt                        >  1:139474/1‑251 (MQ=255)
                                                                                                                                                                                                        cAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCcac                                                                                                  >  2:276696/1‑174 (MQ=255)
                                                                                                                                                                                                        cAATCCCAGGCCACCCATATGGGGCATGTTCCAGTCGGAGATAACAAACCCAAAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCcac                                                                                                  <  1:276696/174‑1 (MQ=255)
                                                                                                                                                                                                                ggCCATCCATATTGGGAATGTTCCCGTCGGAGATAACAAATACATAACC‑GCCTGACTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACCCTCCTGCTTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAAc              >  2:206299/1‑250 (MQ=255)
                                                                                                                                                                                                                      ccatATTGGGCATGTTCCAGTCGGAGATAACACATCCACAACC‑GCCTGCCTGCAACTTATTGAGAGCGCCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCt       >  2:204789/1‑251 (MQ=255)
                                                                                                                                                                                                                        atgtgggatttttcccAGTCGGAGA‑AAAAAAACCAAAACCCGCCTGCCTGCAACTTATTGAGGGCGCCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTGTGCCAACATTTGCTgg     <  2:289227/237‑1 (MQ=255)
                                                                                                                                                                                                                           ttGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCtt  >  2:99005/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                              
ATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGC‑TTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGT‑CACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACC‑GCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTT  >  NC_000913/1967036‑1967503

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