Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913119,3960CT24.1% 105.7 / 15.1 59S39F (TCT→TTT) ampEprotein AmpE
Reads supporting (aligned to +/- strand):  ref base C (7/37);  new base T (14/0);  total (22/37)
Fisher's exact test for biased strand distribution p-value = 1.15e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCGGCA  >  NC_000913/119152‑119643
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
gattgattGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCAGCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTtctctc                                                                                                                                                                                                                                                        <  1:606739/249‑1 (MQ=255)
         cTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACCTTTTTTTTTCgggggcccggt                                                                                                                                                                                                                                             >  1:588195/1‑243 (MQ=255)
           aTCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑GAAAATTTTTTTCTCGGGCGCACG‑tt                                                                                                                                                                                                                                            >  1:37099/1‑249 (MQ=255)
            tCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAAAATTTTTTTCTCGGGCGGACG‑TTAgg                                                                                                                                                                                                                                         >  1:785953/1‑251 (MQ=255)
             ccGGATATCGCATAAAAGATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCGATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCGTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGc                                                                                                                                                                                                                                        <  2:6752/251‑1 (MQ=255)
               ggATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCAt                                                                                                                                                                                                                                      <  1:305840/251‑1 (MQ=255)
                   accGCTAAAAACATGACAGGCCATTGTTATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGAc                                                                                                                                                                                                                                   <  2:786410/248‑1 (MQ=255)
                           aaaCATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑GAACATTTTTCTCTCGGGCGCACG‑GTAGGccataccactcttgc                                                                                                                                                                                                                          >  1:35127/1‑237 (MQ=255)
                                  aCGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACATTTTTTTCTCGGGCGCAAG‑TTAGGCATGGACattatttgaatgggg                                                                                                                                                                                                                   >  1:549106/1‑242 (MQ=255)
                                       ccATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGAc                                                                                                                                                                                                             <  1:715219/251‑1 (MQ=255)
                                                   ggcgtcGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTGTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCAGTATTGC‑GATGGGCGTGACTTTTTTACTGtt                                                                                                                                                                                                 <  2:117335/246‑1 (MQ=255)
                                                         ggATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGCAGCGTTCTTTCGGCGGGGGCAAAC‑TTTTTTTCTCGGGGGCCCGGTTTGCCA‑GACCATTAATCC‑CATGGGCGGGACTTTTTTACtttttaccgc                                                                                                                                                                                           >  1:137475/1‑242 (MQ=255)
                                                          gATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑tt                                                                                                                                                                                                                                            <  1:598588/202‑1 (MQ=255)
                                                             cGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTg                                                                                                                                                                                       <  1:693123/251‑1 (MQ=255)
                                                                     cGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGCGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTTTCTCGGGCGCACG‑GTAGGCATGACCCTTATTTA‑GATGGGGGGGACTTTTTTTCTGTTAAgcgcctttaagggagg                                                                                                                                                                               >  1:460992/1‑239 (MQ=255)
                                                                         cccGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑g                                                                                                                                                                          <  1:470119/251‑1 (MQ=255)
                                                                                 tGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACCACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑GAACATTTTTTTCTCGGGCGCCCG‑TTCGGCATTGCCATTATTTCAGTTGGG‑GGGACTTTTTTTCTGTTTCCCGCCTTGCCGGGAGGATT‑GTTTAACGt                                                                                                                                                                  >  2:734024/1‑251 (MQ=255)
                                                                                    aTTTGATTGGGCACGGTTTCGTTTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTcc                                                                                                                                                               <  1:507280/251‑1 (MQ=255)
                                                                                                       cGTGTGTTTGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TGAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGgct                                                                                                                                             <  1:619958/250‑1 (MQ=255)
                                                                                                        gtgtGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACATTTTTTTCTCGGGCGCCCG‑TTAGGCATGACCATTATTTC‑GATGGGCGGGACTTTTTTACTTTTACGCGCATTGCAGGGGGTATT‑GTTCAAAGTTTCCACGCTACTGGTGTGggtgc                                                                                                                                           >  1:484782/1‑251 (MQ=255)
                                                                                                        gtgtGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAAAATTTTTTTCTCGGGCGCACG‑TTAGGCATGGCCCTTTTTGC‑GGTGGGGGGGGCTTTTTTTCTTTTACGCGCACTTCAGGGAGGATTTGTTCA‑CGTTCCCCCGCTACTTGTGGGgctgg                                                                                                                                           >  2:250944/1‑250 (MQ=255)
                                                                                                          gtgCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGctgctg                                                                                                                                         <  1:192412/251‑1 (MQ=255)
                                                                                                            gCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGa                                                                                                                                        <  1:60620/250‑1 (MQ=255)
                                                                                                            gCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGAt                                                                                                                                       <  1:641233/251‑1 (MQ=255)
                                                                                                             cTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCACTATTGC‑GATGGGCGTGGCTTTTTTTCTGTTACGCGCATTGCAGGGGGGATT‑GTTCAACGTTCCCACGCTACTGGTGTGGctgctgttt                                                                                                                                      >  1:241883/1‑248 (MQ=255)
                                                                                                                  cAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGtt                                                                                                                                  <  1:784657/250‑1 (MQ=255)
                                                                                                                      aaGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑Gt                                                                                                                                                                         <  1:587236/207‑1 (MQ=255)
                                                                                                                        ggAGACAACATGACGCTATTTACACCCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCtgtg                                                                                                                           <  1:426350/251‑1 (MQ=255)
                                                                                                                             caacaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTg                                                                                                                      <  1:629003/251‑1 (MQ=255)
                                                                                                                             caacaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAAAATTTTTTTCTCGGGCGCACG‑TTTGGCATGGCCATTATTTC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGGATT‑GTTTAAAGTTCCCACGCTACTGGTGTGGCTGCTGTTTTGTTTGCTgtgtgttt                                                                                                                      >  1:172547/1‑247 (MQ=255)
                                                                                                                               acaTGACGCTAGTGACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGgc                                                                                                                    <  1:798497/251‑1 (MQ=255)
                                                                                                                                    actcTATTTAGAACCTGACTGGTGTTAATTGTCGAGCTCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTGTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCAATATTGC‑GATGGGCGTTACTTGTTCACTGTTACGCGCATTACAGGGAGTGTT‑GTTCAACGTTGCCACGCTACGGGTGTGGCTGCTGATTGTTTTGCTGTGTATTGGCGCAgg                                                                                                               <  2:95283/248‑1 (MQ=255)
                                                                                                                                    aCGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAgg                                                                                                               <  2:484782/251‑1 (MQ=255)
                                                                                                                                       cTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAAAATTTTTTTCTCGGGGGCACG‑GTAGGCATGACCATTATTTC‑GATGGGGGGGGCTTTTTTACTGTTACGCGCATTTCCGGGGGTATT‑GTTTAACGTTTCCAAGCTACTGGGGGGGCTGCTGGTTGGTTTGCTGTGTATTTGCGCCGGTa                                                                                                             >  1:568154/1‑250 (MQ=255)
                                                                                                                                         aatctcAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGGTTGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATTGGCGTGACTTTTTTACTTTTACGCGCATTGCAGGGAGTATT‑GTGCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTaaa                                                                                                           <  2:372323/245‑1 (MQ=255)
                                                                                                                                                aCCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCt                                                                                                   <  2:613761/251‑1 (MQ=255)
                                                                                                                                                  cTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGTGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATTAACCTTATTTC‑GATGGGCGGGACTTTTTTACTGTTACGCGCCTTGCAGGGAGTATT‑TGTCAACGTTTCCCCTCTACTGGGGTGGCTGCTGATTGGTTTGCTGGGTATTGGGCCCGGTAAAGTTCGTCttt                                                                                                 >  1:652362/1‑250 (MQ=255)
                                                                                                                                                      cTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTTTCTCCGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGGGACTTTTTTTCTGTTACGCGCCTTGCCGGGGGGATT‑TTTCAACGTTTCCACGCTAATGGGGTGGCCGCTGATTGGTTTGCTGTGTATTTGCGCAGGGAAAGTTCGTCTTCAttt                                                                                             >  1:468385/1‑250 (MQ=255)
                                                                                                                                                        ggTGTTAATGTTCGAGCGCCTGTTTAAGTTGGGCGAGCTGTGGCAGCTTGATCATCGTTTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTAt                                                                                            <  2:25106/250‑1 (MQ=255)
                                                                                                                                                         gttgtAATTTTCGAGCGCCTGTTTAAGTTGGGCGGTCGCTTGTCGCGTGATCATCGTCCTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTTCGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCa                                                                                          <  2:501221/247‑1 (MQ=255)
                                                                                                                                                            tgaaTTTTCGAGCGAGTGTTTAAGTTTGGCGAGCACTGGTAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AACCATTTTTCGCTCGGGCGCACG‑TTAGGCATGACCATTATGGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCACACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATg                                                                                        <  2:137475/248‑1 (MQ=255)
                                                                                                                                                             tAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGc                                                                                       <  1:255832/250‑1 (MQ=255)
                                                                                                                                                               aTTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACATTTTTTTCTCGGGCGCACG‑TTTGGCATGACCCTTATTGC‑GATGGGCGGGACTTTTTTTCTTTTACGCGCATTGCAGGGAGGATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGGTTGGTTTTCTGTGTATTGGCGCAGGTAAAGTTCGTCTTTATTATCATTCttt                                                                                    >  1:196002/1‑250 (MQ=255)
                                                                                                                                                                       gcgcCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAg                                                                            <  2:234167/251‑1 (MQ=255)
                                                                                                                                                                         gcCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACATTTTTTTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTTCTGTTACGCGCATTGCAGGGGGGATT‑GTTTAACGGTCCCACGCTTCTGGTGTGGCTGCTGGTTGGTTTTCTGTGGTTTGGGGCCGGTAAAGGTCGTCTTCCTTATCATggttttcttacaact                                                                          >  1:137407/1‑237 (MQ=255)
                                                                                                                                                                                tAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCa                                                                    <  2:380670/250‑1 (MQ=255)
                                                                                                                                                                                tAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCa                                                                    <  1:367692/250‑1 (MQ=255)
                                                                                                                                                                                 aaGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCAc                                                                   <  2:373727/250‑1 (MQ=255)
                                                                                                                                                                                          gAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAg                                                         <  2:37099/251‑1 (MQ=255)
                                                                                                                                                                                          gAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAAAATTTTTATCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGGGGATT‑TTTCCACGTTCCCACGCTACTGGGGTGGCTGCTGATTGGTTTTCTGTGTATTGGGGCAGGGAAAGTTCGTCTTCATTTTAATGCTTTTTTGACAGCTGCtttccggaattata                                                          >  1:349667/1‑238 (MQ=255)
                                                                                                                                                                                             cACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGG‑AAACATTTTTCTCTCGGGCGCACG‑TTAAGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTTCTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGGTTGGTTTGCTGTGTATTGGCGCCGGGAAAAGTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGccc                                                      >  1:799159/1‑250 (MQ=255)
                                                                                                                                                                                                               cgtattGAAGCGTTCTTTCGGCGGGTG‑AAACAGTTTTCTCGCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑TTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGAtgg                                     <  2:453471/246‑1 (MQ=255)
                                                                                                                                                                                                                gTCTGGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGAtggct                                   <  2:785953/251‑1 (MQ=255)
                                                                                                                                                                                                                      aagcggtttttCGGCGGTTT‑AAACATTTTTCTCTCGGTCGCACG‑TTAGGCATGACCATTAGTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATTGCTGGCGaa                             <  2:519261/243‑1 (MQ=255)
                                                                                                                                                                                                                                tCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATg                    <  1:20502/250‑1 (MQ=255)
                                                                                                                                                                                                                                         g‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCCCCATGATTCACGGcc          >  2:364926/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aaCATTTTTCTCTCGGGCGCACG‑TTAGGCATTACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTCAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGCAAGCCATGCCCGTGCCACGGTGTCTGGCGAACTCAACATGATTCCCGGCGTc        >  2:277046/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 ttttCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCgg    <  1:177087/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttttCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCGGc   >  1:536120/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ttCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCGGCa  <  2:309672/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
GATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTG‑AAACATTTTTCTCTCGGGCGCACG‑TTAGGCATGACCATTATTGC‑GATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATT‑GTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGCGTCCCGGCA  >  NC_000913/119152‑119643

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