Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913119,3960CT27.1% 121.0 / 18.5 70S39F (TCT→TTT) ampEprotein AmpE
Reads supporting (aligned to +/- strand):  ref base C (12/39);  new base T (19/0);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 2.22e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCG  >  NC_000913/119152‑119640
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
gattgattGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCg                                                                                                                                                                                                                                                      <  1:407713/250‑1 (MQ=255)
    gattgCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGgcgc                                                                                                                                                                                                                                                 <  1:569386/251‑1 (MQ=255)
       tgCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGCTCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGCGCGCCTGTTTATGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAAAAGTTTTTTCTCGGGCGCACg                                                                                                                                                                                                                                              >  2:355632/1‑251 (MQ=255)
        gCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGCGGAACATTTTTTTTTCgggggcgcgt                                                                                                                                                                                                                                             >  2:69472/1‑244 (MQ=255)
           aTCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTGCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCGGTTTAAGTTGTGCGAGCACTGGCAGCTTGATCAGCGTCTTGAAGCGTTCTTTCGACGGGTGAAACATTTTTCTCTCGGGCGCACGtt                                                                                                                                                                                                                                            <  2:553028/249‑1 (MQ=255)
               ggATATCGCTAAAAACATGACGGGCCATTGTGATATTGCTCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTTGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCTGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAt                                                                                                                                                                                                                                     <  1:80606/251‑1 (MQ=255)
                   atCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGAc                                                                                                                                                                                                                                  <  1:166980/250‑1 (MQ=255)
                      gCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCatt                                                                                                                                                                                                                              <  1:25689/251‑1 (MQ=255)
                       cTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTTTTGGGGGCACGTTAGGCCTGAcccttt                                                                                                                                                                                                                             >  1:165861/1‑247 (MQ=255)
                         aaaaaCATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTGGAAGCGTTCTTTCGGCGGGGGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCattatt                                                                                                                                                                                                                           <  1:93970/251‑1 (MQ=255)
                         aaaaaCAGGTCGTGCCATTGTGATATTTCGCCTGATCGGAAAACCGAGCCCGGTCCTGCATGTGATTAGGCACGGTTTCTTGTGCTGGTCAGCAAGTAGACAACATGACGCTATTTACAAGCTTTCTGGTTTTAATGTTCGATCTCCTTTTTAAGATGCGCGAGCACTGGCAGCTTGAGCATCGGCTGGAAGGGTTGTTTCGGCTGGTTAAGCATTTTTCTCTCGGGCGCACGTTAGGCATGACCattat                                                                                                                                                                                                                            <  1:242285/250‑1 (MQ=255)
                            aaCATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGGCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGTTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGGAGCGTTCTTTCGGCGGGTGAAACATGTTTCTCTCGGGCTCACGTTAGGCATGACCATTATTGCg                                                                                                                                                                                                                        <  2:547031/251‑1 (MQ=255)
                             ccATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCgg                                                                                                                                                                                                                                                     <  2:17441/221‑1 (MQ=255)
                              cATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAATATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGAt                                                                                                                                                                                                                     <  1:187514/251‑1 (MQ=255)
                                      gCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTAGTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGATGCGTTCTTTCGGCGGGTGAGACACTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGa                                                                                                                                                                                                             <  1:534603/251‑1 (MQ=255)
                                         aTTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGCAACATTTTTTTCTCGCGCGCCCGCTATGCATGCCCATTATTTCGCTGGGGGGGACtt                                                                                                                                                                                                          >  2:264599/1‑251 (MQ=255)
                                               atatTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGGAAACTTTTTCTTTCCGGCGGGCGTTTGGCCTGACCATTATTGCGCTGGGCGTGGCTTTTTTat                                                                                                                                                                                                    >  1:214580/1‑250 (MQ=255)
                                                    gcgcCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTa                                                                                                                                                                                              <  1:432258/251‑1 (MQ=255)
                                                    gagcCGGATCGGAAAACCGAGCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTGCTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGtt                                                                                                                                                                                               <  2:5414/248‑1 (MQ=255)
                                                             cGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTg                                                                                                                                                                                     <  1:355632/251‑1 (MQ=255)
                                                              ggAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGc                                                                                                                                                                                    <  1:457783/251‑1 (MQ=255)
                                                              ggAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGCCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGAGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTCTTTTGGGGGCCCGTTATGCCAGGACCTTTTTGCGTTGGGCGGGGCTTTTTTTCTGTTTAGCGCATTGc                                                                                                                                                                                    >  1:365729/1‑251 (MQ=255)
                                                                  aaCCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCCCGTTAAGCCAGACCATTATTGTGCTGGGGGGGGCTTTTTTACTGTTACGCGCATTGCAggg                                                                                                                                                                                >  1:321504/1‑251 (MQ=255)
                                                                      gATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTAc                                                                                                                                                                                                    <  2:141135/228‑1 (MQ=255)
                                                                      gATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCCTGACCCTTATTGCGATGGGCGGGGCTTTTTTAc                                                                                                                                                                                                    >  1:141135/1‑228 (MQ=255)
                                                                         cccGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCGGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTGACTGGTGTTAAGTTTCGAGCGCCTGGTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTAGTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTg                                                                                                                                                                         <  1:341514/251‑1 (MQ=255)
                                                                          ccGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAAAATTTTTCTCTCGGGCGGACGTTAGGTCTGACCATTTTTGCGATGGGCGGGCCTTTTTTACTGTTTCGCGCATTGCCGGGGGAATTGt                                                                                                                                                                        >  1:302526/1‑251 (MQ=255)
                                                                               ccTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCTGGTGAAACATTTTTCTCTCGGGTGCACTGTAGGCAGGACCATTATTGCGATTGGCTTGACTTTATTACTGTTACGCGCATTGCAGGGAGTATTGTTCAAc                                                                                                                                                                  <  1:50158/251‑1 (MQ=255)
                                                                                      ttGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTccc                                                                                                                                                            <  1:476688/250‑1 (MQ=255)
                                                                                           tGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCCCGTTAGGCATGACCATTATTGCGATGGGGGGGGCTTTTTTACTGTTACGCCCATTGCCGGGGGTATTGTTCAACGTTCCCACGCtc                                                                                                                                                      >  1:204847/1‑250 (MQ=255)
                                                                                                aCGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGGACCATTTTTTTCTCTGGCGGCCGGTAAGCATGGACATTTTTGCGGTGGGCGTGACTTTTTTTTTGGTTAGCGCATTGCAGGGGGTATTTTTCAACGTTcccccggtcctggg                                                                                                                                                 >  1:503709/1‑240 (MQ=255)
                                                                                                      tCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGGGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATGGTTCAACGTTCCCACGCTACTGGTGTGgc                                                                                                                                            <  2:165861/250‑1 (MQ=255)
                                                                                                               ggTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTTTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTgg                                                                                                                                  >  1:412724/1‑251 (MQ=255)
                                                                                                                gTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGCTCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTTTCTCGGGCGCCCGTTAGGCCTGACCATTATTGCGATGGGCGTGGCTTTTTTTCTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTTCCCCGCCACTGGTGTGGCTGCTGAtttgt                                                                                                                                 >  1:73418/1‑248 (MQ=255)
                                                                                                                     cAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTAGAGCGCCTGTTTAAGTTGGGCGAGCACTGACAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGCAACATTTTTTTCTCAGGCGCACCTTAGGGACGGCCATTATTGCGATGGGCGGGGCTTTTTTATTGTTTCGAGCATTGCAGGGGGGATTGTTTCAACTTTCCACACCAATGGTGTGgctggtgctttggttgct                                                                                                                            >  1:565691/1‑237 (MQ=255)
                                                                                                                           aacaacaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTGCGGCGGGTGAAACATTTTTCTCTCGGGCGCTCGTTAGGCATGACCATTATGGCGATGGGCGTGACTTTTTTGCTGTTACGCGCATTGCAGGGAGTATGGTTCAACGTGCCCACGCTACTGGTGTGGCTGCTGATGGGTTTGCTGTGTAt                                                                                                                     <  2:503709/250‑1 (MQ=255)
                                                                                                                            acaacaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCGGATTGGTTTGCGGTGTAtt                                                                                                                    <  1:443622/251‑1 (MQ=255)
                                                                                                                                       cTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGGCTTTTTTACTGTTTCGCGCATTGCAGGGAGGATTGTTCAAAGTTCCCACGCTACTGGGGTGGCTGCTGATTGGTTTTCTGTGTATTGGCGCAGGTa                                                                                                          >  1:325253/1‑250 (MQ=255)
                                                                                                                                             aCAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGTAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACTCGCATTGCAGGGAGTATTGGGCAACGTTCCCACGCTACTGGTGTGGCGGCTGATTGGTGTGCTGGGTAGTGGCGCATGTAAAGTTCg                                                                                                   <  2:487644/251‑1 (MQ=255)
                                                                                                                                              cAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCAATGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAAACTTTTTCTCTCCGGCGCACGGTTGGTATGACCATTATTTCGATGGGCGAGAATTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGCTACCACGCTACTGGTGTGGGTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGttagt                                                                                                  >  1:290347/1‑248 (MQ=255)
                                                                                                                                                aCCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGGAAAATTTTTCTCTCCGGGGGCCGGTTGGCATTGCCATTTTTTCGGATGGCGTGACTTTTTTTCTGTTAAGCGCCTTGTAGGGAGGATTTTTCAACGTTCCCACGCTACTGGGGGGGCTTGTGATTGGGTTTCTGGGTATTGGTGCAGGGAAAgatacttt                                                                                                >  1:15231/1‑244 (MQ=255)
                                                                                                                                                aCCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTTTCTCGGGCGCACGTTAGGCATGACCCTTTTTGCGGTGGGCGTGACTTTTTTTCTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCCCGCTACTGGTGGGGCTGCTGTTTTGTTTGCTGTGTATTGGCGCAGGGAAAGTTCGTCt                                                                                                >  1:460594/1‑251 (MQ=255)
                                                                                                                                                 cctgaCTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCtt                                                                                               <  2:426230/247‑1 (MQ=255)
                                                                                                                                                  cTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGAGAAACATTTTTTTCTCGGGCGCACGGTAGGCATGACCCTTATTGCGATGGGCGTGACTTTTTTTCTTTTTCGCGCATTGCAGGGGGTATTGTTCAACGTTCCCACGCCACTGGTGTGGCTGCTGGTTGGTTTTCTTTGTTTTTGCGCAGGTAAAGTTCGTCttt                                                                                              >  1:236513/1‑250 (MQ=255)
                                                                                                                                                     aCTGGTGTTAATTTTCGAGCTCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCGTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTTTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGAGTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCAtt                                                                                           <  1:523392/251‑1 (MQ=255)
                                                                                                                                                        ggTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATc                                                                                        <  1:454251/251‑1 (MQ=255)
                                                                                                                                                        ggTGTTAATGTTCGAGCGCCTGTTGAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCTCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTGACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTGCCCACGCTACTTGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGGCTTCATTa                                                                                          <  2:550933/249‑1 (MQ=255)
                                                                                                                                                                ttttCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGGAACATTTTTTTCTCCGGCGCACGGTAGGCGTGGCCCCTATTTCGGAGGGGGTGGCTTTTTTTCTGTTACGCGCCTTGCAGGGGGGATTTTTTAACCTTCCCCCGCTACTGGGGGGGGTGCTGATTGGTTTTTTGTGGTTTTGCGCCGGGAAAGTTCGTCTTctttttatttctttt                                                                                >  2:366406/1‑237 (MQ=255)
                                                                                                                                                                  ttCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTTCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATc                                                                               <  1:69472/250‑1 (MQ=255)
                                                                                                                                                                    cGAGCGCCTGGTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGGCTTCATTATCATGCTTATCTg                                                                             <  2:302526/250‑1 (MQ=255)
                                                                                                                                                                     gAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGCAGCGTTCTTTCGGCGGGGGAAACCTTTTTTTCTCGGGCGCACGCTAGGCATGACCATTTTTGCGATGGGAGGGACTTTTTTACTGTTACGCGCATTGCAGGGAGGATTTTTCAAAGTTCCCACGCTAATGGGGTGGCTGCTGATTGGTTTTCTGTGTATTGGGGCCGGTAAAGTTTGTCTTTATTTTTATGCTTATCttac                                                                           >  1:54045/1‑248 (MQ=255)
                                                                                                                                                                     gAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGgc                                                                                                                 <  2:14214/213‑1 (MQ=255)
                                                                                                                                                                           cTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACCGTAGGCATGACCATTATTGCGGCGGGCGGGACTTTTTTACTGTTACGCGCATTGCAGGGGGGATTGTTCAACGTTCCCACGCTACTGGGGTGGCTGCTGGTTGGGTTGCTGGGTATTGGGGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAgctgc                                                                     >  1:19008/1‑251 (MQ=255)
                                                                                                                                                                           cTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAAAATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGGGTGGGTGCTGATTGGTTTGCTGTGTATTGGCGCAGGGAAAGTTCGTCTTTCTTATCATGCTTATCTGCAAgctgc                                                                     >  1:447292/1‑251 (MQ=255)
                                                                                                                                                                                tAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATGTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCAGATTGGTTTGCTGTGTATTTGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCa                                                                 <  1:515951/250‑1 (MQ=255)
                                                                                                                                                                                 aaGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCAc                                                                <  1:496874/250‑1 (MQ=255)
                                                                                                                                                                                   gTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTTTCTCGGGCGCACCGTAGGCATGACCATTATTTCGAGGGGCGTGACTTTTTTTCTGTTTCGCGCATTGCAGGGAGGATTTTTCAACGTTTCCACGCTACTGGTGGGGCTGCTGATTGGTTTGCTGTGTGTTTGCGCCGGGAAAGTTCGTCTTCATTATCCTGCTTATCTTACAGCTGCTTCCCGt                                                              >  1:316783/1‑250 (MQ=255)
                                                                                                                                                                                   gTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCGCGGGCGCACGTTAGGCATGACCATTAGTGCGAGGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGATTATTGTTCAACGTTCCCACGCTACTGGTGTGgctgct                                                                                                                                        <  1:36002/177‑1 (MQ=255)
                                                                                                                                                                                           aGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGGAAAATTTTTTTTTCGGGCGCACGTTAGGGCTTACCCTTTTTGTGGTGGGGGTGGCTTTTTTTCTGTTACGCGCATTGCAGGGGGTATTTTTCAACGTTTCCACGCTACTGGTGTGGGTGCTGATTGGGTTGCTGTGGATTGGGGCCGGGAAAGGTCGTCTTTATTATTCTGCTTTTCTggtaggtgcttcccggtatgat                                                        >  2:59773/1‑227 (MQ=255)
                                                                                                                                                                                           aGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTTCTGTTACGCGCATTGCAGGGAGGATTGTTCAACGTTCCCCCGCTACTGGTGTGGCTGCTGATTGGGTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGc                                                     >  1:339809/1‑251 (MQ=255)
                                                                                                                                                                                            gCACTGGCAGCTTTATCATCGTCCTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCCATGGGCGTGACCTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGcc                                                    <  1:418162/251‑1 (MQ=255)
                                                                                                                                                                                                tGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTCTCTCGGGCGCACGGTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCTGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGTTTGGTTTGCTGTGTATTGGCGCAGGTAACGTTCGTCTTCATTTTCATGCTTATCTGACAGCTGGTTCACGTAATGATAGCCATGc                                                >  1:416384/1‑251 (MQ=255)
                                                                                                                                                                                                    aGCTTGATCATCGGCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGt                                            <  2:339809/251‑1 (MQ=255)
                                                                                                                                                                                                          atcatcGTCTTGAAGCGTTCTTTCGGCGGGGGGAACATTTTTCTTTCGGGCGCCCGGTAGGCATGACCCTTATTGCGATGGGCGGGGCTTTTTTACTGTTACGCGCAATTCAGGGAGTATTTTTTAACGTTTCCACGCTACTGGTGTGGGTGCTGGTTGGTTTTCTGGTGTTTGGCGCAGGGAAAGTTCGTCTTCATTTTTCTGCTTTTTTGACAGCTGCTTCACGGAATGATAGGCATCTCCGGGCCACg                                      >  2:394221/1‑251 (MQ=255)
                                                                                                                                                                                                                       aGCGGTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAAc                         <  2:204847/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cggcggGTGAAAACTTTTTCTCTCGGGCGCACGTTAGGCCTGACCATTATTGCGGTGGGCGTGACTTTTTTTCTGTTTCGCGCATTGCAGGGGGTATTGTTCAACGTTCCCACGCTACTGGTGTGGGTGCTGATTGGTTTGCTGTGGATTTGCGGCGGGAAAGTTCGTCTTCCTTATCATTCTTATCTGACCGCTGCTTCACGGAATGATAACCATGCCCGTGCCCCGACGGCTGGCggaaccacccctat               >  2:262993/1‑238 (MQ=255)
                                                                                                                                                                                                                                  ggcggGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGAtt              <  1:366406/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggcggGGGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCAGTCTAGCGATGGGCGTGACTTTCTTAATGTTACGCGCATTGCAGGGCGCATTGTTCAACGTTCCCACGCTACTGGTGTGGATGCTGGTTGGATTGCTGTGTATTCGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGATGCTTCACGTAATGATAGCCATGACCGTGCCACGTTGGCTGGCGAACCCACCATGAtt              >  1:349005/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGCCTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCCGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAACGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCGCCATGATTCACGGCGTccc   >  1:299398/1‑251 (MQ=255)
                                                                                                                                                                                                                                               agttttCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCAGTGCAGGGAGTAGTGTGCAACGTTCCCACGCTACTGGTGTTGCTGCTGAGTGGGTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCGCGTAATGATAGCCTTTCCCGTGCCACGAGAGCTGGCGAACGCACCATGATTCACGGCGTCCCg  <  2:227378/249‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
GATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCG  >  NC_000913/119152‑119640

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