Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,592,2900TG25.5% 68.4 / 10.4 47F212V (TTC→GTC) tsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base T (8/27);  new base G (12/0);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 2.41e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCAC  >  NC_000913/4592041‑4592513
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                  
tccGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCTGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGt                                                                                                                                                                                                                                    <  1:450103/248‑1 (MQ=255)
         tATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑a                                                                                                                                                                                                                          <  1:89903/250‑1 (MQ=255)
            cacaATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCATCGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGTCTCCATGATATCGCCGTCAGCGAGAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGCGGGCGGGAGGATCGTCGGACTGGGGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGc                                                                                                                                                                                                                      <  1:21508/251‑1 (MQ=255)
              caATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGGAtttaaaccc                                                                                                                                                                                                                     >  1:197852/1‑244 (MQ=255)
               aaTGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTc                                                                                                                                                                                                                   <  2:347157/251‑1 (MQ=255)
               aaTGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTc                                                                                                                                                                                                                   <  1:191454/251‑1 (MQ=255)
                aTGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCg                                                                                                                                                                                                                  <  2:33122/251‑1 (MQ=255)
                  gcgcTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGc                                                                                                                                                                                                                 <  1:324854/250‑1 (MQ=255)
                                gATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAAt                                                                                                                                                                                                    <  1:257540/249‑1 (MQ=255)
                                      aCTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGGATCGcc                                                                                                                                                                                            <  1:121436/251‑1 (MQ=255)
                                           tAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGGCCTCTTTGCCGGCTGGGTCGGGATT‑TAAGACTCGCTGGGAGCGGCAAAGaaacggctggtt                                                                                                                                                                                       >  2:144586/1‑241 (MQ=255)
                                                  aGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCCTCTGGGTCGGTATT‑AAAGCCTCGCCGGTAGCGCCCATGTAACGCCTGGTTGGCAGcc                                                                                                                                                                                >  2:128718/1‑250 (MQ=255)
                                                           cAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCAGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTATGGGTCGGGATTTAAAGCC‑CGCTGGGAGGGACAATTAAGCGCCTGTTTGACAGGA‑ATCGTCaat                                                                                                                                                                      >  2:36844/1‑249 (MQ=255)
                                                               gAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCGTCTGGGTCGGGATT‑AAAACCCCGCTGGGAGCGCCCATGAAACGCCCGGTTGGCAGCA‑TTTGTCATATTGc                                                                                                                                                                  >  1:83412/1‑251 (MQ=255)
                                                                  ttCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGGCTGGGTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATAGCCTGCTTTAAAGCA‑ATCGTCATATTGGAgg                                                                                                                                                               >  1:256809/1‑251 (MQ=255)
                                                                     tttGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTTGGGATT‑AAAGCCTCGCTGGGAGCGCCCATGAATCGCCTGATTGGAAGCA‑TTCGTCATATTGCAggcgg                                                                                                                                                            >  1:259904/1‑251 (MQ=255)
                                                                           cAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATc                                                                                                                                                       <  1:144586/250‑1 (MQ=255)
                                                                             gCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTgg                                                                                                                                                    <  2:450641/251‑1 (MQ=255)
                                                                                     aGGGATATCAGGACGGTGGCGAGAAGCAGTATGTGGCGTACAGGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGTGTGATGAGCGTCGTACTCGCGGTCAGCTGCGCCTTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCTCCTGATTGGCAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCg                                                                                                                                            <  1:424736/251‑1 (MQ=255)
                                                                                        gatagtatgACGGTTTCGAGAAGCAGTATGTGTCTTACATGGAGCAAAACGATCGGCTCCATGATACCGGCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAGCCTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTTGTAGCGCCAATGAAACGCCTGATTGACAGCA‑ATCGTCATATTGCAGGCGGCGATCTGGTGAAACCGAtt                                                                                                                                         <  2:350386/243‑1 (MQ=255)
                                                                                           atCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGCGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATTTTCGCCGTCTGGTTCGGTATT‑AAAAACTCGCTGGTAGCGCCAATGAAACGCCTGATTGGCAGCA‑TTCGTCATATTGCCGGGGGCGATCTGGTGAAACCGATTggg                                                                                                                                      >  1:32221/1‑249 (MQ=255)
                                                                                                  cGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATg                                                                                                                               <  1:306623/251‑1 (MQ=255)
                                                                                                             aGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCAAATgag                                                                                                                    <  1:128718/251‑1 (MQ=255)
                                                                                                                aGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑Tg                                                                                                                <  2:83412/251‑1 (MQ=255)
                                                                                                                       ggCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCAAc                                                                                                         <  1:278456/251‑1 (MQ=255)
                                                                                                                            aCATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCAACTTCGCCGGCTGGTTCGGGATT‑TAAGACTCGCTGGGAGCGCCCAAGAAACGCCCGATTGACAGCA‑CTCGTTCTATTTCAGGCGGCGAACTTGGGGAACCGGTTGAGGGGGATggggccaaagagaggggggaaacggc                                                                                                      >  2:397062/1‑226 (MQ=255)
                                                                                                                              aTGGAGCCAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTCCGGATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTTCAGGCGGCGATATGGGGGAAACGCTTGAGGGGGGTGGCTCTAATGAGA‑TGGGGGAAcagggagg                                                                                                  >  1:175482/1‑244 (MQ=255)
                                                                                                                                 gAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCAACTGGCagagag                                                                                               <  1:55389/251‑1 (MQ=255)
                                                                                                                                  ggCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCAACTGGCAGAGAGt                                                                                              <  1:79959/250‑1 (MQ=255)
                                                                                                                                                      gATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TTGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGgag                                                                           <  1:415337/250‑1 (MQ=255)
                                                                                                                                                       atatCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGCTTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCCTCTTTGACGTCTGGTTCGGTATT‑AAAGCCTCGCTGGGAGCGCCAAAGAATCGCCTGATTGACA‑CACATCGTCATATTGCAGGCGGCGATCAGGTGCAACCGATTGAGGTGGATGGCTCTAATTAGA‑TGGGGGAAATGGCAGAGAGGTTGCGCCACATGCAGGgag                                                                           >  1:166402/1‑249 (MQ=255)
                                                                                                                                                         atCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCt                                                                       <  1:11265/251‑1 (MQ=255)
                                                                                                                                                         atCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCt                                                                       <  1:170892/251‑1 (MQ=255)
                                                                                                                                                             ccGTCAGCGATAACAATGCCTCCTACAGCCAGGAGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCAACTTTGCCGTCTGGGTCGGACTT‑CAAGCCTCGCTGCTAGCGCAAATGAATCGCCTGATTGACAGCA‑TTCGTCATAATGCAGGGGGCGGTCTGGTGGAACCGATTGAGGGGGCTGGCTCTAAAAAGACTGGGG‑TACAGGCAGAGCGATTGGGCCCCATGCAGGGAGAGCTGATg                                                                   >  1:145427/1‑251 (MQ=255)
                                                                                                                                                              cGTCAGCGATAACAATGCCTCCTCCAGCCAGGCGATGTGGCTTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTCGGTCGGGATT‑TAAGCCTCGCTGGTGAAGCCAATGGTACGCCGGATTGACAACA‑TTGGTCATATTGTCGGCGGCGATCGGGTGCAACCCATTGAGGGGGATGGCTCTAATGAAA‑TGGGGCAACTGGCAGAGAGATTGCGCCCTATGCCggggggactatagt                                                                  >  2:96005/1‑240 (MQ=255)
                                                                                                                                                              cGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGc                                                                  <  1:291303/251‑1 (MQ=255)
                                                                                                                                                                gcAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTTATGCGt                                                                <  2:197852/250‑1 (MQ=255)
                                                                                                                                                                  aGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTAc                                                              <  1:291247/251‑1 (MQ=255)
                                                                                                                                                                               cctcctACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGGAGCGCCAATGAATCGCCTGATTGACCGCC‑TTCGTCATATTGCAGGCGGCGGTCTGGTGAAACCGATTGAGGGGGATGGCTCTAATTAGA‑TGGGGCCAATGGCAGAGAGTTTGCGCCCTATGCAGGGAGAGGTGATGCGTACCGTCGGTGAtggg                                                 >  1:75661/1‑249 (MQ=255)
                                                                                                                                                                                 tcctACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCCCGCTGGTAGCGCCAATGAAACGCCTGAGTGACAGCA‑ATCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCCACTGGCAGAGAGTTTGCGCCATATGCCGGGGGAGCTGATGCGTACCGTCgggggtgggcg                                               >  1:323558/1‑242 (MQ=255)
                                                                                                                                                                                 tcctACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTCGGGATT‑AAAACCTCGCTGGGAGCGCCAATGAATCCCCTGGTTTACAGCA‑TTCGTCCTATTGCAGGCGGGGATCTGGGGAAACCGCTTTTGGGGGGTGGCGCTAAAGAGA‑AGGGGCAACTTGCAGAGGGTTTTGGCCCTATGCAGGGGGAGGGGATGCGGAcccgccgtggattggg                                               >  2:310740/1‑237 (MQ=255)
                                                                                                                                                                                     aCAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCGGCTGGGTAGGTATT‑AAAACCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCCGGCGGCGATCTGGTGAAACCGATTTAGGGGGGTGGCTCTAATGAGA‑AGGGGCAACTGGCAGAGAGTTTGCGCCCTATGCAGGGAGAGCTGGTGCGTACCGTCGGTGATGTGGGCAAc                                           >  1:357255/1‑251 (MQ=255)
                                                                                                                                                                                     aCAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGAt                                                                                                                                <  2:249274/167‑1 (MQ=255)
                                                                                                                                                                                          cAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGGATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCCTATGTCAGGCGGCGATCTGGTGACACCGAGTGGGGGGGATGGCTCTAATGAGA‑GGGGGCAACTGGCAGAGAGGTTGGGCCCTATGCCGGGGGAGCTGATGCGTACCGTCGGGGATGTGGGCAACggagg                                      >  2:89925/1‑248 (MQ=255)
                                                                                                                                                                                                   tgGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATATTTGACGTCTGGGTTGGTATT‑AAAGCCTCGCTTGTAGCGCCAATGAAACGCCTGATTGAACGAA‑TTCGTCATATTTCAGGCGGCGATCTGGTGAAACCGAATGAGGTGGATGGGTCTTATGAGA‑TGGGGCAACTGGCCGAGAGGTTTCGCCCTATGCCGGGAGAGCTGATGAGTACCGTCGGGGATGTGCGCAAAGGGGGCAAAGCCac                             >  2:273202/1‑250 (MQ=255)
                                                                                                                                                                                                               ggCGTGATGATCGTCGTACTGGCGGTTATCTTCGCCGTCTGGGTTGGTATT‑AAAGCCTCGCTGGGAGCGCCCATGAAACGCCCGGTTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGGGAAACCGATTGAGGTGGGTGGCTCTAATGAAG‑TGGGGCAAACGGCCGGGAGGTTGCGCCATATGCAGGGAGAGCTGAAGCGTACCGTCGGGGGTGTGGGCAACGGGGCCAATGCCATCTATAGCggggc                 >  1:72952/1‑248 (MQ=255)
                                                                                                                                                                                                                        agagtcgtACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAAGGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGTTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAAt        <  1:107208/248‑1 (MQ=255)
                                                                                                                                                                                                                              gtACTGGCGGTCATCTTCGCCGTCGGGTTCGTGATT‑AAAGCCTCGCTGGTAGCGACAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCAACTGGCAGAGAGTGTGCGCCATATGCGGGGAGAGCTGATGCGTACCGTCGGTGATGTTCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCAc  <  2:128045/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                  
TACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATT‑AAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCA‑TTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGA‑TGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCAC  >  NC_000913/4592041‑4592513

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