Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,416,444 A→T 15.5% I472F (ATC→TTC)  acrF → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,416,4440AT15.5% 146.3 / 2.4 66I472F (ATC→TTC) acrFmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (9/43);  new base T (9/2);  total (21/45)
Fisher's exact test for biased strand distribution p-value = 8.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTCGGGATGGTGCTTGCCATCGGGCTGCTCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCGGATCCACAGGACGATATTTACTGATCTATGCGCT  >  NC_000913/3416214‑3416667
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                  
tGTTCGGGATGGTGCTTGCCATCGGGCTGCTCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTAAGCGGTATTTATTCCGATGGCATCCTTCGGCGGTTCTCA‑TTGGGCCATTTTTCGCCCGTTTTCTTTCACCAt                                                                                                                                                                                                                             >  2:813032/1‑238 (MQ=255)
  ttCGGGATGGTGCTTGCCATCGGGCTGCTCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGGGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TTGGGGCATTTTTCGCCAGTTCTCTTTCCCCCTCGt                                                                                                                                                                                                                          >  2:627305/1‑239 (MQ=255)
          ggTGCTTGCCATCGGGCTGCTCGTGGATGATGCGATTGTGGTGGTGGAGAACGTCGAGCGGGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCa                                                                                                                                                                                                                   <  2:515514/238‑1 (MQ=255)
          ggTGCTTGCCATCGGGCTGCTCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTTTAC‑TGGGGCCATTTATTGCCCGTTCTCTTTTTCCATCGTTTTGGCaa                                                                                                                                                                                                                  >  2:513067/1‑239 (MQ=255)
            tGCTTGCCATCGGGCTGCTCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGTTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTGTTAC‑TGGGGGAATTTTTCTCCCGTCTTCTTTCCCCACACTTTCTGCAA‑tt                                                                                                                                                                                                               >  2:715145/1‑238 (MQ=255)
            tGCTTGCCATCGGGCTGCTCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TTGGGCAATTTTTTCCCAGGTTTCTATCACCATCGGTTCGGCAA‑Tg                                                                                                                                                                                                               >  2:515593/1‑239 (MQ=255)
            tGCTTGCCATCGGGCTGCTCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGAACTGGCGGGGATCGCGATGGGGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGTCGGTTCTAC‑TGTGGGAATTTATCGCCCGGTCTCTCTCACCCTCGCTTCGGCAA‑Tg                                                                                                                                                                                                               >  2:511962/1‑239 (MQ=255)
               ttGCCATCGGGCTGCTCGTCGATGGTGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGgc                                                                                                                                                                                                             <  2:642387/238‑1 (MQ=255)
                tGCCATCGGGCTGCTCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGGGCTGTCAGCGGTATTTATTCCGATGGCATTTTTCGCCGGTTCTTC‑TGGGGCAACTTTTCTCCCGGTCTCTTTCACCCACGTTTCGGCAA‑GGgc                                                                                                                                                                                                             >  2:349129/1‑237 (MQ=255)
                      cGGGCTGCTCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGTTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TTGGGGAATTTATCGCCAGGTTTCTACCTCCCTCGTTTCGGCAA‑TGGGGCCttcgg                                                                                                                                                                                                     >  2:214895/1‑237 (MQ=255)
                      cGGGCTGCTCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TTGGGCAATTTATCGCCAGTTCTCTTTTACCATCGTTTTGGCAA‑TGGCCCTttctt                                                                                                                                                                                                     >  2:760356/1‑238 (MQ=255)
                      cGGGCTGCTCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGGGGGTATCGCGATGGGGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTACTTGGGGAATTTTTTCCCCAG‑TCTCTCTCACCACCGGTTCGGACA‑AGGCGCtt                                                                                                                                                                                                         >  2:570780/1‑235 (MQ=255)
                             ctCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTTTAC‑CGGGGGACATTTTCGCCCGTTTTTTTTCACCCCCGTTTTGGCaaaaggggcttttggtttggt                                                                                                                                                                                               >  1:185062/1‑220 (MQ=255)
                                 tggagatgggtaaaaggaccagggAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCa                                                                                                                                                                                          <  2:421165/217‑1 (MQ=255)
                                  cgATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAt                                                                                                                                                                                         <  2:746048/239‑1 (MQ=255)
                                   gttgaggCGATAGTGGGGGTGGAGAACGTCGAGCGAGGGATGAGGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAt                                                                                                                                                                                         <  2:79136/236‑1 (MQ=255)
                                   gattgagCGATAGGGGTGGGGGAGAACGTCGACCGAGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAt                                                                                                                                                                                         <  2:792468/232‑1 (MQ=255)
                                    agggtgCGATAGTGGTGGTGGAGAACGGCGAGCGCGTGATGAGGGGGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAt‑t                                                                                                                                                                                       <  2:93608/234‑1 (MQ=255)
                                     tgatgCGATAGTGGTGGTGGAAAACGTGGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑Tg                                                                                                                                                                                      <  2:603243/238‑1 (MQ=255)
                                            aTAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATTGCCCGGTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTTTTCTGGTGGCCT‑TGATTCTTa                                                                                                                                                                               >  2:204801/1‑238 (MQ=255)
                                             tggtgggggggGAAAACGTCGAGCGCGTGAGGATGGAGGATAAGCCCCCGCCAAAAGAAGCGGCGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTTGTGGGTTTCGCGGTGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTAcc                                                                                                                                                                             <  1:794450/237‑1 (MQ=255)
                                              ggtggtgggcGAGAACGTCGCGCGCGGTAGGATGGAGGATAATCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTAccc                                                                                                                                                                            <  2:722174/238‑1 (MQ=255)
                                                      tggAGAACGTCGGGCGCGGGATGATGGAGGATAAGCCCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑Gtt                                                                                                                                                                    <  2:135021/238‑1 (MQ=255)
                                                      gggAAAACGTCGAGCGCGTGAGGATGGAGGAAAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑Gtt                                                                                                                                                                    <  2:129670/237‑1 (MQ=255)
                                                             cGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGGATCGCGATGGAGCTGTCAGCGGTATTTATTCCGATGGAATTCTTCGGCGGTTCTAC‑TGGGGCAATTTTTCTCCCGTTTCTTATCCCCCTCCTTTCGGGAAGTGGCCGTTTCTTTTCT‑GTGGCAT‑TGTGTTTTAccccccgtgtttggggaa                                                                                                                                                             >  1:651967/1‑225 (MQ=255)
                                                                cgggggcgTGATGATGGAGGATAAGCCCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACg                                                                                                                                                          <  2:310440/233‑1 (MQ=255)
                                                                cggggcggtgtTGAGGGGGGCTAAGCCCCCGCCAAAAGAAGCGACGGAAAAAGCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCGGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACgc                                                                                                                                                         <  1:620799/232‑1 (MQ=255)
                                                                cggagggggtgtgatgGAGGATAAGCCCCCGCCCAAAGAAGCGAAGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACg                                                                                                                                                          <  2:839638/227‑1 (MQ=255)
                                                                cgcacgggTGAGGATGGGGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAAACGAGGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACg                                                                                                                                                          <  2:682627/231‑1 (MQ=255)
                                                                cGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCCGTTTTCTATCACCATCGTTTTGGCAA‑TGGCGCTTTCTGTTCTGGTTGCAT‑TGATTCTTACCCCCGC‑Gtttttggcaaac                                                                                                                                                          >  2:478238/1‑228 (MQ=255)
                                                                           ggttggcgataaaactctCGCCAAAAGAGGCGCCGGAAAAATTGTTGTCGCAAATTCT‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTCCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAAcc                                                                                                                                               <  1:349129/224‑1 (MQ=255)
                                                                                     tAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTCC‑GTTATGTTCAACGCTGCTTAAACCCCtccttgcgt                                                                                                                                    >  2:230836/1‑232 (MQ=255)
                                                                                       accgcccGCCAAAAGAGCGGAAGGAAAAATTGCGGCCGCAAATTCA‑GGGCGCGCTGGTGGGTATCGCGATGGGTCTGGCAGCGGTATTTATTCCGATGGCATTCCTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCTTCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAAGCCGTCTCTGCTGAg                                                                                                                                  <  1:290245/235‑1 (MQ=255)
                                                                                         ccccccccAAAAGGGGCCGCGGAAAAATCGATGTCGCAAATGCA‑GGGCGCACTGGTGGGTTCCGCGATGGTGCTTGCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCa                                                                                                                                <  1:513067/237‑1 (MQ=255)
                                                                                             cGCCAAAAGGAGCCACGGAAAAAGCGATGTCGAAAATTCA‑GGGCGCACTGGTGGGTATCGGGTTGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCAc                                                                                                                            <  1:575515/239‑1 (MQ=255)
                                                                                                   aagaagCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCt                                                                                                                                                                                 <  2:223786/181‑1 (MQ=255)
                                                                                                          cgcGGAAAAATCGGTTTCGCAAATTCAGGGGCGC‑CTTGTGGGTATCGGGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCg                                                                                                               <  1:459623/237‑1 (MQ=255)
                                                                                                           aCGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCg                                                                                                               <  2:185062/238‑1 (MQ=255)
                                                                                                                     cGATGGCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCggt                                                                                                    <  2:366720/239‑1 (MQ=255)
                                                                                                                      gATGTCGCAAATTCA‑GGGGCCACTGGTGGGTGTCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑GGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATATGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCggtt                                                                                                   <  1:541229/239‑1 (MQ=255)
                                                                                                                                   tctgggCGCCTTGGTGGGTTTCGGGATGGTGTGGTGCGGGGTATTTTTTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATacca                                                                                       <  1:515593/236‑1 (MQ=255)
                                                                                                                                   tccgggCGCACTGGTGGGTAGCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATacca                                                                                       <  2:450826/236‑1 (MQ=255)
                                                                                                                                            aCTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATTCACCGTTTTCTTTCACCCACGTTTCGGCCA‑TGGCGCTTTCTGTTCTTGTGGCATGTGTTTTTTCCCCCT‑C‑CTTCTGTTCCACGCCGCTTAAAACCCGCTTTTCTTGGACTCCCGGAAAAAAGGGGGGGGTCTTTGGGTGTTTTATTACCAACTTccgtc                                                                              >  1:824548/1‑235 (MQ=255)
                                                                                                                                                     tATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATAACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCCT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCCGCTCTTCTGAGCATCACGGAAATAAAGGCGGGTTTTTCGGTTGGTTTAATACCACCTTCCATCAAAACGTTa                                                                     >  2:647125/1‑239 (MQ=255)
                                                                                                                                                        cggtgtggtgcttgcatgtgtttttactCCCATGGCATGCTTGGGGGGTTCTGC‑GGGGGCCATTTACGTGCAGTTCTCTTTCACCTTCGGTTCGGAAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTTCCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAAcc                                                                  <  1:536674/217‑1 (MQ=255)
                                                                                                                                                        cgcgATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTGTAC‑TGGGGCAATTTTTTGCCAGTTCTCTATCACCATCGTTTCGGCCA‑TGGCGCTTTTTTTTTTGGTGGAGT‑TGATTCTTTACCCCGC‑GTTGTTTTCCACACTGCTTATACCCGTCTCTTCTGAGCATCACGCAAATAAAGGCGGGTTTTTCGGTGTGTTTATTAACACCCTCgctataagcggtaaac                                                                  >  1:378732/1‑224 (MQ=255)
                                                                                                                                                        cgcgATGGTGCTGTCAGCGGAATTTATTCCGATGGCATTCTTCGGCGGTTTTAC‑TTGGGCCATTTTTCGCCCGTTCTCTCTCCCCATCCGTTCGGCAA‑TGGCCCTTTCTGTTCTGGTGGGCT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCCTTTCTGCTTAGCATCACGGAAAAAAGGGGGGGTTCTTCTGTTGGTTTTATA‑AACCTTCGACCaaaagggtaacac                                                                 >  1:488214/1‑227 (MQ=255)
                                                                                                                                                          cgtttgtgtTGTCAGTGTTTTTTTTTTCGTTGGCTTTCTTCGGCGGTTCTAC‑TGGGGCAATTTAACGCCAGTTCTCTAGCACCATCGTTGCGGCAA‑TGGCGTTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCAc                                                                <  1:46433/234‑1 (MQ=255)
                                                                                                                                                                  cTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTTTCACCATCGTTTCGGCAA‑TGGGGCTTTCTGTTCTGTGGGCAT‑TTATTTTTTCCCCTGC‑GTTATTTGCAACGCTGCTTAAACCCGTCTCTGCTGGGCATAACGAAAATAAGGGGGGTTTTTTCGGTTGGTTTAATACAACCTCTGATCAAACGGTTAACcaccaccacaa                                                        >  1:226068/1‑234 (MQ=255)
                                                                                                                                                                                ttttCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCaaaa                                          <  1:627305/237‑1 (MQ=255)
                                                                                                                                                                                  ttgcGATGGCATTCTTCGGGGGTTTGTC‑TGGGGCAATTTAGCGCCCGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTCTTTTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAGCCACTACACCAACAGCGTCGGCAAAATc                                        <  2:442955/236‑1 (MQ=255)
                                                                                                                                                                                  ttCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCCCCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATc                                        <  2:287814/239‑1 (MQ=255)
                                                                                                                                                                                       gtgtctttcttcGGCGGTTCTTC‑GGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTGCTGGTGGCAT‑TGATTCTTTCCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCgg                                   <  1:487691/233‑1 (MQ=255)
                                                                                                                                                                                       gtGGCATTTTTCGGCGGTTCCAC‑TGGGGCAATTTATCGGCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCgg                                   <  2:719605/238‑1 (MQ=255)
                                                                                                                                                                                            gtttttgttgtgTTTTAC‑TGGGGAAATTTATTCCCGGTTCCCCATCCCCTTCGTTCTGGCAA‑GGGCGCTTTCTGTTCTGGTGGCAA‑TTATTCTTACCCCTGC‑CGTATGGTCAACGCTGCTTAAACCCTGCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCGGATCca                              <  1:354344/228‑1 (MQ=255)
                                                                                                                                                                                            aTTTTTCGGGGGTTTTAC‑TGGGGCAATTTATCGCCCGTTCCTTATCCCCCGCCTTTGGGCAA‑TGGCGCTTCCTGTTCTGGTGGCAT‑TGATTCTTCCCCCTGC‑GTTATGTGCAAGGGTGCTTAAACCCGCCTCAGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCACCAGCGTCGGCAAAATCCTCGGATCca                              <  1:97348/239‑1 (MQ=255)
                                                                                                                                                                                                ttcGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCGGATCCACAgg                          <  2:226068/239‑1 (MQ=255)
                                                                                                                                                                                                ttcGGCGGGTCTTC‑TGGGGCAATTTTACGGCAGTTCTCTATCCCCATGGTTGCGGCAA‑TGGCGCTTTCCGTCCTGGTGGCAT‑TGATTGTTACCCCTGC‑GGTATGTGCAACGCGGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCGGATCCACAgg                          <  1:511057/239‑1 (MQ=255)
                                                                                                                                                                                                 tcGGGGGTTTTCC‑TGGGGCAATTTATCGCCAGTTCTTTAGCCCCCTCGTTTCGGCAA‑TGGCGTTTTCTGTTCCGGTGGCAT‑TGATTCTTTCCCTTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCGGATCCACAGGa                         <  1:230836/239‑1 (MQ=255)
                                                                                                                                                                                                      ggTTCTAC‑TGGGGCAATTTTGCGTCAGTTCTCTAGCGCCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCTT‑TGATTCGTACCACTGC‑GTTATGTGCAACGCTTCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATCCCACCTTCGATCATAGCGTTAACCACTACTCCAACAGCGTCGGCAAAATCCTCGGATCCACAGGACGat                     <  1:570780/238‑1 (MQ=255)
                                                                                                                                                                                                            tc‑TGGGGCAATTTAACGGCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCGGATCCACAGGACGATATTTACt              <  2:106994/238‑1 (MQ=255)
                                                                                                                                                                                                            tc‑TGGGGAAATTTATCCCCTGTTTTTTTTCACCGTCGTTTCGGCAA‑GGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCGGATCCACAGGACGATATTTACt              <  1:776904/238‑1 (MQ=255)
                                                                                                                                                                                                            aC‑TGGGGCAATTTACCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCGGATCCACAGGACGATATTTACt              <  2:880800/239‑1 (MQ=255)
                                                                                                                                                                                                                   gCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCGGATCCACAGGACGATATTTACTGATCTa        >  2:251762/1‑239 (MQ=255)
                                                                                                                                                                                                                       tttctccccctttttctATCCCCATCGTTTGGGCTA‑TGGCGCTTTCTGTTCTGGGGGCAT‑TGAGTCTTCCCCCGGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCGGATCCACAGGACGATATTTACTGATCTATgcg    <  1:251762/225‑1 (MQ=255)
                                                                                                                                                                                                                         ttgcGCCATTTCTCTATCGCCCGCGTTTGGGCAA‑TGGCTCTTTCTTTTTTGGGGGCAT‑TGAGTCTTACCCCTGCTGTTA‑GTGCAACTCTGGTTCACCCCGTGTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCGGATCCACAGGACGATATTTACTGATCTATGCGCt  <  1:476841/236‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                  
TGTTCGGGATGGTGCTTGCCATCGGGCTGCTCGTCGATGATGCGATAGTGGTGGTGGAGAACGTCGAGCGCGTGATGATGGAGGATAAGCTCCCGCCAAAAGAAGCGACGGAAAAATCGATGTCGCAAATTCA‑GGGCGCACTGGTGGGTATCGCGATGGTGCTGTCAGCGGTATTTATTCCGATGGCATTCTTCGGCGGTTCTAC‑TGGGGCAATTTATCGCCAGTTCTCTATCACCATCGTTTCGGCAA‑TGGCGCTTTCTGTTCTGGTGGCAT‑TGATTCTTACCCCTGC‑GTTATGTGCAACGCTGCTTAAACCCGTCTCTGCTGAGCATCACGAAAATAAGGGCGGTTTCTTCGGTTGGTTTAATACCACCTTCGATCATAGCGTTAACCACTACACCAACAGCGTCGGCAAAATCCTCGGATCCACAGGACGATATTTACTGATCTATGCGCT  >  NC_000913/3416214‑3416667

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