Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913119,3960CT20.3% 121.8 / 11.3 59S39F (TCT→TTT) ampEprotein AmpE
Reads supporting (aligned to +/- strand):  ref base C (13/34);  new base T (12/0);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 4.65e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCC  >  NC_000913/119156‑119639
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       
gattgCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGgcg                                                                                                                                                                                                                                                <  1:562599/250‑1 (MQ=255)
    gCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACg                                                                                                                                                                                                                                            <  1:556083/250‑1 (MQ=255)
            gATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGTTCCTGCATTTGATTGGGCCCTGTTTCGTGTGCTGGTCAGCAAGGAGACAACACGACGCTATTTACAACCTTACTGGTGTTCATTTTCGATCCCCTTTTTAAGTTGGGCGTGCACTGGCAGCTTCAGCATCGTCTCCAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑g                                                                                                                                                                                                                                  <  1:252369/251‑1 (MQ=255)
                   cTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGTAAAATTTTTTTCTCCGgcgccagttgggcaagtgcactt                                                                                                                                                                                                                            >  2:202033/1‑232 (MQ=255)
                    tAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTCTTTCGGGCGCCCGTTAGGCAT‑GCCCattat                                                                                                                                                                                                                          >  1:99860/1‑251 (MQ=255)
                       aaaCATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGGAACATTTTTTTTTCGGgcgccccgtaggcattgccatttttgg                                                                                                                                                                                                                        >  1:531619/1‑228 (MQ=255)
                       aaaCATGACGGGCCATTGTGATAGTGCGCCGGATCGGAAAACCGATCCCGGGCCTGCATTTGATTGGGCACGGTTTCGTGTGCGGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTGAGGCAT‑GACCATTATTGc                                                                                                                                                                                                                       <  2:96759/251‑1 (MQ=255)
                               cGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCg                                                                                                                                                                                                               <  1:441088/251‑1 (MQ=255)
                                        gtgATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑Cttttt                                                                                                                                                                                                     <  1:412727/251‑1 (MQ=255)
                                               tGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑Gt                                                                                                                                                                                              <  1:385151/250‑1 (MQ=255)
                                                  gcCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGGAAAAATTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTTCGGTGGGCGGGG‑GTTTTTTTCT‑Gtttcg                                                                                                                                                                                          >  1:137183/1‑248 (MQ=255)
                                                     ggATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACgcgc                                                                                                                                                                                       <  1:18855/251‑1 (MQ=255)
                                                       aTCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCAt                                                                                                                                                                                     <  1:149404/251‑1 (MQ=255)
                                                         cGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCAtt                                                                                                                                                                                    <  1:372840/250‑1 (MQ=255)
                                                                  gATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTTTCTCCGGCGCCCGGTAGGCCT‑GACCATTATTGCGATGGGCGGGC‑CTTTTTTTCT‑GTTACGCGCATTGC‑AGGGAgg                                                                                                                                                                          >  1:66611/1‑249 (MQ=255)
                                                                    tCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCCT‑GACCCTTATTGCGGTGGGCGGGA‑CTTTTTTTCT‑GTTACGCGCATTGC‑CGGGGGTAtt                                                                                                                                                                       >  1:512144/1‑251 (MQ=255)
                                                                      ccGGTCCTGCGTTTGATTGGGCACTGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTGTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGTGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATGGC‑AGGGAGTACTGt                                                                                                                                                                     <  2:11841/251‑1 (MQ=255)
                                                                       cGGTCCTGCATTTGATGGGGCACGGTTTCGTGGGCGGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATGTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGtt                                                                                                                                                                    <  2:66611/251‑1 (MQ=255)
                                                                            cTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTTTCTCGGGCGCCCGGTAGGCAT‑GACCATTATTTCGGTGGGCGTTG‑CTTTTTTACT‑GTTTCGCGCAATTC‑AGGGGGTATTGtttca                                                                                                                                                                 >  1:552040/1‑246 (MQ=255)
                                                                              gCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGtt                                                                                                                                                             <  1:457904/251‑1 (MQ=255)
                                                                                   tGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTCTATCTAGCGCACGTTAGGGAT‑GACCATTATTGCGCTGGGGGGGAACTTTTTTA‑T‑GTTAAGCGCATTGC‑AGGGCGTTTTGTTCCACGTTACCAc                                                                                                                                                        >  1:370225/1‑251 (MQ=255)
                                                                                    gATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCAc                                                                                                                                                        <  1:281140/250‑1 (MQ=255)
                                                                                     aTTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGTCAACATGACGTTATTTACAACCTTACTGGTGTTAATTTTCGAGCGGTTGGTTAAGGTGTGCGAGCACTGGCAGCTTTATCATCGTCTTGAAGTGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCTCGTTAGGCAT‑GACCATTATTGCGATGGGCGTGT‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGc                                                                                                                                                      <  1:39553/251‑1 (MQ=255)
                                                                                           cACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACAATTTTTTTTCGGGCGCCCGTTAGGGATTGACC‑TTTTTGCGGTGGGCGTGG‑CTTTTTTTCT‑GTTTCGCGCCTTTC‑AGGGGGGATTGtttcacggtcccaacccccctg                                                                                                                                                >  2:380544/1‑231 (MQ=255)
                                                                                              ggTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGgtgt                                                                                                                                             <  1:196099/251‑1 (MQ=255)
                                                                                                            gTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGGAACATTTTTCTTTCGGGCGCACGTTAGGGCT‑GACCATTTTTTCGGTGGGGGGGA‑CTTTTTTACT‑GTTACGCGCATTGCCAGGG‑GTATTGTTTAAACTTCCCCCGCTACTgggggggcggctggtttgt                                                                                                                               >  1:46709/1‑234 (MQ=255)
                                                                                                               aGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTg                                                                                                                            <  1:46839/251‑1 (MQ=255)
                                                                                                                  aaGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑Ct                                                                                                                                                                                                         <  1:81852/173‑1 (MQ=255)
                                                                                                                  aaGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCCT‑GACCCTTATTTCGGTGGGCGTGA‑Ct                                                                                                                                                                                                         >  2:81852/1‑173 (MQ=255)
                                                                                                                   aGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCCGGCGCACGTTTGGCAT‑GACCATTATTTCGATGGGCGTGC‑CTTTTTTTCTTGTT‑CGCGCATTAC‑CGGGGGTATTGTTCAACGTTTCCACGCTTCCGGTGTGGGTGGTGATTGGTTTGCtgg                                                                                                                        >  1:116419/1‑250 (MQ=255)
                                                                                                                           acaTGACGCTAGTTACAACCTTACTGGTGTTAGTTTTCGAGCGCCGGTTGAAGTGGGGCGAGCACTGGCAGCGTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGgc                                                                                                                <  2:60516/251‑1 (MQ=255)
                                                                                                                            aaTGACGCTATTTACAAACTTACTGGTGTTAATTTTCGAGCGCCTGGTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGGCTTGAGGCGGTCTTTCGGGGGGTGAAACATTTTTCTCGCGGGCGCTCGTTAGGCAT‑GTCCATTATTGCTATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCTTTGC‑ATGTAGTGTTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGgcg                                                                                                               <  2:217348/250‑1 (MQ=255)
                                                                                                                                     aTTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTaa                                                                                                        <  1:461281/249‑1 (MQ=255)
                                                                                                                                     aTTTACAACCTTACGGGTGTTAATTTTCGAGCTCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCTCACGTTAGGCAT‑GACCATTATTGCGATGTGCGGGA‑CTTTTTTACT‑GTGACGCGCAGTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTGCTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTTGCGCAGGTaa                                                                                                        <  2:47035/249‑1 (MQ=255)
                                                                                                                                       ttACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTACGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAAAAATTTTTTCTCTGGCGCACGTTAGGCAT‑GACCATTAATGCGCTGGGGGGGA‑CTTTTTTTCT‑GTTTCGCGCATTGC‑AGGGGGGAGTTTTCAACGGTTCCACGCTACTGGGGGGGCGTCTGCTTGGGTTTCTTTGTATTGTCGCAGGGAAAggt                                                                                                    >  2:142635/1‑249 (MQ=255)
                                                                                                                                           aaCCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCAT‑GACCCTTATTGCGGTGGGCGGGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGGGGATTGTTCAACGTTCCCCCGCCACTGGTGGGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGGAAAGTTCGTc                                                                                                >  1:462096/1‑251 (MQ=255)
                                                                                                                                            aCCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTTCGATGGGCGTGA‑CTTTTTTTCT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCt                                                                                               >  1:6498/1‑251 (MQ=255)
                                                                                                                                                       gtTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCCCGTTAGGCAT‑GACCATTATTTCGGTTGGGGTGA‑CTTTTTTATT‑GTTACTCGCATTTT‑AGGGGGTATTTTTCAACGTTCCCACGCTACTTGTGTGGCTGCTGTTTGGTTTGTTGTGTATTGGCGCCGGGAAAGTTCGTCTTCAtttttatg                                                                                    >  1:390216/1‑245 (MQ=255)
                                                                                                                                                          aaTTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCtt                                                                                 <  1:515619/251‑1 (MQ=255)
                                                                                                                                                           aTTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACCTTTTTCTCTCCGGCGCACGTTAAGCAT‑GACCCTTATTGCGGTGGGCGGGA‑CTTTTTTTTT‑TTTACCCGCCTTGC‑AGGGGGGATTTTTTCACGTTTCCACCCTACCGGTGTGGCTGCTGATTGGTTTTTGGTGTATTGGGTCCGGGAAAGGTCGTCTTCATTTTCATggctt                                                                                >  2:120249/1‑247 (MQ=255)
                                                                                                                                                           aTTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACAATTTTTTTTCGGGGGCACGCTAGGCAT‑GACCATTATTTCGATGGGGGGGA‑CTTTTTTACT‑GTTTCGCGCATTGG‑AGGGGGTATTGTTCAACGTTCCCACGCTACTGGTGGGGCTTCTGGTTGGGTTGCTGTGTATTAGCGCAGGGAAAGTTCGTCTTCATTTTCTTGCttt                                                                                >  1:566834/1‑250 (MQ=255)
                                                                                                                                                                cGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTg                                                                            <  1:548553/250‑1 (MQ=255)
                                                                                                                                                                      ccTGTTTAAGTTGGGCGGGCCCTGGCAGCTTGATCATCGTCTTGAAGCGTTCGTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAgctg                                                                     <  2:33578/251‑1 (MQ=255)
                                                                                                                                                                           ttAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTACGGCGGGGGAAACATTTTTTTCTCGGGCGCACGCTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑CGGGAGGATTGTTCAACGGTCCCACGCTACTGGGGTGGCTGCTGATTGGTTTGCTGTGTTTTTGCGCAGGGAAAGGTCGTCTTCATTATTCTGCTTATCTTACAGCCGCtttc                                                                >  1:316401/1‑249 (MQ=255)
                                                                                                                                                                               gTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTTCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTa                                                            <  1:515610/251‑1 (MQ=255)
                                                                                                                                                                                                   ttGATCATCGTCTTGAAGCGTTCTTTCGGCGGGAGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGAGATGGGCGTGA‑CTTTTTTTCT‑GTTTCGCGCATTGC‑AGGGGGGATTGTTCAACGTTCCCCCGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTTTTTGGGCAGGGAAAGTTCGTTTTCATTTTCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGcc                                        >  1:443756/1‑251 (MQ=255)
                                                                                                                                                                                                     gATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCa                                       <  2:390216/250‑1 (MQ=255)
                                                                                                                                                                                                     gATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCAc                                      <  1:127193/251‑1 (MQ=255)
                                                                                                                                                                                                     gATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCAc                                      <  2:531619/251‑1 (MQ=255)
                                                                                                                                                                                                        cgtcGTCTTGAAGCGTTCTTGCGGCGGGGGANACATTTTTCTCTCTGGCGCACGGTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGTCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGAt                                   <  2:260388/249‑1 (MQ=255)
                                                                                                                                                                                                                 gAAGCGTTCTTGCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCCTCATGATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGa                          <  1:184682/251‑1 (MQ=255)
                                                                                                                                                                                                                        tCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGCGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTTTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGGCTTCATTATCATGCTTATCTTGCAGCTGCTTCACGTAATGATAGCCCTGCCCGTGCCACGATGGCTGGCGAACTCAcc                   >  1:265619/1‑251 (MQ=255)
                                                                                                                                                                                                                             cggcggGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGAt              <  1:121737/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggcggGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGAtt             >  1:41740/1‑251 (MQ=255)
                                                                                                                                                                                                                                cggGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCa           <  2:62697/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gggagAAACATTTTTCTCTCGGGCGCACGATATGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTCCT‑CTTACGCGCATCGC‑AGGGAGTATTGTTCCACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCAGGCTTATCTGACAGCTGCTTCACGTAATGATATCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCAc          >  2:543766/5‑251 (MQ=255)
                                                                                                                                                                                                                                   gTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACgg        >  1:79848/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aaaCATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGt     <  2:137183/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cctttttCTCTCGGGCGCACGGTAGGCAT‑GACCATTATTACGCTGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGA‑AGGGAGTATTGTTCAACGTTCACACGCTACTGGAGTGGCTGCTGATTGGATTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTCTCTGACAGCTGCTTCACGTAATGATAGCAATGCCCGTGCCACAGTGGCTGGCGACCTACCCATGATTCACGGCGGccc  >  2:79245/3‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       
GATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAT‑GACCATTATTGCGATGGGCGTGA‑CTTTTTTACT‑GTTACGCGCATTGC‑AGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCC  >  NC_000913/119156‑119639

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: ‑