Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,548,0330TG29.5% 129.6 / 11.5 61V454G (GTG→GGG) murPN‑acetylmuramic acid permease, EIIBC component, PTS system
Reads supporting (aligned to +/- strand):  ref base T (9/34);  new base G (17/1);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 8.45e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCCAGAGAAAGCAGGGTTTAACGTCGAACGGCTTAACCAGATGGATCGCTGGATTAGCCAGCAAGTTGATGTCGGTTATCCC  >  NC_000913/2547799‑2548269
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
tgGGCGGTGGTGAACCGCTGATTTACGGTGTCACCCGGCCCCGAATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCa                                                                                                                                                                                                                                            <  1:1083287/239‑1 (MQ=255)
             cccGCTGATTTACGGTGTCCCCCTGCCCCGCATGAAGCCGTTTGTTCCCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCggg                                                                                                                                                                                                                              <  1:1288129/237‑1 (MQ=255)
                 cTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCCTGGCGGCTTATGCTTGTGGGGTGTTTGGGGGCATGGGTTGCGGGGTTAt                                                                                                                                                                                                                         >  2:1243790/1‑239 (MQ=255)
                 cTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGGTTTTGCTGGCGGGGTTTTGGTGGCCTGGGGTTGCGGGttttt                                                                                                                                                                                                                         >  1:488161/1‑237 (MQ=255)
                          gggtgcccccTGCCCCGAATGAGGCCGTTGGTTCCCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGGTAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTcaccac                                                                                                                                                                                                                <  2:320736/232‑1 (MQ=255)
                              tCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGGTTTTGGGGGGCTTGGGTTTCGGGGTTTTTTTTCCCACGc                                                                                                                                                                                                              >  1:320736/1‑237 (MQ=255)
                                      ccccGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCCTGGCGGTTTTTGCTTGCGGGGTTCTTGGGGCCTGGGGTTGCGGGGTTTTTTTCTCCCCCCTCtttt                                                                                                                                                                                                        >  1:157330/1‑234 (MQ=255)
                                      ccccGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCCTGGCGGGTTTTGCTGGCGGGGGTCTTGTGGGGTTGGGTTTTGGGGTTTTTTTTACCACGCTCcttt                                                                                                                                                                                                        >  2:73909/1‑234 (MQ=255)
                                           cgtggcgCCGTGTGTTCCCGCCTGGTTGGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGt                                                                                                                                                                                                <  2:1220382/232‑1 (MQ=255)
                                           cATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGGTTATGCTGGCGGGGTTCTGGGGGCATGGGGTTTGGGGGTTTTTTTTACCACGCCCTTTGGGTGCCgca                                                                                                                                                                                               >  1:1220382/1‑237 (MQ=255)
                                           cATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGGCCGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGGCTGGTAGCGCTGCCGCTGAGGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGGTTATGATGGGGGGGGTGTTGGGGGGATGGGGTTCGGGGGTTTTTGTCACCACGCCCTTTggtgg                                                                                                                                                                                                    >  2:280459/1‑232 (MQ=255)
                                                        gtgtcggCCTGTTGAGGCGTCGTCGCGGGCGTTTTGTTTATCTGTCTGATAGCGTGGTGGGGTCGACCGATGGGCTTAAACAGCGCGTTTGGGCCGGCTTGTCTGGTAGCGCTGCGGCTGATTGCTTCCGCACACGGCTTCCTTCCGGCCATGGCGATTTATGCTGGCGGGATTCTGTTGGCATGGGTTTGCGGGTTTATTTTCACCACCCTATTTGGCTGCCGTAACGTCAATCTGGa                                                                                                                                                                                  <  2:198382/233‑1 (MQ=255)
                                                         ttcccGCCTGTTTGGGGGGGGCGGGGGGGGGTTTTTTTATCGGCCTGGTAGCCTGGGGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGa                                                                                                                                                                                  <  2:1381426/235‑1 (MQ=255)
                                                           aCCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGGTTCTTGTGGGATGGGGTTGGGGGGTTTTTTTTACCCCGCCCCTTTGGTGCCCGAACGCCCATTTgg                                                                                                                                                                                   >  1:794739/1‑235 (MQ=255)
                                                                       ggcggcGCGGCGGGCGGTTTGTTTTTCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAAc                                                                                                                                                                   <  1:280459/239‑1 (MQ=255)
                                                                        ggggcggggcggGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAAc                                                                                                                                                                   <  1:210526/236‑1 (MQ=255)
                                                                            cgcggcgGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGAc                                                                                                                                                               <  1:73909/238‑1 (MQ=255)
                                                                             ggggcgggcgGTTTGTTTATCGGCCTGAGTGCCGGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTGTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACa                                                                                                                                                              <  1:800001/236‑1 (MQ=255)
                                                                             gcggggGGGGGTTTGTTTATCGCCTTGTTAGCCTGGTGGGGTCGGCCGAGGGGCTTAAACAGCGCCTTGGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACa                                                                                                                                                              <  2:794739/238‑1 (MQ=255)
                                                                                     cggTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGTTTTATTTTCACCACTCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTa                                                                                                                                                      <  1:862440/238‑1 (MQ=255)
                                                                                                        aTAGCCTGGTGGGGTCTGCCGATGGTCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTAtt                                                                                                                                                                                                                        <  1:1130658/153‑1 (MQ=255)
                                                                                                                  ggggTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGGAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGGTTGTTGGGGGATTGGGTTGGGGGGTTTTTTTTACCCCGc                                                                                                                                                                                                              >  2:1130658/1‑153 (MQ=255)
                                                                                                                        tGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGGTTCTGGGGGCATGGGGTTGCGGGTTTTTTTTCACCACGCTCTTTTGCTGCCCTAACGTCAATCTGGGCTGATATCATTAAAAGTGCAAAGCCCCTTCTTTTTCTTCTGGGTTGTAGTTTTAGGGttatgg                                                                                                                   >  2:916387/1‑233 (MQ=255)
                                                                                                                        tGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGa                                                                                                                                                                <  1:1005516/193‑1 (MQ=255)
                                                                                                                                 gCTTAAACAGCGCCTTTGGGCTGGCTGTTCTGGTCGCGCTGCCGCTGATGCCTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAgcg                                                                                                                   <  1:1348928/229‑1 (MQ=255)
                                                                                                                                 gCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACg                                                                                                                                                                                            <  1:1403888/156‑1 (MQ=255)
                                                                                                                                 gCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCa                                                                                                                                                                         <  1:1242618/175‑1 (MQ=255)
                                                                                                                                 gCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATc                                                                                                                                                    <  1:1021445/196‑1 (MQ=255)
                                                                                                                                  cTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGGTTTTTGGGGCATGGGTTTGCGGGGTTATTTTTACCCCGCTCTTTGGTTGCCGGAACGTCCATCTTGGGTGATATTATGAAAAGGGCAATGctctt                                                                                                                                                      >  2:1005516/1‑192 (MQ=255)
                                                                                                                                  cTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCAt                                                                                                                                                                        <  1:1094812/175‑1 (MQ=255)
                                                                                                                                   ttAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATAcc                                                                                                       <  1:1243790/239‑1 (MQ=255)
                                                                                                                                    tAAACGGCGCCTTTGGGCCGTCTGGTCGGGTAGCGCTGCCGCTGATGACTGCCGCCCAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATAccc                                                                                                      <  1:405/239‑1 (MQ=255)
                                                                                                                                      aaCAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCAGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATg                                                                                                                                                           <  1:1002203/184‑1 (MQ=255)
                                                                                                                                           cgcCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGGTTCTGGGGGGATGGGGTTGCGGGGTTTTTTTTTCCCACCCCTTTGGCTGCCGTCACGTCCATTTGGGCTGAGATCATTAAACgggc                                                                                                                                                               >  2:1242618/1‑173 (MQ=255)
                                                                                                                                           cgcCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGGTTCTGGGGGGATGGGGTTGCGGGGGTATTTTTACCCCCCTCTTTGGGTGCCGTAACGGCAATCTggg                                                                                                                                                                                  >  2:1403888/1‑155 (MQ=255)
                                                                                                                                           cgcCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTTGTGGCATGGGTTTGCGGGGTTATTTTTACCACGCTCTTTGGCTGCCGGAACGTCAATCTGGGCTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTTGCGTCAGCGCCGCAAAAAACCCTgg                                                                                                   >  2:226372/1‑234 (MQ=255)
                                                                                                                                           cgcCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGAGTCTGGGGGAATGGGGTTGGGGGGTTATTTTCAACAAGCTCTTTGGCGGCCGGAACGTGAATATGGGTTGATATCATGTAACAGACAAAGCTCTTTCTTTCTCTGCGGGCCGTTAGTTGTAGGGTGAGAGCCGCCAAATa                                                                                                         >  2:1348928/1‑229 (MQ=255)
                                                                                                                                           cgcCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGGTTATGCTGGCGGGGTTCTTGGGGCATGGGGTTTCGGGGTTTTTTTTCCCACGCCCTTTTGGTGCCGGAACGGCAATTTGGGCTTGTATCATTAAACGGAAAAAGCTCTATCTTTTTCTTCTCGCTGTTTGTTTTAGAGTCAGCGCCGCGAAATATCCCGttttt                                                                                               >  2:1268039/1‑236 (MQ=255)
                                                                                                                                              cTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTcacc                                                                                                                                                                                                                  <  1:444130/121‑1 (MQ=255)
                                                                                                                                                        tctggtctggTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGGTTCTGGGGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCCAAATACCCCGTTCTGACAGAAAGCTCg                                                                                   >  1:1159856/1‑238 (MQ=255)
                                                                                                                                                        tctggtctggTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGGTTCTGGGGGCATGGGGTTGCGGGGTTTTTTTTACCACGCCCTTTg                                                                                                                                                                                                        >  2:444130/1‑121 (MQ=255)
                                                                                                                                                        tctggtctggTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGGGGCATGGGTTTGCGGGGTTATTTTTACCACGCTCCTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACACTGTTCTTACCGACAGCTCg                                                                                   >  2:142314/1‑238 (MQ=255)
                                                                                                                                                               tggtAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATAc                                                                                                        <  1:1341205/210‑1 (MQ=255)
                                                                                                                                                                 gtAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTAtttt                                                                                                                                                                                                                      <  1:914092/98‑1 (MQ=255)
                                                                                                                                                                   ggcggTGCCGCTGATGGCTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCCAGAGAAAGCa                                                                       <  1:1226451/238‑1 (MQ=255)
                                                                                                                                                                         gCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGGTTCTGGTGGCATGGGGTTGCGGGGTTATTTTTACCACGCCCTTTGGCTGCCGTAACGGCCATCTGGACTGATATCATTAAAAGGACAAAGCTCTCTCTTTTTCTTCTGGGTGTTAGTTTTAGCGTCAGCGCCGCCAAATACCCCGTTCTGa                                                                                              >  2:1341205/1‑210 (MQ=255)
                                                                                                                                                                           cGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGGTTTTGGGGGCATGGGGTTGGGGGTTTTTTTTTACCACGctc                                                                                                                                                                                                            >  2:914092/1‑98 (MQ=255)
                                                                                                                                                                                 tgaCTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCCAGAGAAAGCAGGGTTTAACGTCGa                                                         <  1:1268039/239‑1 (MQ=255)
                                                                                                                                                                                    cTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTTGTGGCATTGGGTTGCCGGTTTATTTTTACCAGCCGCGTTGGCCGACGTTACGTCCATTTGGACTGTTATCCTGAACAGGACAATGCTCCATCTTTCTCTGATGGCTGCTAGTATTAGCGTCAGCGCCGCAAAATACCCTGTTCTGTCAGAAAACTCTCCGAAAAAAACAGGGGTTAACgcgcaaac                                                      >  2:133667/1‑232 (MQ=255)
                                                                                                                                                                                        cGCACAAGGAATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATGCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCCAGAGAAAGCAGGGTTTAACGTCGAACGGCtt                                                  <  2:1321934/239‑1 (MQ=255)
                                                                                                                                                                                               ggCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGa                                                                                                                                                                                  <  2:823415/104‑1 (MQ=255)
                                                                                                                                                                                               ggCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGa                                                                                                                                                                <  2:479518/122‑1 (MQ=255)
                                                                                                                                                                                                  gtCCTTCCGGCCGTGGCGGTTTATGCTGGCGGGATTCTGGGGGCGGGGGTGTGCGGGTTTATTTTCCCCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCCAGAGAAAGCAGGGTTTAACGTCGAACGGCTTAACCAGATg                                         <  2:157330/237‑1 (MQ=255)
                                                                                                                                                                                                    ccTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCgggg                                                                                                                                                                                                                             >  1:366932/1‑55 (MQ=255)
                                                                                                                                                                                                    ccTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGt                                                                                                                                                                                                                             <  2:366932/56‑1 (MQ=255)
                                                                                                                                                                                                    ccTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCg                                                                                                                                                                                                 <  1:950036/84‑1 (MQ=255)
                                                                                                                                                                                                      ttCCGGCAAGGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCCAGAGAAAGCAGGGTTTAACGTCGAACGGCTTAACCAGATGGATCg                                    <  2:488161/239‑1 (MQ=255)
                                                                                                                                                                                                         cGGCCATGGCGGTTTATGCTGGCGGGGTTCTGGTGGCCTGGGGTTGCGGGTTTTTTTTCACCACGCTCTTTGGCTGCCGGAACGTCCATCTGGGCTGATATCAt                                                                                                                                                                        >  1:823415/1‑104 (MQ=255)
                                                                                                                                                                                                         cGGCCATGGCGGTTTATGCTGGCGGGATTTTGGTGGGATGGGGTTGGGGGGTTATTTTCACCACGCTCTTTGGGTGCCGTAAAGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTa                                                                                                                                                      >  1:479518/1‑122 (MQ=255)
                                                                                                                                                                                                              aTGGCGGTTTATGCTGGCGGGGTTCTGGGGGGATGGGGTTGCGGGGTTTTTTTTACCACGCTCTTTGGCGGCCGTAAAGTCAAt                                                                                                                                                                                       >  2:950036/1‑84 (MQ=255)
                                                                                                                                                                                                                                        cTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCCAGAGAAAGCAGGGTTTAACGTCGAACGGCTTAACCAGATGGATCGCTGGATTAGCCAGCAAGTTGATGTCGGTTATc    <  1:226372/239‑1 (MQ=255)
                                                                                                                                                                                                                                         tggtggCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCCAGAGAAAGCAGGGTTTAACGTCGAACGGCTTAACCAGATGGATCGCTGGATTAGCCAGCAAGTTGATGTCGGTTATcc   <  1:916387/239‑1 (MQ=255)
                                                                                                                                                                                                                                         tggtggCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCCAGAGAAAGCAGGGTTTAACGTCGAACGGCTTAACCAGATGGATCGCTGGATTAGCCAGCAAGTTGATGTCGGTTATcc   <  1:192350/239‑1 (MQ=255)
                                                                                                                                                                                                                                            gtggCATGGGTTTGCGGGGTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCCCGCCAGAGGAAACAGGGGTTAACGTCGAACGGCTTCACCAGAAGGATTGTTGGATTAGCCACGAAGTTGATGTCGGtttaccc  >  2:651197/1‑233 (MQ=255)
                                                                                                                                                                                                                                             tggCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCCAGAGAAAGCAGGGTTTAACGTCGAACGGCTTAACCAGATGGATCGCTGGATTAGCCa                        <  1:920795/215‑1 (MQ=255)
                                                                                                                                                                                                                                             tggCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCCAGAGAAAGCAGGGTTTAACGTCGAACGGCTTAACCAGATGGATCGCTGGATTAGCCa                        >  2:920795/1‑215 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
TGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCCCGCATGAAGCCGTTTGTTACCGCCTGTTTAGGCGGCGCGGCGGGCGGTTTGTTTATCGGCCTGATAGCCTGGTGGGGTCTGCCGATGGGCTTAAACAGCGCCTTTGGGCCGTCTGGTCTGGTAGCGCTGCCGCTGATGACTTCCGCACAAGGCATCCTTCCGGCCATGGCGGTTTATGCTGGCGGGATTCTGGTGGCATGGGTTTGCGGGTTTATTTTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCCAGAGAAAGCAGGGTTTAACGTCGAACGGCTTAACCAGATGGATCGCTGGATTAGCCAGCAAGTTGATGTCGGTTATCCC  >  NC_000913/2547799‑2548269

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