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→GGGcheYchemotaxis 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.

ATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTT  >  NC_000913/1967036‑1967503
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                              
atAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCAt                                                                                                                                                                                                                            >  1:345051/1‑251 (MQ=255)
          ccTCACATGCCCAGTTTCTCCAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCgcctgc                                                                                                                                                                                                                 >  1:91873/1‑251 (MQ=255)
            tCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAAAACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCgcctgcct                                                                                                                                                                                                               >  1:307919/1‑251 (MQ=255)
                   cccAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAACGTCGCGGCGGTAAATGGCTTCAACACATAGACACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCCCCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACtt                                                                                                                                                                                                        >  1:306145/1‑251 (MQ=255)
                    ccAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACtt                                                                                                                                                                                                        >  1:411191/1‑250 (MQ=255)
                     cAGTTTCTAAAAGATTTGTTTGAGTTTTTCCCCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACGGGCCCCCGCTTGCGCCGCAGCAATGATGTCCTCTTTCTTCGCTTCTGCAGTCACCATAACCACTGGCAATCCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGCCATGTTCCAGTCGGAGATAACAAATCCAGAACCGCCGGCCTGCAACTTAt                                                                                                                                                                                                      <  1:361394/251‑1 (MQ=255)
                        tttCTCAAAGATTCTCTTGCGTTTTTACTCCAGCGTCGCGGCGGTCAATGCCTTAACCAAATAGCCACTGGACCCCGCCTTGACCGCCGCAATGTCGTTCTCTTTCTTCGCGTCTGCAGTCAACATTCACACTGGCAATGCCCACATCTCGCCCTCCGCACGAATTGTTTTCAGCAATTCCAGGCCCTCCATAATCGGCATGTTCCCGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTCTTg                                                                                                                                                                                                    >  2:404905/1‑250 (MQ=255)
                              aaaGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGGGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAAGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCATCAACTTATTGAGAGCgg                                                                                                                                                                                             >  1:160576/1‑250 (MQ=255)
                                  aTTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGAc                                                                                                                                                                                         >  1:374234/1‑251 (MQ=255)
                                    ttttttgAGTTTTCCCCCCCGCGCCGCGGCGGTAAAGGGCTACCCCACATAGCCCCGGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTCCGCTTCTGCAGTCCCCATTAACCCTGGCAATGCCGACATCGCGCCACCCCCACAATTTGTTTTCAGCAATCCCAGGCCATCCATATGGGGCATGTTCCAGTCGGAGATAACAAATCCAAAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACg                                                                                                                                                                                        <  2:377686/246‑1 (MQ=255)
                                              ttttCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAACCACTGGCCCCCGCTTGCGCCGCATCAATGATGTTCTCTTTCTTCGCTTCTGCAGACACCATTCACACTGGCACTGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCCATTCCAGGCCATCCATATTGGGCATGTTCCAGTAGGAGATAACACATCCATAACCGCCTGCCTGCCACTTATTGCGAGCGTCGACGcccacctccgc                                                                                                                                                                             >  1:32235/1‑242 (MQ=255)
                                                  cctccAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTcc                                                                                                                                                                         >  2:331494/1‑251 (MQ=255)
                                                       aGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCCTCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCCCTTCCTCacc                                                                                                                                                                    >  2:4922/1‑249 (MQ=255)
                                                             gcggcgGTAAAGGGCTTCCCCCAATAGCCCCTGCCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACGGGAAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACattatt                                                                                                                                                              <  1:29163/251‑1 (MQ=255)
                                                                gcgatAAATGCCTTCACCACATAGCCTCTTGCCACAGATTTCCCCGCAGCAAGTATGTTCTTTTTCTTCGCTTCCGCAGTCACCATTAACACGGCCAATCCCCACACCGCGCAATCCGCACGAATTGTTTTCAGCAATACCAGGCAATCCACATTGGGCATGTTCAAGTCGGAGTTACCAAATCCAAAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCGCAACATAATTGaa                                                                                                                                                           <  1:88829/247‑1 (MQ=255)
                                                                     aaagggCTTCACCACATAGCCCCTGGCCCCCGCTTGCCCCGCAGCAATGATGTCCTCTTTCTCCGCTCCTGCAGTCCCCATTAACCCTGGCAATCCCGACATCGCCCCATCCGCACGAATTGTTTTCAGCAATCCCAGCCCATCCATATTGGGCATGTTCCAGTCGGAGAAAACAAACCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCa                                                                                                                                                      <  1:265513/247‑1 (MQ=255)
                                                                       aTGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGc                                                                                                                                                    >  1:208851/1‑251 (MQ=255)
                                                                             tccccacATAGCCACGGCCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCCCCATTAACACTGGCAATGCCGCCATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCAAGCCCCCCCATATTGGGCATGTTCCATTCGGAGATAACAAACCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATGATTGAATCCCAGCTCtt                                                                                                                                                <  1:373735/246‑1 (MQ=255)
                                                                                    aTAGCCACTGGCCCCCGCTTGCGCCGCACCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCACCCGCACGAATTTTGTTCAGCAATTCCAGGCCATCCATATTGGGCATGTCCCAGCGGGAAAAAACAAATCCAAAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGt                                                                                                                                       <  1:181650/251‑1 (MQ=255)
                                                                                     gaGCCCCTGCCCCCCCCGTGCGCCGCAGCAATGATTTTCCCTTTCTTCGCTTCTGCATTCCCCATTACCACTGGCAATGCCGCCATCGGGCCGGCCGCACAAATGGTTTTCGGAAATTCCAGGCCACCCATATGGGGCAGGTTCCGGCGGGAGAAAACAAACCCAAAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTGTCAGCAGGtt                                                                                                                                      <  2:363931/250‑1 (MQ=255)
                                                                                       cccACTGGCCCCCGCTGGCGCGGCAGCAATGATGTTCTCTTTCTTCGCTTCGGCAGTCACCATGACCACTGGCAATGCCGACATCGCCCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGAAAACAAACCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTAc                                                                                                                                    <  1:186894/250‑1 (MQ=255)
                                                                                          aCTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCTGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGc                                                                                                                                  >  2:134705/1‑250 (MQ=255)
                                                                                           cggccccccGCTTACGCCCAAGCAATGATTTTCTCTTTCTTCGCTTTCGCAGTCCCCATTAACACGGGCAATGCCAACATCGCGCCACCCGCACGAATTTTTTTCAGCAATTCCAGGCCACCCATATGGGGCATGTCCCAGTCGGAGATAACAAACCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATGATTGAATCCCAGCTCTTTCAGCAGGTTACGCAc                                                                                                                                <  2:80666/247‑1 (MQ=255)
                                                                                              gCCCCCGCTTGCGCCGCAGCAATGATGATCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTAt                                                                                                                             >  2:256395/1‑251 (MQ=255)
                                                                                                   cGCTTGCGCCGCAGCATTGATGTCCCCTTTCTTCGCTTCTACAGTCCCCATTAACCCGGGCAACGCCGACATCGCCCCACCCGCACAAATTGTTTGCACCAATGCCAGGCCATCCATATTGGGCATTTTCCAGCGGGAGATAACAAACCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCATGTTACGCACTATGCGTc                                                                                                                        <  1:136702/251‑1 (MQ=255)
                                                                                                              gcagcaATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCAtggtgg                                                                                                               >  1:281264/1‑250 (MQ=255)
                                                                                                                agcaATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACCCGGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGaaa                                                                                                            <  1:120222/251‑1 (MQ=255)
                                                                                                                    atgatgTTCTCTTTCTTCGCCTCTCCTGTCACCAATAACACTGGAAATGCCGACATCGCGCCATCCCCAAGAATTGTTTTCAGCAATCCCAGGCCATCCATATTGGGCATGTCCCAGTCGGAGATAACAAACCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTc                                                                                                        <  2:99475/251‑1 (MQ=255)
                                                                                                                       atgTTCTCTTCCTCCGCTTCTGCATTCACCATTAACACGGGCAATGCCGCCACCGCGCCATCCGCACGAATTGTTTGCAGCAATCCCAGGCCATCCATATGGGGCATGTTCCAGCCGGAGAAAACAAATCCAAAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGGCGCATGGTGGAAAAGTCATc                                                                                                     <  1:151636/251‑1 (MQ=255)
                                                                                                                           cctctTTCTTCGCTTCTGAGTTCACCATTACCCCTGCCAAGCCCGACACCGCCCCATCCGCCCGAATTGTTTTCAACAATTCCAGGCCATCCATATTGGGCAGTTTCCAGTCGGAGATAACAAACCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCcac                                                                                                  <  2:407590/249‑1 (MQ=255)
                                                                                                                                ttcttcGCTTCTGCAGTCACCATTAACCCGGTCAATGCGGACATCGCGCCTTCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGGCGGAGATAACAAATCCAGAACCGCCTGCCTGCAACGTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTGCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAAc                                                                                               <  1:434683/248‑1 (MQ=255)
                                                                                                                                ttctccgcttctcccttCACCATTAACACGGGAAAAGCCGACATCGCCCCATCCGCACGAATTGTTTCCATCAATTCCAGGCCATCCATTTGGGGCATGTTCCAGTCGGAGATACCAAACCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCa                                                                                             <  2:160576/234‑1 (MQ=255)
                                                                                                                                    tcGCTTCTGCATTCACCATTAACACTGGCAATGCTGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAACCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCaaa                                                                                           <  1:134705/248‑1 (MQ=255)
                                                                                                                                    tcGCTTCTGCAGTCCCCATAAACCCGGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGAAATTCCAGCCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCaaa                                                                                           <  1:434110/248‑1 (MQ=255)
                                                                                                                                    tcGCTTCTGCAGTCACCATTAACACTGCCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCTTTTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCaaaa                                                                                          <  1:331494/249‑1 (MQ=255)
                                                                                                                                             gcagccacaattaaaccTGGCAAAGCCGAAATCGCGCCATCCGCCCGAATTGTTTTCAGCAATTCCGGCCCCACCATATGGGGCATGTTCCAGTGGGAGATAACAAATCAATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTc                                                                                <  1:251254/235‑1 (MQ=255)
                                                                                                                                               gTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCttt                                                                             >  2:69592/1‑251 (MQ=255)
                                                                                                                                                           cacTGGCAATGCCGCCACCGCTCCATCCGCACGAATTGTTGGCAGAAATTCCAGGCTATCCATTTTGGGCATGTTCCAGTCGGAGATAACAAATCCATATCCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTTGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCAttt                                                                  <  1:409321/251‑1 (MQ=255)
                                                                                                                                                                    cgccgccATCGCCCCACCCGCACGATTTGTTTCAAACAATCCCAGCCCACCCATATTGGGCATTTTCCAGTCGGAGATAACAAACCAATACCCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCt                                                         <  1:4922/250‑1 (MQ=255)
                                                                                                                                                                        gACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCTCCATTTCACACTCCTgctt                                                     >  2:111279/1‑248 (MQ=255)
                                                                                                                                                                           aTCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAACAATTTAAGTTCTTTATCCGCCCTTTCACACTCCTGATTTa                                                   >  2:299245/1‑250 (MQ=255)
                                                                                                                                                                                gcCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACg                                             >  1:356711/1‑251 (MQ=255)
                                                                                                                                                                                             aTTGTTTCCAGAAATCCCGGCCCACCCATATTGGGCATGTTCCAGTCGGAGAAAACAAACCCAAACCCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCg                                <  1:69592/251‑1 (MQ=255)
                                                                                                                                                                                              tttttttcagaaatctcgggCCATCCATTTTGGGCATTTTCCTGCCGGAAAAAAAAAATGCAAAACCTCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGTTTCCTCAACAGTATTGAAGCCCAGCGCTGTCAGCAGGTTACTCACTATGCGTCGCATGGTGGAAAAGTCATCCACACCCAAAAATTTAAGTTCTTTATCCGCCATTTCACACGCCTGATTTAAATACGTATCGCCTGTCCg                                <  1:386792/232‑1 (MQ=255)
                                                                                                                                                                                                     cAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAAt                        >  1:139474/1‑251 (MQ=255)
                                                                                                                                                                                                        cAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCcac                                                                                                  >  2:276696/1‑174 (MQ=255)
                                                                                                                                                                                                        cAATCCCAGGCCACCCATATGGGGCATGTTCCAGTCGGAGATAACAAACCCAAAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCcac                                                                                                  <  1:276696/174‑1 (MQ=255)
                                                                                                                                                                                                                ggCCATCCATATTGGGAATGTTCCCGTCGGAGATAACAAATACATAACCGCCTGACTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACCCTCCTGCTTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAAc              >  2:206299/1‑250 (MQ=255)
                                                                                                                                                                                                                      ccatATTGGGCATGTTCCAGTCGGAGATAACACATCCACAACCGCCTGCCTGCAACTTATTGAGAGCGCCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCt       >  2:204789/1‑251 (MQ=255)
                                                                                                                                                                                                                        atgtgggatttttcccAGTCGGAGAAAAAAAACCAAAACCCGCCTGCCTGCAACTTATTGAGGGCGCCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTGTGCCAACATTTGCTgg     <  2:289227/237‑1 (MQ=255)
                                                                                                                                                                                                                           ttGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCtt  >  2:99005/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                              
ATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTT  >  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: