Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913119,3960CT28.6% 61.8 / 10.6 43S39F (TCT→TTT) ampEprotein AmpE
Reads supporting (aligned to +/- strand):  ref base C (8/22);  new base T (12/0);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTT‑CACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCG  >  NC_000913/119170‑119640
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          
tatCGCTAAAAACATGACGGGCCATTGTGCTATTGCACCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTACTTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑GGGAAACATTTTT‑TTTCTGGCGCAC‑G‑tatggggtgtc                                                                                                                                                                                                                                  >  2:521484/1‑241 (MQ=255)
                                    gcCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTg                                                                                                                                                                                                                          <  2:414029/223‑1 (MQ=255)
                                    gcCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑GG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTg                                                                                                                                                                                                                          >  1:414029/1‑223 (MQ=255)
                                      cGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑GG‑AAACATTTTTCTCCCGGGCGCACCG‑TT‑GGCAAGGCCCTTATTGCGGTGGGCGGGG‑Ctttttttttgttaagcg                                                                                                                                                                                           >  1:270364/1‑240 (MQ=255)
                                        gATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGTG‑AA‑CATTTTTTTTTCTGGCGCAC‑G‑TTATGCATGACCATTATTGCGGTGGGGGTGA‑Ccttttttctttttttccgc                                                                                                                                                                                         >  2:78383/1‑238 (MQ=255)
                                            ggAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑GG‑GAAAATTTTTTTTTCCGGGGCAC‑C‑TTTGGCCAGTCCCTTATTTCGGTGGGCGTGA‑Cttttttttcgttacccgcatttc                                                                                                                                                                                     >  2:240623/1‑234 (MQ=255)
                                             gAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTTCAACCTTACTGGAGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGACGCGTTCTTTCGGCGGC‑AG‑AAACCTTTTTTTCTCCGGCGCGC‑GATGAGGAAGGGCC‑TTTTTGCGATAGGCGTGG‑CTTTTTTACTGTTAggcggcctgc                                                                                                                                                                                     >  1:325744/1‑244 (MQ=255)
                                               aaaCCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGATACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCTCTGTCAGCGTGAGCAGCGGCTTGAAGGGGTCTATCGGCGTT‑TG‑AAACATTGTGCTCGCTGGCGCAC‑G‑TTAGCCAGGGCCATTACTGGTATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAgg                                                                                                                                                                                  <  1:245633/251‑1 (MQ=255)
                                                 aCCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑GG‑GAACATTTTTCTTTCGGGCGCAC‑G‑TCAGGGATGACCCTTATTTGGATGGGCGTGA‑TTTTTTTCCTGTTACCCGCCTTGCCgggg                                                                                                                                                                                >  1:405498/1‑250 (MQ=255)
                                                   cGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTACTTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑CG‑AAAAATTTTTCTCTTCGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGGGTGA‑CTTTTTTTCTGTTACGCGCATTGCCGGGAgg                                                                                                                                                                              >  1:193080/1‑250 (MQ=255)
                                                    gATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTa                                                                                                                                                                             <  1:201710/251‑1 (MQ=255)
                                                          ggTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTc                                                                                                                                                                       <  1:153352/251‑1 (MQ=255)
                                                               ggCATTTGATTGGGCTCGGTTTCGTGTTCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTGAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCTGCGGG‑TG‑AATCATTTTTCTCTCGGGCGCAG‑G‑TTAGGCATGACCATTTGTGCGAGGGGCGTGA‑CATTTTTACTGTTACGCGCATGGCAGGGAGTATTGTGCAACGt                                                                                                                                                                  <  2:467977/250‑1 (MQ=255)
                                                                gCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGt                                                                                                                                                                  <  1:52495/250‑1 (MQ=255)
                                                                  aTTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑GG‑AAACATTTTTTTCTCGGGCGCAC‑G‑TTAGGCCTGACCATTATTTCGATGGGGGTGG‑CTTTTTTTCTGTTTCGCGCCTTGCAGGGAGTGTTGTTTAActtttc                                                                                                                                                               >  2:134639/1‑246 (MQ=255)
                                                                  aTTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑GG‑AAACATTTTTTTCTCGGGCGCAC‑G‑GTAGGCATGACCATTATTGCGATGGGGGGGG‑CTTTTTTTCTGTTACGCGCATTGCAGGGGGTATTGTTTAACGtttc                                                                                                                                                               >  1:493616/1‑249 (MQ=255)
                                                                  aTTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGAAAACATGACGCTATTTACAACCTTACGGGTGTTAATGTTCGAGCGCCTGGTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGCATGTTTCAACGTTcc                                                                                                                                                               <  2:433118/251‑1 (MQ=255)
                                                                      gATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑Ac                                                                                                                                                           <  1:373196/250‑1 (MQ=255)
                                                                            gCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑GG‑AAAAATTTTTCTCTCGGGCGCAC‑G‑TTTGGCATGACCCTTATTGCGATGGGGGTGA‑CTTTTTTTCTGTTACGCGCATTTCAGGGGGTATTTTTCAACGTTCCC‑ACGCTACTg                                                                                                                                                    >  1:236784/1‑251 (MQ=255)
                                                                                     cGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGgct                                                                                                                                            <  1:113040/250‑1 (MQ=255)
                                                                                      gtgtGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑GG‑GAACATTTTTTTTTCTGGCGCAC‑G‑TTAGGCATGACCCTTATTTCGATGGGGGTGG‑CTTTTTTTTTGTTTTGCGCCTTGCCGGGAGGATTGTTCAAAGTTTCC‑CCCCTACTTGTGTGgctgg                                                                                                                                          >  2:5539/1‑250 (MQ=255)
                                                                                        ttgCTGGTCAGCAAGGAGCCAACATGACGCTATTTACAGCCTTACTGGTGTTAAGGTTCGAGCGCCGGTTTAAGGTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCGCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGgctgct                                                                                                                                         <  2:399467/249‑1 (MQ=255)
                                                                                                        agaCAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTGCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCAt                                                                                                                                                                                                                                     <  1:91748/144‑1 (MQ=255)
                                                                                                           caacaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTg                                                                                                                    <  1:521484/251‑1 (MQ=255)
                                                                                                           caacaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑GG‑AAACATTTTTTTCTCGGGCGCAC‑C‑TTAGGCATGACCCTTATTGCGATGGGCGGGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTAttt                                                                                                                    >  1:403335/1‑250 (MQ=255)
                                                                                                             acaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCCC‑G‑TTAGGCATGACCATTATTTCGATGGGCGTGA‑CTTTTTTTCTGTTACGCGCCTTGCAGGGGGGATTGTTCAACGTTCCC‑ACGCTACTGGTGGGGCTGCTG‑ATTGGTTTGCTGTGTATTggg                                                                                                                  >  1:162675/1‑250 (MQ=255)
                                                                                                             acaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTGgc                                                                                                                  <  1:404985/251‑1 (MQ=255)
                                                                                                               aTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTTTCCCTGGCGCCC‑G‑TTAGGCATGACCCTTTTTTCGATGGGCGTGA‑CTTTTTTTCTGTTACGCGCATTTCAGGGGGGATTGTTTAACGTTTCC‑ACGCTACTGTGGTGGCTGCTG‑GTTTGTTTTCtgggttatgtggc                                                                                                                >  2:167009/1‑240 (MQ=255)
                                                                                                                        tttACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGACGCGTTCTTTCGGCGGG‑GG‑AAACATTTTTTTTTCTGGCGCAC‑G‑TTTAGCATGGACACTATTGGGGTGGGCGTGAACTTTTTT‑CTTTTATGCGCATTGCTGGGGGGGTTGTTCCACGTTCCCCAC‑CCACTTGTGTGGGTGCTG‑ATTGTTTTGCTGTGTATTgggggcgggacagg                                                                                                       >  2:480987/1‑239 (MQ=255)
                                                                                                                          tACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCACCGTCTTGAAGCGTTCTTTCGGCGGG‑TT‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACTATTATTGCGATGGGCGTGT‑GTGTGTGACTGTTACGCGCATTTGAGTTAGTATTGTTCAACGTTCCC‑ACGCTTCTGGTGGGGCTGCTG‑ATTGGTGTGCTGTGTATTGGCGCAGGTAAAGtt                                                                                                      <  1:240244/250‑1 (MQ=255)
                                                                                                                           aCAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCg                                                                                                    <  1:225328/251‑1 (MQ=255)
                                                                                                                              aCCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCGCTCGGGCGCAC‑G‑TTAGGCATGACCATTAGTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTGCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCt                                                                                                 <  1:518921/251‑1 (MQ=255)
                                                                                                                                   aCTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCAt                                                                                             <  1:5539/250‑1 (MQ=255)
                                                                                                                                        ggtTAATTTTCGAGCGCCTGTTGAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGGTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCAt                                                                                       <  1:80894/250‑1 (MQ=255)
                                                                                                                                            aaTTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCt                                                                                    <  1:327977/250‑1 (MQ=255)
                                                                                                                                            aaTTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑CG‑AAACATTTTTTTCTTGGGCGGCC‑G‑TTAGGGCTGAACATTATTGCGATGGGCGGGA‑CTTTTTTACTGTTACGCGCATTGCAGGGGGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑GTTGGTTTGCTGTGTATTGGCGCAGGGAAAGTTCGTCTTCCTTATCATGCtt                                                                                   >  1:498395/1‑251 (MQ=255)
                                                                                                                                                           gTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATGGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCtt                                                                    <  2:493616/251‑1 (MQ=255)
                                                                                                                                                                        gAGCACTAGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑CG‑AAACATTTTTTTCTCGGGCGCAC‑G‑ATAGGCCTGACCATTTTTGCGATGGGGGTGA‑CTTTTTTTCTTTTACGCGCCTTGCCGGGGGTATTGTTCAACGTTCCC‑CCGCTACTGGTGTGGCTGCTG‑GTTGGTTTGCTGTGTATTGTCGCAGGGAAAGTTCGTCTTCCTTTTCATTCTTTTCTGACAGCTGCTT‑TACGTAAtaaaa                                                       >  2:239988/1‑246 (MQ=255)
                                                                                                                                                                              tGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAAAATTTTTTTCTCCGGCGCAC‑G‑GTTGGCATGTCCATTTTTGCGGTGGGCGTGA‑CTTTTTTTCTTTTACGCGCATTGGCGGGAGGATTGTTTAAAGTTCCC‑ACGCTACTGGTGTGGGTGCTGGTTTGGTTTG‑TGTGTTTTTGCGCCGGGAAAGTTCGTTGTCCATATTATGCTTATCTGCCAGACGCTTTCA‑GTAATGATAGcccttg                                                >  2:396367/1‑247 (MQ=255)
                                                                                                                                                                                                                        g‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CGTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTGGCGCAGGTAAAGGTCGTCTTCATTATCATGCTTATCTGACAGCTTCTT‑CACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGAGTCACGGc        <  2:193080/250‑1 (MQ=255)
                                                                                                                                                                                                                          aaaCATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTCTTCCTTTCACGCGCATTGCAGGGACTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTT‑CACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGACTTCACCATGATTCACGGCGt      >  2:410542/1‑251 (MQ=255)
                                                                                                                                                                                                                            aCATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTT‑CACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTcc    <  1:134639/251‑1 (MQ=255)
                                                                                                                                                                                                                              aTTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTT‑CACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCg  >  1:379054/1‑251 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          
TATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGG‑TG‑AAACATTTTTCTCTCGGGCGCAC‑G‑TTAGGCATGACCATTATTGCGATGGGCGTGA‑CTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCC‑ACGCTACTGGTGTGGCTGCTG‑ATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTT‑CACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCG  >  NC_000913/119170‑119640

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