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.

GATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCC  >  NC_000913/119156‑119639
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       
gattgCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGgcg                                                                                                                                                                                                                                                <  1:562599/250‑1 (MQ=255)
    gCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACg                                                                                                                                                                                                                                            <  1:556083/250‑1 (MQ=255)
            gATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGTTCCTGCATTTGATTGGGCCCTGTTTCGTGTGCTGGTCAGCAAGGAGACAACACGACGCTATTTACAACCTTACTGGTGTTCATTTTCGATCCCCTTTTTAAGTTGGGCGTGCACTGGCAGCTTCAGCATCGTCTCCAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATg                                                                                                                                                                                                                                  <  1:252369/251‑1 (MQ=255)
                   cTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGTAAAATTTTTTTCTCCGgcgccagttgggcaagtgcactt                                                                                                                                                                                                                            >  2:202033/1‑232 (MQ=255)
                    tAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTCTTTCGGGCGCCCGTTAGGCATGCCCattat                                                                                                                                                                                                                          >  1:99860/1‑251 (MQ=255)
                       aaaCATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGGAACATTTTTTTTTCGGgcgccccgtaggcattgccatttttgg                                                                                                                                                                                                                        >  1:531619/1‑228 (MQ=255)
                       aaaCATGACGGGCCATTGTGATAGTGCGCCGGATCGGAAAACCGATCCCGGGCCTGCATTTGATTGGGCACGGTTTCGTGTGCGGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTGAGGCATGACCATTATTGc                                                                                                                                                                                                                       <  2:96759/251‑1 (MQ=255)
                               cGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCg                                                                                                                                                                                                               <  1:441088/251‑1 (MQ=255)
                                        gtgATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACttttt                                                                                                                                                                                                     <  1:412727/251‑1 (MQ=255)
                                               tGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGt                                                                                                                                                                                              <  1:385151/250‑1 (MQ=255)
                                                  gcCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGGAAAAATTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTTCGGTGGGCGGGGGTTTTTTTCTGtttcg                                                                                                                                                                                          >  1:137183/1‑248 (MQ=255)
                                                     ggATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACgcgc                                                                                                                                                                                       <  1:18855/251‑1 (MQ=255)
                                                       aTCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCAt                                                                                                                                                                                     <  1:149404/251‑1 (MQ=255)
                                                         cGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCAtt                                                                                                                                                                                    <  1:372840/250‑1 (MQ=255)
                                                                  gATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTTTCTCCGGCGCCCGGTAGGCCTGACCATTATTGCGATGGGCGGGCCTTTTTTTCTGTTACGCGCATTGCAGGGAgg                                                                                                                                                                          >  1:66611/1‑249 (MQ=255)
                                                                    tCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCCTGACCCTTATTGCGGTGGGCGGGACTTTTTTTCTGTTACGCGCATTGCCGGGGGTAtt                                                                                                                                                                       >  1:512144/1‑251 (MQ=255)
                                                                      ccGGTCCTGCGTTTGATTGGGCACTGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTGTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGTGATGGGCGTGACTTTTTTACTGTTACGCGCATGGCAGGGAGTACTGt                                                                                                                                                                     <  2:11841/251‑1 (MQ=255)
                                                                       cGGTCCTGCATTTGATGGGGCACGGTTTCGTGGGCGGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATGTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGtt                                                                                                                                                                    <  2:66611/251‑1 (MQ=255)
                                                                            cTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTTTCTCGGGCGCCCGGTAGGCATGACCATTATTTCGGTGGGCGTTGCTTTTTTACTGTTTCGCGCAATTCAGGGGGTATTGtttca                                                                                                                                                                 >  1:552040/1‑246 (MQ=255)
                                                                              gCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGtt                                                                                                                                                             <  1:457904/251‑1 (MQ=255)
                                                                                   tGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTCTATCTAGCGCACGTTAGGGATGACCATTATTGCGCTGGGGGGGAACTTTTTTATGTTAAGCGCATTGCAGGGCGTTTTGTTCCACGTTACCAc                                                                                                                                                        >  1:370225/1‑251 (MQ=255)
                                                                                    gATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCAc                                                                                                                                                        <  1:281140/250‑1 (MQ=255)
                                                                                     aTTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGTCAACATGACGTTATTTACAACCTTACTGGTGTTAATTTTCGAGCGGTTGGTTAAGGTGTGCGAGCACTGGCAGCTTTATCATCGTCTTGAAGTGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCTCGTTAGGCATGACCATTATTGCGATGGGCGTGTCTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGc                                                                                                                                                      <  1:39553/251‑1 (MQ=255)
                                                                                           cACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACAATTTTTTTTCGGGCGCCCGTTAGGGATTGACCTTTTTGCGGTGGGCGTGGCTTTTTTTCTGTTTCGCGCCTTTCAGGGGGGATTGtttcacggtcccaacccccctg                                                                                                                                                >  2:380544/1‑231 (MQ=255)
                                                                                              ggTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGgtgt                                                                                                                                             <  1:196099/251‑1 (MQ=255)
                                                                                                            gTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGGAACATTTTTCTTTCGGGCGCACGTTAGGGCTGACCATTTTTTCGGTGGGGGGGACTTTTTTACTGTTACGCGCATTGCCAGGGGTATTGTTTAAACTTCCCCCGCTACTgggggggcggctggtttgt                                                                                                                               >  1:46709/1‑234 (MQ=255)
                                                                                                               aGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTg                                                                                                                            <  1:46839/251‑1 (MQ=255)
                                                                                                                  aaGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACt                                                                                                                                                                                                         <  1:81852/173‑1 (MQ=255)
                                                                                                                  aaGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCCTGACCCTTATTTCGGTGGGCGTGACt                                                                                                                                                                                                         >  2:81852/1‑173 (MQ=255)
                                                                                                                   aGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCCGGCGCACGTTTGGCATGACCATTATTTCGATGGGCGTGCCTTTTTTTCTTGTTCGCGCATTACCGGGGGTATTGTTCAACGTTTCCACGCTTCCGGTGTGGGTGGTGATTGGTTTGCtgg                                                                                                                        >  1:116419/1‑250 (MQ=255)
                                                                                                                           acaTGACGCTAGTTACAACCTTACTGGTGTTAGTTTTCGAGCGCCGGTTGAAGTGGGGCGAGCACTGGCAGCGTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGgc                                                                                                                <  2:60516/251‑1 (MQ=255)
                                                                                                                            aaTGACGCTATTTACAAACTTACTGGTGTTAATTTTCGAGCGCCTGGTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGGCTTGAGGCGGTCTTTCGGGGGGTGAAACATTTTTCTCGCGGGCGCTCGTTAGGCATGTCCATTATTGCTATGGGCGTGACTTTTTTACTGTTACGCGCTTTGCATGTAGTGTTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGgcg                                                                                                               <  2:217348/250‑1 (MQ=255)
                                                                                                                                     aTTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTaa                                                                                                        <  1:461281/249‑1 (MQ=255)
                                                                                                                                     aTTTACAACCTTACGGGTGTTAATTTTCGAGCTCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCTCACGTTAGGCATGACCATTATTGCGATGTGCGGGACTTTTTTACTGTGACGCGCAGTGCAGGGAGTATTGTTCAACGTTCCCACGCTGCTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTTGCGCAGGTaa                                                                                                        <  2:47035/249‑1 (MQ=255)
                                                                                                                                       ttACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTACGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAAAAATTTTTTCTCTGGCGCACGTTAGGCATGACCATTAATGCGCTGGGGGGGACTTTTTTTCTGTTTCGCGCATTGCAGGGGGGAGTTTTCAACGGTTCCACGCTACTGGGGGGGCGTCTGCTTGGGTTTCTTTGTATTGTCGCAGGGAAAggt                                                                                                    >  2:142635/1‑249 (MQ=255)
                                                                                                                                           aaCCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCATGACCCTTATTGCGGTGGGCGGGACTTTTTTACTGTTACGCGCATTGCAGGGGGGATTGTTCAACGTTCCCCCGCCACTGGTGGGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGGAAAGTTCGTc                                                                                                >  1:462096/1‑251 (MQ=255)
                                                                                                                                            aCCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTTCGATGGGCGTGACTTTTTTTCTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCt                                                                                               >  1:6498/1‑251 (MQ=255)
                                                                                                                                                       gtTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCCCGTTAGGCATGACCATTATTTCGGTTGGGGTGACTTTTTTATTGTTACTCGCATTTTAGGGGGTATTTTTCAACGTTCCCACGCTACTTGTGTGGCTGCTGTTTGGTTTGTTGTGTATTGGCGCCGGGAAAGTTCGTCTTCAtttttatg                                                                                    >  1:390216/1‑245 (MQ=255)
                                                                                                                                                          aaTTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCtt                                                                                 <  1:515619/251‑1 (MQ=255)
                                                                                                                                                           aTTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACCTTTTTCTCTCCGGCGCACGTTAAGCATGACCCTTATTGCGGTGGGCGGGACTTTTTTTTTTTTACCCGCCTTGCAGGGGGGATTTTTTCACGTTTCCACCCTACCGGTGTGGCTGCTGATTGGTTTTTGGTGTATTGGGTCCGGGAAAGGTCGTCTTCATTTTCATggctt                                                                                >  2:120249/1‑247 (MQ=255)
                                                                                                                                                           aTTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACAATTTTTTTTCGGGGGCACGCTAGGCATGACCATTATTTCGATGGGGGGGACTTTTTTACTGTTTCGCGCATTGGAGGGGGTATTGTTCAACGTTCCCACGCTACTGGTGGGGCTTCTGGTTGGGTTGCTGTGTATTAGCGCAGGGAAAGTTCGTCTTCATTTTCTTGCttt                                                                                >  1:566834/1‑250 (MQ=255)
                                                                                                                                                                cGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTg                                                                            <  1:548553/250‑1 (MQ=255)
                                                                                                                                                                      ccTGTTTAAGTTGGGCGGGCCCTGGCAGCTTGATCATCGTCTTGAAGCGTTCGTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAgctg                                                                     <  2:33578/251‑1 (MQ=255)
                                                                                                                                                                           ttAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTACGGCGGGGGAAACATTTTTTTCTCGGGCGCACGCTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCCGGGAGGATTGTTCAACGGTCCCACGCTACTGGGGTGGCTGCTGATTGGTTTGCTGTGTTTTTGCGCAGGGAAAGGTCGTCTTCATTATTCTGCTTATCTTACAGCCGCtttc                                                                >  1:316401/1‑249 (MQ=255)
                                                                                                                                                                               gTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTTCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTa                                                            <  1:515610/251‑1 (MQ=255)
                                                                                                                                                                                                   ttGATCATCGTCTTGAAGCGTTCTTTCGGCGGGAGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGAGATGGGCGTGACTTTTTTTCTGTTTCGCGCATTGCAGGGGGGATTGTTCAACGTTCCCCCGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTTTTTGGGCAGGGAAAGTTCGTTTTCATTTTCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGcc                                        >  1:443756/1‑251 (MQ=255)
                                                                                                                                                                                                     gATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCa                                       <  2:390216/250‑1 (MQ=255)
                                                                                                                                                                                                     gATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCAc                                      <  1:127193/251‑1 (MQ=255)
                                                                                                                                                                                                     gATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCAc                                      <  2:531619/251‑1 (MQ=255)
                                                                                                                                                                                                        cgtcGTCTTGAAGCGTTCTTGCGGCGGGGGANACATTTTTCTCTCTGGCGCACGGTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGTCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGAt                                   <  2:260388/249‑1 (MQ=255)
                                                                                                                                                                                                                 gAAGCGTTCTTGCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCCTCATGATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGa                          <  1:184682/251‑1 (MQ=255)
                                                                                                                                                                                                                        tCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGCGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTTTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGGCTTCATTATCATGCTTATCTTGCAGCTGCTTCACGTAATGATAGCCCTGCCCGTGCCACGATGGCTGGCGAACTCAcc                   >  1:265619/1‑251 (MQ=255)
                                                                                                                                                                                                                             cggcggGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGAt              <  1:121737/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggcggGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGAtt             >  1:41740/1‑251 (MQ=255)
                                                                                                                                                                                                                                cggGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCa           <  2:62697/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gggagAAACATTTTTCTCTCGGGCGCACGATATGCATGACCATTATTGCGATGGGCGTGACTTTTTTCCTCTTACGCGCATCGCAGGGAGTATTGTTCCACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCAGGCTTATCTGACAGCTGCTTCACGTAATGATATCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCAc          >  2:543766/5‑251 (MQ=255)
                                                                                                                                                                                                                                   gTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACgg        >  1:79848/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aaaCATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGt     <  2:137183/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cctttttCTCTCGGGCGCACGGTAGGCATGACCATTATTACGCTGGGCGTGACTTTTTTACTGTTACGCGCATTGAAGGGAGTATTGTTCAACGTTCACACGCTACTGGAGTGGCTGCTGATTGGATTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTCTCTGACAGCTGCTTCACGTAATGATAGCAATGCCCGTGCCACAGTGGCTGGCGACCTACCCATGATTCACGGCGGccc  >  2:79245/3‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       
GATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCC  >  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: