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

GATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCGGCAG  >  NC_000913/119152‑119644
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           
ggttgatggCTATCCGGATAGCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCgg                                                                                                                                                                                                                                                       <  2:180617/249‑1 (MQ=255)
gattgattGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCATGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCgg                                                                                                                                                                                                                                                       <  1:63673/251‑1 (MQ=255)
 attgattgCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAAAATTTTTTTTTCggg                                                                                                                                                                                                                                                      >  1:409492/1‑251 (MQ=255)
   tgattgCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAGACATTTTtttctcttgcg                                                                                                                                                                                                                                                    >  2:177981/1‑246 (MQ=255)
   tgattgCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCggggg                                                                                                                                                                                                                                                    >  1:163526/1‑249 (MQ=255)
                gATATCGCTAAAAACATGACGGGCCATTGTGATATTGCTACGGTGCGGAAAACCGCTCCCGGTCTGTCATTTGATTGGTCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGGTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGGTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCaa                                                                                                                                                                                                                                        <  2:242491/250‑2 (MQ=255)
                         aaaaaCATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCatta                                                                                                                                                                                                                               <  1:437695/249‑1 (MQ=255)
                            aaCATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGGGACATTTTTTTCTCGGGCGCGCGTTAGGCATAGACATTATTgtg                                                                                                                                                                                                                          >  1:430323/1‑249 (MQ=255)
                            aaCATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCg                                                                                                                                                                                                                          <  1:298173/251‑1 (MQ=255)
                                 gACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATGTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATgg                                                                                                                                                                                                                      <  1:210305/250‑1 (MQ=255)
                                    gggCCATTGTGATATGGCGCCGGATCGGAAAACCGATCCTGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGt                                                                                                                                                                                                                  <  2:154791/251‑1 (MQ=255)
                                     ggCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGTGGAAACCATTTTTTTTTCGGTGGACGTTAGGCGTGACCGTTATTTTGATggggggg                                                                                                                                                                                                                 >  2:375772/1‑247 (MQ=255)
                                     ggCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCATGACACTTATTggggtgggggtg                                                                                                                                                                                                                 >  1:35071/1‑240 (MQ=255)
                                                 atTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTg                                                                                                                                                                                                     <  1:550902/251‑1 (MQ=255)
                                                      gcCGGTTCGGAAAAACGATTCCGGTTCTGCATTTGATTGGGCTCGGTGTCGTGTGCTGGGTCGCAAGGAGACAACATGACGCTATTTACATCCTTACTGGTGTTAATTTGCGAGCGCCTGGTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGTGCGTGACTTTTTTACTGTTAcg                                                                                                                                                                                                <  2:415165/251‑1 (MQ=255)
                                                             cGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGAGGGGCGTGACTTTTTTACTGTTACGCGCAtt                                                                                                                                                                                          <  2:109445/250‑1 (MQ=255)
                                                                  aaCCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGGAACATTTTTTTTTCGGGCGCACGTTAGGCATTACCCTTTTTTCGGGTGGCGTGACTTTTTTTCTTTTTCgcgcctttccggg                                                                                                                                                                                    >  1:564323/1‑242 (MQ=255)
                                                                   aCCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGTCCATTATTGCGATGGGCGTTACTTTTTTACTGTTACGCGCATTGCAGGGa                                                                                                                                                                                  <  1:457553/251‑1 (MQ=255)
                                                                     cgagccctttCCTGCATTTGATTGGGCACGGTGTCGTGGTCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCTTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCTCGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAga                                                                                                                                                                                <  2:145668/242‑2 (MQ=255)
                                                                          ccGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGt                                                                                                                                                                          <  1:373211/251‑1 (MQ=255)
                                                                           cGGTCCTGAAGTTGTTTGGGCAATGTTGCGTGTGCTGGTCAGCAAGTAGACAACATGACGCTATTTACAACCTTACTGGGGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCAGTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGtt                                                                                                                                                                         <  2:512995/251‑1 (MQ=255)
                                                                                       ttGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTccc                                                                                                                                                               <  1:314240/250‑1 (MQ=255)
                                                                                              ggCACGGTTTCGTGGGCTGGGCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTGCTCTCGGGCTCACGTTAGGCATGTCCATTATTGCGATGGGCGTCACTTTTTTa                                                                                                                                                                                                        <  2:83962/204‑1 (MQ=255)
                                                                                                 cactgtgtcggttgctgtgcAGCAAGGAGACAACATGATGTTAGGTACAACCTTACTGGTGTTAATTTTCGAGCGCCCGTTTAAGTTGGGCGATCACTGGCAGCTTGATCATCTTCTTGAAGCGTTTTTTCGGCGGGGGAAATATTTTTCTCTCGGTCGTACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTGACTGTTACGCGCAGTGCAGGGAGTATTGTTCAACGTTCCCACGCTACttg                                                                                                                                                     <  2:54092/232‑3 (MQ=255)
                                                                                                   ggTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCgg                                                                                                                                                                                                                                                       <  1:395404/153‑1 (MQ=255)
                                                                                                        cGTGTGCTGGTCAGTAATGATACAACATGACGTCATTTCCAACCTTACTTGTGTTTATGTTCGATCGCCGTGTTAAGTTGTGTGCTCTCTGGCAGCTTGATCCTCGGCTTGAAGCGTTATTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACTGTAGGCATTACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCTCATTGCAGGGTTCATGGTTCAACGTTCCCACGCTACTGGTGTggtt                                                                                                                                              <  2:523858/250‑3 (MQ=255)
                                                                                                           gtgCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGgctgct                                                                                                                                           <  2:163526/250‑1 (MQ=255)
                                                                                                             gCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTTTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGATCCTGGGGTGGCTGCTGAt                                                                                                                                        >  1:177862/1‑251 (MQ=255)
                                                                                                             gCTGGGCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCGCTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTTACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGAt                                                                                                                                        <  1:551716/251‑1 (MQ=255)
                                                                                                                           gagaCAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTgtgt                                                                                                                           <  1:110217/251‑1 (MQ=255)
                                                                                                                            agacaacagggagtattttACAACCGTACTGGTGTTAATTTTCGAGGTCCTGTTTAAGTTGGGCGAGCACTGGCAGTTTGATCGACGTCTTGAAGCGTTCTTTCGGCGTGTGAAACATTTTTCTCTCGGGCGCTCGTTAGGCATGACCATTATTGCGATGGGCGTGCCTTTTTTACTGTTACGCGCATGGCAGGGAGTAGTGTTCAACGTTCCCACGCTACTGGTGTGGTTTCTGATTGGTTTGCTGTGTa                                                                                                                          <  2:475492/235‑1 (MQ=255)
                                                                                                                             gaCAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTTCGATGGGCGGGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTTTTCAACGTTCCCACGCTACTGGGGTGGCTGCTGATTGGTTTGCTGTGTAt                                                                                                                         >  1:219263/1‑251 (MQ=255)
                                                                                                                                 acaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTTTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTTCTGTTACGCGCCTTGCAGGGGGGATTGTTCAACGTTCCCCCGCTACTGGTGTGGGTGCTGTTTGGGTTTCTGTGTAtttgg                                                                                                                     >  1:303672/1‑248 (MQ=255)
                                                                                                                                     gACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAg                                                                                                                 <  1:177981/251‑1 (MQ=255)
                                                                                                                                        gCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTa                                                                                                              <  1:115836/251‑1 (MQ=255)
                                                                                                                                            tttACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGGGCCCGTTAGGGACGACCATTATTGCGGTGGGCGGGGCTTTTTTACTGTTAAGCGCATTTGCGGGAGGATTGTTCAACGTTCCCCCGCCACTGGGGGGGCTGCTGTTTGGTTTGGTGTGGATTGGCGCAGGTAAAGt                                                                                                          >  1:397783/1‑251 (MQ=255)
                                                                                                                                             ttACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACCATTTTTTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGGGGCTTTTTTTATGTTTCGCGCATTGCAGGGAGTATTTGTTAACGTTCCCACGCTACTGGTGGGGCTGCTGGTTGGTTTGCTGGGTGTTGGCGCAGGTAAAGtt                                                                                                         >  1:178927/1‑251 (MQ=255)
                                                                                                                                                       aCTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCAt                                                                                                <  1:7595/250‑1 (MQ=255)
                                                                                                                                                         tGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGGATTGTTCAAAGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGGAAAGTTCGTCTTtattat                                                                                             >  1:502578/1‑251 (MQ=255)
                                                                                                                                                               tAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGc                                                                                        <  1:512030/250‑1 (MQ=255)
                                                                                                                                                                         gcgcCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTATGCATGACCATTATTTTGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTTTTCAACCTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTTCTGTGTATTTGCGCAAGTAAAGTTCTTCTTCCTTATCATGCTTATCTGAc                                                                               >  1:302634/1‑249 (MQ=255)
                                                                                                                                                                          cgcCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAgc                                                                            <  1:439426/251‑1 (MQ=255)
                                                                                                                                                                                             gAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGGAACATTTTTTTCTCGGGCGCCCGCTAGGCATGACCCTTATTTCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTTAACGTTCCCACGCTACTGGTGTGGCTGCTGGTTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCTTTTTCATTCTTTTCTTACAACTACTTCACGAAATGAt                                                            >  1:480074/1‑249 (MQ=255)
                                                                                                                                                                                                cACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCACATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACGGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTGGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCa                                                       <  1:473888/251‑1 (MQ=255)
                                                                                                                                                                                                 aCTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCAt                                                      <  1:61268/251‑1 (MQ=255)
                                                                                                                                                                                                    ggCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGGAACATTTTTTTCTCCGGCGCCCGGTAGGGATTACCCTTATTGCGGTGGGCGTGACTTTTTTTCTTTTTCGCGCCTTGCCGGGGGGATTTTTCAACGGTTCCACGCTACTTGTGGGGCTGCTGATTGGTTTGGTGTGTATTGGCGCAGGTAAAGTTTGTCTTCCTTATTTTTTTTTTCTGGCAGCGTCTTCACGTAATGATAACCAttcc                                                   >  2:290362/1‑248 (MQ=255)
                                                                                                                                                                                                    ggCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGTGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTTACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGcc                                                   <  1:549684/251‑1 (MQ=255)
                                                                                                                                                                                                             atcatcGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTACTGTGTATTGGCGCAGGTTAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACg                                          <  1:180894/251‑1 (MQ=255)
                                                                                                                                                                                                                             gTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCGCGCTAGTCATTACACTTATTGGGATGGGAGTGTATTTTTTACTGGTACGCGCATTGCAGGGTGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCAGGTTGGGTTGTTGTGGTTTGGGGCAGGGAAAGTTCGGCTTCGTTTTCATGCTTATCTGGCAGCTGCTTCACGTAATGATACGCATGCCCGTGTCAAGAtggctggaggactc                           >  2:166682/1‑243 (MQ=255)
                                                                                                                                                                                                                                    gtgcggGTGTAACATTTTTCTCGTGGGTGCAAGTTAGTCATGACACTTATTGCGATGGTCGTGACTTTTTTTCTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCCCGCTACTGTTGTGGCTGCTGATTGGTTGGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCAGGCTTATCTTACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGAt                   <  2:249829/249‑1 (MQ=255)
                                                                                                                                                                                                                                        gggTTAAACATTTTTCGCTCTGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTAGTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGGTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCACTATCATGCTGATCTGGCAGCTGCTTCACGGAATGATAGCCATGCCCGTGCCACTATGGCTGGCGATCTCACCAGGATTCAc               <  2:513249/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACgg             >  1:20782/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      ttCTCTCGGGAGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGACATGCCCGTGCCACGATGGCTGGCGAACTCCCCCTGATTCACGGCGTCCCGGccg  >  2:523107/1‑249 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           
GATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCGGCAG  >  NC_000913/119152‑119644

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: