Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,1050AG31.1% 65.0 / 12.9 45I441V (ATC→GTC) mplUDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase
Reads supporting (aligned to +/- strand):  ref base A (4/27);  new base G (14/0);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 1.83e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑CGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTA‑CCCTTCGTTCTCCGCTAACTCGCGCAGATAC  >  NC_000913/4456879‑4457347
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
gAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATgggg                                                                                                                                                                                                                                  >  1:390341/1‑250 (MQ=255)
gAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGc                                                                                                                                                                                                                                  <  2:146801/251‑1 (MQ=255)
 aaaaTGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGGCCATCCGAAAATGCTgggaggggg                                                                                                                                                                                                                                 >  2:20358/1‑244 (MQ=255)
    aTGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGCTggg                                                                                                                                                                                                                              >  1:401358/1‑250 (MQ=255)
             tGCAAAGAAGATCTGGCACCTTCATTAGGTCGTGCCCATGAAGTCTTACTGCTGCAACCGGCGCATATCCCGTGGCAGGTCGCCGCAGTGGCAGAAGACTGCGTTCAGCCTGCCCACTGGAGTGGCGCTGTGGATACGCTGGCCGCTATGGTGGTGAAAAACGCCCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGCGGGATCCCTCCGAAACTGCTGGCTGGGCTGGCggagaacgc                                                                                                                                                                                                                     >  1:34113/1‑248 (MQ=255)
              gCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGc                                                                                                                                                                                                                    <  2:448629/250‑1 (MQ=255)
              gCAAAGACGATCTGGCACCTTCATTAGGTCGGGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTTGCgaaga                                                                                                                                                                                                                         <  1:358673/246‑1 (MQ=255)
              gCAAAGACGATCTGGCACCTTCAGTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATTAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATTGGCTGGCGAAGAA‑GGCg                                                                                                                                                                                                                   <  1:392996/251‑1 (MQ=255)
                    aCGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCGGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGTCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGACTGGCGAAGAA‑GGCGGAAGcc                                                                                                                                                                                                             <  1:161060/251‑1 (MQ=255)
                       aTCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCTCAGAAACTGCTGGGTGGGGGGGGGAAGAA‑GGCGGGAGCcgcg                                                                                                                                                                                                          >  1:291691/1‑251 (MQ=255)
                           ggCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCa                                                                                                                                                                                                        <  2:118566/249‑1 (MQ=255)
                               ccTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTTCCTCAGAAACTGCTGGGTGGGCTGGCGGAGAA‑GGCGGGAGCCGCGCAGTAAtt                                                                                                                                                                                                  >  1:227453/1‑251 (MQ=255)
                               ccTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAAtt                                                                                                                                                                                                  <  1:205293/251‑1 (MQ=255)
                                  tCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCg                                                                                                                                                                                                <  2:167710/250‑1 (MQ=255)
                                       aGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCtt                                                                                                                                                                                          <  1:104950/250‑1 (MQ=255)
                                        ggTCGTGCCGCTGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCCGAAACTGCTGGATGGGCCGGGGGAGAA‑GGCGGAAGCCGCGCAGGAAtttgggctt                                                                                                                                                                                           >  1:384317/1‑242 (MQ=255)
                                                   tGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGGGGGGTTCCTCAGAAAACGCTGGGTGGGGTGGGGGAGAA‑GGGGGGAGCCGGCCAGTAATTAGG‑GCTTA‑GGCAGataatat                                                                                                                                                                            >  1:460856/1‑247 (MQ=255)
                                                    gAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGt                                                                                                                                                                            <  2:384317/250‑1 (MQ=255)
                                                              tgctgcAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGATCCCTCAGAAAATGCTGGATGGGGTGGCGAAGAA‑GGCGGGAGCCGCGCAGTAATTCGG‑CCTTA‑GCCCGATATAGGATTCGTAGGcc                                                                                                                                                                 >  1:142279/1‑251 (MQ=255)
                                                                     aCCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAAATCGG‑CCTTA‑GCCCGATATAGTATGCGTAGGCCGGATa                                                                                                                                                            >  1:368170/1‑249 (MQ=255)
                                                                      ccGGCGCATCTTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTTCACCCGAAAACGCTGGGTGGGGTGGGGGAGAA‑GGCGGGAGACGCGGAGGTATTCGG‑CCTTAAGCCGGA‑ATAGTTTGCGTTGGCCGGGTaa                                                                                                                                                           >  1:84796/1‑249 (MQ=255)
                                                                         gcgcATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGGAGAA‑AGCGGAAGCCGCGCAGGAATTCGG‑CCTTT‑GCCCGATATAGTATGCGTAAGCCGGATAAGGCgg                                                                                                                                                      >  1:171240/1‑250 (MQ=255)
                                                                             atatTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAc                                                                                                                                                  <  1:430838/251‑1 (MQ=255)
                                                                               atTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                                <  1:327530/251‑1 (MQ=255)
                                                                                     tggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGAGCCATCAGAAACGGCTGGATGGGCTg                                                                                                                                                                                                                                <  2:467705/168‑1 (MQ=255)
                                                                                         aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCCCCAAAACTGCTGGGTGGGCTGGCGAAGAAAGGGGGAAGC‑GCGCCGGAAATCGGGCCTT‑‑GCCCGGAATAGGATGCGTAAGCCGGATAAGGccgtccccccccctcccgc                                                                                                                                      >  2:79508/1‑233 (MQ=255)
                                                                                           gtggcagaagtggcagaaGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCTCAGAAACTGCTGGATGGGCTGGCGGAGAA‑AGCGGCAACCGCGCCGTAATTCGG‑CCCTA‑GCCCGATATAGTATGCGTAGGCCGGATAAGGCGTTTCCGCCGCATCCGGCAt                                                                                                                                    >  1:70700/1‑251 (MQ=255)
                                                                                              gcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGTCGTTCACGCCGCATCCGGCATTTg                                                                                                                                 <  1:394906/251‑1 (MQ=255)
                                                                                              gcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGGGTTCACGCCGCATCCGGCATTTg                                                                                                                                 <  1:166576/251‑1 (MQ=255)
                                                                                                    gtggcagaagCCTGCGTTCAGCCGGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGGACATCAGAAACAGCTTGAAGGGGGGGCGGAGAA‑GGCGGAGACCGCGCAGAAATTCGGGCCTT‑‑GGCAGCTCTAGGAAGGGTAGGCCGGAAAAGGCGGTGAcgccgccacaggcctttgagcct                                                                                                                            >  2:270065/1‑233 (MQ=255)
                                                                                                                            gCACACGGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑Cg                                                                                                 <  2:429584/251‑1 (MQ=255)
                                                                                                                            gCACACGGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑Cg                                                                                                 <  2:266466/251‑1 (MQ=255)
                                                                                                                                        ggCGATGTGGATACGCTGGCAGATATGGGGGTGAAAACCGCTCAGCCTGGCGACCATAGTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑CGTCTTATCAGGcc                                                                                     <  2:242774/251‑1 (MQ=255)
                                                                                                                                           gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑CGTCTTATCAGGCCTag                                                                                  <  2:340063/251‑2 (MQ=255)
                                                                                                                                           gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑CGTCTTATCAGGCCTAc                                                                                  <  1:468540/251‑1 (MQ=255)
                                                                                                                                                    acgatGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑CGTCTTATCAGGCCTACAGGTGCAcc                                                                         <  2:242456/247‑1 (MQ=255)
                                                                                                                                                     cGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑CGTCTTATCAGGCCTACAGGTGCACCg                                                                        <  2:165538/251‑1 (MQ=255)
                                                                                                                                                       cTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGGAGGCCGGGTAAGGCGTTCACGCCGCATCCGGCATTTTAGCATAGTGCCTGATGCG‑ACGCTTGATG‑CGTCTTATCAGGCCTACAGGGGCACCGt                                                                       >  1:127787/1‑250 (MQ=255)
                                                                                                                                                            aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGAGGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑CGTCTTATCAGGCCTACAGGTGCACCGTATCCgg                                                                 <  2:248534/251‑1 (MQ=255)
                                                                                                                                                                        gAAAACCGCTAAGCCTGGAGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGCGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAAACGCGCAGTAATTCGG‑ACTTA‑ACCCGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCATCCGGCATTTGGGCATAGTGCCTGATGCG‑ACGCCTGATT‑CGTCTTATCAGGCCTACAGGGGCACCGAATCCGGCAAAACATACTc                                                     >  1:87521/1‑251 (MQ=255)
                                                                                                                                                                            aCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGTTATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATATTGCCTGATGCG‑ACGGTTGAGG‑CGTGTTATCAGGCGTACAGGTGCACCGTATCCGGCAAACCATACTCACAt                                                 <  1:279811/251‑1 (MQ=255)
                                                                                                                                                                                        ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑CGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGata                                      <  1:77921/250‑3 (MQ=255)
                                                                                                                                                                                                  ttCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑CGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTA‑CCCTTc                           <  2:356305/250‑1 (MQ=255)
                                                                                                                                                                                                        tgatgaGCAACGGCGGTTTTGGGGGGATCCCTCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGAGCAGGAATTCGG‑CCTTA‑GCCAGCTATAGTAAGCGTAGGCCGGATAAGGGGTTCCCGCCGCATCCCGCATTTGAGCATAATGCCTGATGCGCAC‑CCTGATGACG‑CCTATCAGGCCCAACCGTGCACCGCATCCGGCAAACCCTACCCACCTCAACAACGAAAATTAACCCTTcttcttc                     >  2:400073/1‑245 (MQ=255)
                                                                                                                                                                                                                           tggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑CGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTA‑CCCTTCGTTCTCCGCTAACTCGCGCAGATAc  <  1:148676/250‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
GAAAATGGGGATCTGCAAAGACGATCTGGCACCTTCATTAGGTCGTGCCGATGAAGTCTTCCTGCTGCAACCGGCGCATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAA‑GGCGGAAGCCGCGCAGTAATTCGG‑CCTTA‑GCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCG‑ACGCTTGATG‑CGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTA‑CCCTTCGTTCTCCGCTAACTCGCGCAGATAC  >  NC_000913/4456879‑4457347

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