Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913119,3960CT23.1% 122.0 / 13.9 65S39F (TCT→TTT) ampEprotein AmpE
Reads supporting (aligned to +/- strand):  ref base C (13/37);  new base T (15/0);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.81e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCGGCAGGC  >  NC_000913/119159‑119646
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                             
tgCTATCCGGATATCGATAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAc                                                                                                                                                                                                                                                   <  2:446355/250‑1 (MQ=255)
  cgaTCCGGATATCGTTAAAAACATGACTGGCCATTGTGATATTGCGCCGGACCGGAAAACCGATCCCGGTCCTGCATTTGATTGGG‑ACGGGTTTCGTGTGCTGGTCAGCAATGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑Gtt                                                                                                                                                                                                                                               <  2:360893/249‑1 (MQ=255)
  cTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACAATTTTTTCTCgggggccagtt                                                                                                                                                                                                                                                >  1:38232/1‑243 (MQ=255)
     tCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACCTTTTTTTCTCgggggcccgggttg                                                                                                                                                                                                                                             >  2:181895/1‑240 (MQ=255)
      ccGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑Gttcgg                                                                                                                                                                                                                                            >  1:359864/1‑247 (MQ=255)
          atatCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTTTCTCCGGCGCCC‑GTTAGGC‑ATgg                                                                                                                                                                                                                                      >  1:214594/1‑250 (MQ=255)
          atatCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGGAAAAATTTTTTTTTCGGGGCCAC‑GTTAgcatc                                                                                                                                                                                                                                         >  2:483635/1‑246 (MQ=255)
                      aCATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑GAACATTTTTTTCTCGGGCGCCC‑GGTTGGC‑ATGATCCTTATTGCgg                                                                                                                                                                                                                          >  1:4944/1‑250 (MQ=255)
                                ccATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAAAATTTTTCTCTCGGGGGCAC‑GGTAAGC‑ATGACCATTATTGCGGTgggggggg                                                                                                                                                                                                                 >  1:190175/1‑245 (MQ=255)
                                    tgtgATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACttt                                                                                                                                                                                                            <  2:116172/250‑1 (MQ=255)
                                    tgtgATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACttt                                                                                                                                                                                                            <  1:466842/250‑1 (MQ=255)
                                        atatTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTAc                                                                                                                                                                                                       <  1:448310/251‑1 (MQ=255)
                                          atTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACt                                                                                                                                                                                                      <  1:65328/250‑1 (MQ=255)
                                            gtcgcCGGATCGGAAAACCGCTTCCGGTATTGCATTTGATTGGGCATTG‑TTGCGTGTGCTGGTCAGCACGGAGACTAC‑TGTCCGCTATGTACCACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCGTGAAGCTTGCTTTCGGCGGTTG‑AATCATGTTGCGCTCGGGCGCAC‑GTTAGGC‑TTGAGCATTATTGCGATGGGCTTG‑ACTTTTTTACTgc                                                                                                                                                                                                    <  2:266467/248‑2 (MQ=255)
                                              cgcCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACATTTTTCTCTCTGGCGCCC‑GTTAGGC‑ATGACCATTATTGCGAAGGGCGGG‑GCTTTTTTTCTGTTAc                                                                                                                                                                                                 >  1:356249/1‑251 (MQ=255)
                                                    aTCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑GAACATTTTTTTCTTCGGCGCCC‑GTTAGGC‑ATGACCCTTATTGCGATGGGCGGG‑GCTTTTTTTCTGTTACGCGCAt                                                                                                                                                                                           >  1:321675/1‑251 (MQ=255)
                                                     acGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCAtt                                                                                                                                                                                          <  2:148948/250‑1 (MQ=255)
                                                        gAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑GCTTTTTTTCTGTTTAGCGCATTGc                                                                                                                                                                                        >  1:135189/1‑250 (MQ=255)
                                                        gAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCa                                                                                                                                                                                       <  2:190175/251‑1 (MQ=255)
                                                             ccGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAg                                                                                                                                                                                  <  1:318673/251‑1 (MQ=255)
                                                               gATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTTGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTa                                                                                                                                                                                <  1:216983/251‑1 (MQ=255)
                                                               caTCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACGGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTa                                                                                                                                                                                <  2:483659/250‑1 (MQ=255)
                                                                      gTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACATTTTTTTCCCGGGCGCAC‑GGTAAGC‑CAAACCCATATTTCGCTGGGGGGG‑GCTTTTTTATTGTTACGCGCCTTGCAggggatattttttc                                                                                                                                                                         >  1:410052/1‑240 (MQ=255)
                                                                            cATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTTCGATGGGCGGG‑ACTTTTTTTTTTTTACGCGCCTTTCAGGGGGTATTGTTTAACGttt                                                                                                                                                                   >  1:372537/1‑250 (MQ=255)
                                                                              tttGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGGTGGGCGGG‑GCTTTTTTTCTGTTACGCGCCTTGCAGGGGGTATTGTTCAACGtttcc                                                                                                                                                                 >  1:295951/1‑248 (MQ=255)
                                                                                 gATTGGGCACGG‑TTTCGTGTGCTGGTCAGCNAGGAGACAACATG‑ACGATATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCAc                                                                                                                                                               <  2:234734/250‑1 (MQ=255)
                                                                                    tGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCt                                                                                                                                                            <  1:301376/250‑1 (MQ=255)
                                                                                     gggCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACAATTTTCTCTCCGGCGCCCCGTT‑GGC‑ATGACCCTTATTGTGATGGGGGCG‑CCTTTTTTTCTGTTACGCTCCTTTCCGGGGGGTTTGTTCACCATTCCCACGCTAc                                                                                                                                                          >  2:129126/1‑251 (MQ=255)
                                                                                         aCGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACCGGCAGCTTGATCATCGTCTTGAAGCGGTCTTTCGGCGGGTG‑AAACATTTTTCTCGCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGgt                                                                                                                                                      <  1:26561/251‑1 (MQ=255)
                                                                                                 cGTGTGCTGGTCAGCAAGGAGACCACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGgct                                                                                                                                                <  2:18641/250‑1 (MQ=255)
                                                                                                   ggtgCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGgctgct                                                                                                                                             <  1:262405/250‑1 (MQ=255)
                                                                                                    gtgCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGgctgct                                                                                                                                             <  1:65203/250‑1 (MQ=255)
                                                                                                    gtgCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGgctgct                                                                                                                                             <  1:198337/250‑1 (MQ=255)
                                                                                                         ggTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTgg                                                                                                                                       <  2:295951/251‑1 (MQ=255)
                                                                                                               cAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACATTTTTCTCTCGGGGGCCC‑GTTAGGC‑ATGACCCTTATTGCGGTGGGGGGG‑GCTTTTTTATTTTTACGCGCAGTGCAGGGGGGATTGTTTACCGTTCCCACCCTCCTGGGGGGGCTGCTGAtttgtttgtt                                                                                                                                 >  2:389320/1‑249 (MQ=255)
                                                                                                                   gagaCAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑GAACAATTTTCTCTCGGGCGCAC‑GGTAGGC‑CGGAACATTATTTCGGTGGGCGTGGATTTTTTTTCTGT‑ACGCGCCTTGCAGGGGGGATTGTTTTACGTTTCCACGATACTGgttgggggggtgatttgttttctgtgg                                                                                                                             >  1:182784/1‑226 (MQ=255)
                                                                                                                                      ttACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGt                                                                                                            <  1:138718/250‑1 (MQ=255)
                                                                                                                                       tACATCCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCAGCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATTGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGtt                                                                                                           <  2:38232/250‑1 (MQ=255)
                                                                                                                                        gcAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑GCAAATTTTTCTCTCCCGCGCCA‑GTTAGGC‑ATGGCCATTATTTCGATGGGCGTG‑TCTTTTTTACTGTTTCGCGCCTTTCAGGGGGGATTTTTTCACGTTCCCCCCCTTCTGGTGTGGCTGCTGATTGGGTTGGTGTGTATTGGCGCAggggtagtttc                                                                                                         >  2:265527/2‑242 (MQ=255)
                                                                                                                                        aCAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑GTGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTc                                                                                                          <  1:392407/250‑1 (MQ=255)
                                                                                                                                        aCAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAAAATTTTTTTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTTCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCCTTGCAGGGGGTATTGTTCAAAGTTTCCACGCTACTGGGGTGGCTGCTGGTTGTTTTTCTgtgtttttgggccgggaaaatttt                                                                                                         >  1:387116/1‑231 (MQ=255)
                                                                                                                                         cAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGAGTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGt                                                                                                        <  1:181895/251‑1 (MQ=255)
                                                                                                                                          aaCCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑GAACATTTTTTTCTCGGGCGCAC‑GTTAGGC‑ATGCCCATTATTGCGATGGGCGGG‑GCTTTTTTACTGTTACGCGCATTGCAGGGGGTATTGTTCAAACTTCCCACGCTACTGGGGGGGCTGCTGATTGGTTTGCTGTGTATTGGGGCAGGGAAAGttttgc                                                                                                       >  1:32025/1‑247 (MQ=255)
                                                                                                                                             cTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAAAATTTTTCTTTTGGGCGCAC‑GCTAGGC‑ATTTCCATTATTTGGGTGGGCGTG‑ACTTTTTTACTGTTACGCGCATTTCAGGGAGGATTGTTCCACGTTCCCCCGCTACTGGTGTGGGTGGTGGTTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTTGTCTTc                                                                                                    >  1:476437/1‑251 (MQ=255)
                                                                                                                                               tACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAAAATTTTTTTCTCGGGCGCAC‑GTTTGGC‑ATGACCCTTATTTCGATGGGGGTG‑ACTTTTTTTCTGTTACGCGCATTGCAGGGAGTATTGTTCTAAGTTCCCACGCTCCTGGTGGGGGTGCTGATTTGTTTGCTGTGTATTGGCGGCGGGAAAGTTCGTCTTcct                                                                                                  >  2:170849/1‑249 (MQ=255)
                                                                                                                                                actcgtggtAATTTTCGAGCGCCTGTTTAAGTTGGGCGATCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTTGGCGTGTG‑AATCATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACTATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCAt                                                                                                  <  2:478858/242‑1 (MQ=255)
                                                                                                                                                aCTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTTTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTTTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGGTTGGTTTGCTGTGTATTTGCGCAGGTAAAGTTCGTCTTCAtt                                                                                                 >  1:374903/1‑251 (MQ=255)
                                                                                                                                                   ggTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATc                                                                                              <  1:117420/251‑1 (MQ=255)
                                                                                                                                                    gtgtTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCa                                                                                             <  1:147129/251‑1 (MQ=255)
                                                                                                                                                         aaTTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGGCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCGCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCtt                                                                                        <  2:939/251‑1 (MQ=255)
                                                                                                                                                                      cTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GGTAGGC‑ATGACCCTTATTTCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGGATTGTTCAACGTTTCCACGCTACTGGTGTGGCGTCTGATTGGTTTGCTGTGGATTGGCGCAGGTAAAGTTAGTCTTCATTATCATGCTTATCTGACAgctgc                                                                           >  1:152587/1‑251 (MQ=255)
                                                                                                                                                                        gTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAgctgc                                                                           <  1:247521/249‑1 (MQ=255)
                                                                                                                                                                          ttAAGTTGGGCGAGCACTGGCAGCTTTATCATCGTCTTGAATCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑TCTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTc                                                                        <  1:90103/250‑1 (MQ=255)
                                                                                                                                                                           tAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCa                                                                       <  1:366496/250‑1 (MQ=255)
                                                                                                                                                                                         aCTGGCAGCTTGATCGTCGTCTTGAAGCGTGCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCAt                                                        <  2:4944/251‑1 (MQ=255)
                                                                                                                                                                                            ggCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACCTTTTTTTTTCCGGCGCAC‑GCTAGGCCAGGACCAT‑CTTGCGATGGGGGTG‑GATTTTTTAATGGTACCCGCATTGCAGGGGGTGTTGTTCACCGTTCCCCCACTCCTGGTGTGGGCGCTGCTTGGTTTGCTGTGTTTTGGCGCCGGTAAAGGTCGGCTTCCTTTTCCTTCTTTTCTGGTGGCTGCTCCACGTAATGATAACCCTGcc                                                     >  2:160865/1‑251 (MQ=255)
                                                                                                                                                                                             gCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGccc                                                    <  1:426042/251‑1 (MQ=255)
                                                                                                                                                                                               aGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTTTCTCGGGGGCAC‑GTTTGGC‑TTGACCATTATTGAGATGTGCGTG‑GCTTTTTTACTGTTAAGCGCATTTCTGGGGGTATTGTTTACCGTTCCCCAGCTACTGGTGTGGCTGCTGTTTGGTTTGCTGTGTATTTGGGCCGTGAGAGTTTGTCTTCATTTTTATTCTTTTCTTACGGCCGCTTCACGTACAGATATCCATGCCCgc                                                  >  2:357903/1‑250 (MQ=255)
                                                                                                                                                                                               aGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGt                                                  <  1:160865/251‑1 (MQ=255)
                                                                                                                                                                                                   tGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGAAGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGcc                                               <  2:356249/250‑1 (MQ=255)
                                                                                                                                                                                                             cTTGAAGCGTTCTTTCGGCGGGAG‑AAACATTTTTTTCTCGGGGGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTAAGCGCCTTTCAGGGAGTATTGTTCAACGTTCCCACGCTCCTGGTGTGGCTGCTGATTGGTTTTCTGTGTATTGGCGCAGGGAAAGTTCGGCTTCATTATCATGCTTATCTGAAAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACCAtggctgg                                    >  1:90751/1‑251 (MQ=255)
                                                                                                                                                                                                                  aGCGTTCTTTCGGCGGGTG‑AAACATTTTTTTCTCCGGCGCAC‑GGTAGGC‑ATCACCATTCTTTCGATGGGCGTG‑ACTTTTTTTCTGTTACGCGCCTTGCAGGGAGGATTGTTAAACGTTCCCACGCTACTGGGGTGGGTGCTGATTGGTTTTGTGTGTATTTGCGCAGGTCAAGGTCGTCTTCATTATCATGCTTTTTTGACAGCTGCTTCACGTAATGANAACCATGCCCGTGCCACGAtggcggggggac                               >  2:233840/1‑243 (MQ=255)
                                                                                                                                                                                                                     gTTCTTTCGGCGGGGG‑ACACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGGGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTc                             >  1:349522/1‑250 (MQ=255)
                                                                                                                                                                                                                          ttCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCAt                        <  2:135189/250‑1 (MQ=255)
                                                                                                                                                                                                                                        aCATTTTTCTCTCTGGCTCAC‑GTTAGTC‑CTGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGGTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTcc          <  2:72217/251‑1 (MQ=255)
                                                                                                                                                                                                                                               tctctcGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCGGCAg    <  1:397964/250‑1 (MQ=255)
                                                                                                                                                                                                                                                ctctcGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCGGCAGGc  >  1:365836/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                             
TGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGG‑TTTCGTGTGCTGGTCAGCAAGGAGACAACATG‑ACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCAC‑GTTAGGC‑ATGACCATTATTGCGATGGGCGTG‑ACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCGGCAGGC  >  NC_000913/119159‑119646

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