Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,477,207 A→C 16.4% A688A (GCA→GCC ydbD → PF10971 family putative periplasmic methylglyoxal resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,477,2070AC16.4% 165.9 / 5.0 78A688A (GCA→GCCydbDPF10971 family putative periplasmic methylglyoxal resistance protein
Reads supporting (aligned to +/- strand):  ref base A (47/16);  new base C (3/11);  total (50/28)
Fisher's exact test for biased strand distribution p-value = 3.31e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACCCGCTGTTGTTCTCAAGGAATTACAGCAAAACACGCCACGCAATGATATTTGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGTTCCCGTGCGCAGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGACGCGCGCAATACGGGAACAAGCGATGAGAAAACTTGATGAG  >  NC_000913/1476971‑1477441
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
ttACCCGCTGTTGTTCTCAAGGAATTACAGCAAAACACGCCACGCAATGATATTTGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGGTCCCGTGCGCAGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACACCGAGGCGTTAGCACAACTGGATGCACGGGCAAAACAACCGGACGGTGTATGGGGGGGGCAACAGGCGAAAAGAAATCTACCATTTTTCAAATTGGGGGCaaa                                                                                                                                                                                                                                             >  1:1314095/1‑237 (MQ=255)
    ccGCTGTTGTTCTCAAGGAATTACAGCAAAACACGCCACGCAATGATATTTGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGTTCCCGTGCGCAGGTTGACGCAGCGTTACGTGTTTATTATGCACTAGACCCCGGTGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGGATAGTGTGGCGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTacc                                                                                                                                                                                                                                         >  2:736281/1‑237 (MQ=255)
    ccGCTGTTGTTCTCAAGGAATTACAGCAAAACACGCCACGCAATGATATTTGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGTTCCCGTGCGCAGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCACTTTCGGCGCACTTaac                                                                                                                                                                                                                                         >  2:577562/1‑239 (MQ=255)
    ccGCTGTTGTTCTCAAGGAATTACAGCAAAACACGCCACGCAATGATATTTGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGTTCCCGTGCGCAGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGGATATTGTGGAGTACCTTGGCGATAAGCAATCTCACATTTTTCAAGTTCGGCGCCCTTaac                                                                                                                                                                                                                                         >  2:517704/1‑239 (MQ=255)
       ctgatcttCTCAAGGAATTACAGCAAAACACGCCACGCAATGTTATTTGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGTTCACGTGCGCAGGGTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTAATGGCAAAAAAAACAGGACGTATATGGTGGAATACACTGGCGAAAAGCAATCTCACTTTTTTCAAGTTCGGCGCACTTaacaac                                                                                                                                                                                                                                     >  1:1082401/7‑239 (MQ=255)
             ttCTCAAGGAATTACAGAAAAACCCGCCACGCAATGATATTGGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGTTCCCGTGGGCAGGTGGACGAAGCTTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTCCGGCGCACTTAACAACCGCcac                                                                                                                                                                                                                              <  2:586651/239‑1 (MQ=255)
                gcAAGTAATTAACGCAAAACACCCCCCGCAATGATTGTTGGGCGAAGTTTTTGCTGAGGCAGAAACACGTTTCCCGGGCGGGGGTTGCAGAAGCGTTACTGTTTTATGAAGCACAAGACCCCGATGGTTGGGCACAACGGGATGTCCTGGAAAAAAAACCGGTTCGTATATGGGGGAGTACTCTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCCTTTAACAACCGCcacaca                                                                                                                                                                                                                           <  1:319793/238‑1 (MQ=255)
                    ggAATTACAGCAAAACACGCCACGCAATGATATTTGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGTTCCCGTGCGCAGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTCCTGGCAAAACAACCGGACCGTATATGGTGGGGTACACTGCCGCCAAGCAATCTCACATTTTTCAAGTTTGGGGCACTTAACAACCGCCACACCCCg                                                                                                                                                                                                                        >  2:448492/1‑238 (MQ=255)
                            aGCAAAACACGCCACGCAATGATATTTGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGTTCCCGTGCGCAGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAGAAGCAACCTCACATTTTTCCAGTTCCGGCCCCTTAAAAACCCCCCacacacaccagccaaa                                                                                                                                                                                                                >  1:947211/1‑229 (MQ=255)
                                     cgccaggcAATGATTTTTGGGCGAAGTTTTTTCTGAGGAAGGAAAACAGTTCCCGGGCGCAGGTTGACGAAGCGTTACGTTTTTATTATGCACTAGACCCCGATGCTTTAGCACACCTGGATGTACTGGCAAACCAACCGGATCTTATATGGTGGAGTACACTGGCGAAAAGAAATCTCACATTTTTCAAGTTCGGCGCCCTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCa                                                                                                                                                                                                      <  2:527004/239‑1 (MQ=255)
                                      gccacgcAATGATATTTGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGTTCCCGTGCGCCGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGGATATGGTGGAGTAACCTGGCGAAACGCCATCCCACCTTTTTCAAGTTTGGCGCCATTAACAACCCCCACACACCCCCGCCAGCACCggcggccg                                                                                                                                                                                                     >  1:1291200/1‑237 (MQ=255)
                                      gccacgcAATGATATTTGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGTTCCCGTGCGCAGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGCAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACAGCCCCACACCGCCTGCAGACCTGGCGGCAg                                                                                                                                                                                                     >  2:897133/1‑239 (MQ=255)
                                        cacgcAATGATATTTGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGTTCCCGTGCGCAGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGaa                                                                                                                                                                                                   >  2:1106626/1‑239 (MQ=255)
                                                 atatTGGGGGAAAGTTTTTTCTGAGGCAGGAAAACGGTTCCCGGGCGCAGGTTGACGAAGCGTTACGTGTTTATTTGGCACTAGCCCCCGATGCGTTAGCACAACGGGATGTACTGGAAAAACACCCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATccc                                                                                                                                                                                          <  2:189993/239‑1 (MQ=255)
                                                 atatTGGGGCGAATTTTTTTCTGAGGCAGGAAAACGGTTCCCGTGCGCAGGTTGAGGAAGCGTTACGTGTTTATTATGCCCTAGCCCCCGATGCTTTAGCACAACTGGATGTACTGGAAAAACACCCGGATCGTATATGGTGGAGTACACTGGCAAAAAGCAATCTCACATTTTCCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATccc                                                                                                                                                                                          <  2:170676/239‑1 (MQ=255)
                                                  tgtttGGGGAAATTTTTTTTGGGGGCAGAAAACCAGTTCCCGGGCGAAGTTTGACGAAGGGTTACGTTTTTTGTACGCACTAGACCCCGAGTCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGACCGTATAGGGGGGAGTACACTGGCGAAAAGAAATCTCACATTTTTAAAGTTGGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATccc                                                                                                                                                                                          <  1:998341/236‑1 (MQ=255)
                                                      ttgggtaggTTTTTTTTGAGCAAGGAAAACATGCCCTGTGACGGGTTGGACGAAGTGTATCGTTTTTATGAGCCCTTAGCCCCCGATGGTTTAGTCAAACGGGTTGTACTGGAAAAACAACCGTGTCGTTCTGGGGGGAGTCCACGGGCGAAAGCCAATCCCAGATTTTCAAAGTTCGCCGCCCTTACAACCCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAgtg                                                                                                                                                                                     <  1:897133/231‑1 (MQ=255)
                                                                         ggCAGGAAAACAGTTCCCGTGCGCAGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGAACACTGGCGAAAAGCAATCTAACATTTTTCAAGTTCGGCGCACTTAACAACCGCCCCACACCGCCTGCAGTACTGGCGGCAGAAATTAATCACCAAGGGGGGAATGTGGGGACGaaa                                                                                                                                                                  >  1:34102/1‑238 (MQ=255)
                                                                                    aGTTCCCGTGCGCAGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGttt                                                                                                                                                       >  2:943483/1‑239 (MQ=255)
                                                                                      ttCCCGTGCGCAGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCCCACACCGCCTGCAGTCCTGGCGGCATAAATTGATCCCGAGTGGTGGATTGTGCGGATGAATATTCCCCGtttt                                                                                                                                                      >  2:847187/1‑238 (MQ=255)
                                                                                              cgcgGGTTGACGAGGCTTTACGTTTTTTTTAGGCACGAGACCCCGATGCGTTAGCACACCTGGATGTACTGGCAAAACAACCGGATCGTATATGGGGGAGTACACTGGCGAAAAGCAATCTCACATTTTCCAAGTTCGGCCCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGa                                                                                                                                              <  2:686262/238‑1 (MQ=255)
                                                                                                           atgttgttgttttttAATATCCAATGGCCCCCGTTGGGTTGACACAATGGGTTGTACGTGCAAAACAACGGGTTCATATATGGGTGCATCCTAGGGCGAAAGCCAATCTCACATTTTTCAGGTCCGGAGCCATTATCAACCGCCACCCACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCg                                                                                                                                 <  1:1082831/227‑1 (MQ=255)
                                                                                                            aGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGCTTCTGGCGATGGATAATCCCCGTTTTCCAGTTGATGTATTACAGGCgg                                                                                                                                >  2:187829/1‑238 (MQ=255)
                                                                                                              cGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGCAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCTCCACACACCGCCTCCAGTCATGGCGGCAGAAATTAATCCCGGACGGTGGCATCTGGGGTTGAAAAATCCCCGGTTTcaggtgtgtgggtaaaaggggggg                                                                                                                              >  1:689989/1‑216 (MQ=255)
                                                                                                                        ttattaTGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGCATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACAAACCGCCTGCCGTACTGCCGGAGGAAATTTACCCCAAGGGGGGGATTTTGGCCATAAAAAATCCCCGTTTTCAACTTTAGGGGTTAAAGGctgggtgtaaaaggg                                                                                                                    >  1:196538/1‑225 (MQ=255)
                                                                                                                         ttttatGCTCTAGCCCCCGCTGGTTGAGCCAGACTGGTGTACCGGACAAACCAACCGTGTCGTATAGGGGGGGTTACCCTGGCAAAAAGCAATCTCACGTTTTTCAAGTGCGGCGCATTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGa                                                                                                                   <  1:70386/237‑1 (MQ=255)
                                                                                                                          attatGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAAAAATCCCCGTTTTCCAGTTGCAGGATTAAAGGCGAGGCTGAAACGGGa                                                                                                                   >  2:70386/1‑238 (MQ=255)
                                                                                                                          attatGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGAAGAATAATCCCCGGTTTCCCGGTGATTGATTAAAGGCGAGGCTGAAACGTGa                                                                                                                   >  2:810732/1‑238 (MQ=255)
                                                                                                                             atGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCCGTACTGGCGGCAGAAATTGATCCCGGGTGGGGGCTTGTGGCGCTGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGaaggggaaaagagat                                                                                                                  >  1:579462/1‑222 (MQ=255)
                                                                                                                                  gtgaccccgtgtggttttgCAAAACTGGTGTACTTGAAAAAACATCCGGTTTGTATTGGGTGGCTTACCCGGGCAAAACCCAACCCACCTTTTTCCAAGTTGGGCGCCTTTACCAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCt                                                                                                          <  1:1106626/221‑1 (MQ=255)
                                                                                                                                  cctcatatttttataaggaacagaCTGGATCTACTGAAAAAACAACCGGTACGTTTGTGGGGGATTCCACTGGCAAAAAACAATCTCACTTTTTCCAAGGCCGGCCCATTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCt                                                                                                          <  1:1256906/216‑1 (MQ=255)
                                                                                                                                  cTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGGGGTGGATTGTGGCGATGAAGAATCCCCGTTTTCCCGTTGTTGGATTAAAGGCGAGGCTGAAACGTGATCCttt                                                                                                             >  2:91295/1‑236 (MQ=255)
                                                                                                                                     gACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCCGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGCGGAGGCTGAAACCGGATCCTTTGCTGGc                                                                                                       >  2:592681/1‑239 (MQ=255)
                                                                                                                                      aCCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGGTGAAAAATCCCCGTTTTCCCGTTGACTGATTAAAGGGCAGGGTGaaaaggagacctttggggggg                                                                                                      >  1:861665/1‑221 (MQ=255)
                                                                                                                                       ccccgtggtttgaacacaGCTGGTTGTAGGGGAAAAACAACCGGAGTGTATGGGGGGGGTTCCCCGGCCGAAAGGCACTCCCACTTTGTTCAAGTCGGGCGATATTCACACCCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCg                                                                                                      <  1:943483/224‑1 (MQ=255)
                                                                                                                                        cccGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGATTAATCCCCGTTTTCCAGTTGACGTATTAAAGGCAAGGCTGAAACGGGACCCTTTCCTGGGGtt                                                                                                    >  2:1291584/1‑239 (MQ=255)
                                                                                                                                          cGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCGGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGGGATCCTTTGCTGGCGttta                                                                                                  >  2:191806/1‑237 (MQ=255)
                                                                                                                                          cGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGGGGTGGATTGTGGCGGTGAATAATCCCCGTTTTCCAAGTGATGTATTAAAGGGGAGGCTGAAAAGAGAGAATTTGCCGGCGGTaa                                                                                                  >  1:454555/1‑238 (MQ=255)
                                                                                                                                                 ttAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGCGGGGGGGATTGTGGCGCAGAATAATCCCCGGTTTCCAGTTGATTGATTAAAGTCGAGGCGGAAAGATGACCCTTTGCTggaggtctagctggt                                                                                           >  1:888775/1‑226 (MQ=255)
                                                                                                                                                  tAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGGTTTAAAGGCGAGGCTTGACAGGGATCCTTTCCTGGCGTCAGACCTTGtt                                                                                          >  2:516557/1‑239 (MQ=255)
                                                                                                                                                   aGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGGTTTCCCGTTGTTGTATTAAAGGCCAGGCTGCAACGTGGTCCTTTTCTGGCGTTGAAGCTTGTTa                                                                                         >  2:113149/1‑239 (MQ=255)
                                                                                                                                                     cacaACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAGGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCCCCTGCAGTACTGGCGGCCGAAATTGATCCCGAGTGGTGGCTTTCGGCGAGGAAAAATCCCAGTTTTCCACGTGCTGGATTAAAGGCGAGGCTGAAACGGGATCCTTTTTCGGCGTTGAAGCTGTTTaa                                                                                        >  1:1163493/1‑238 (MQ=255)
                                                                                                                                                            ttggtggaacgggaaaaaaaACCGGATCTTAGTGGGGGGATTACACTGGCAAAAACCAATCTCACATTTTTAAAGTTCGGCGCCCTAAACAGCCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGa                                                                                  <  2:541750/225‑1 (MQ=255)
                                                                                                                                                             gagttACGGGAAAAACAACCGGATGGTATATGGTGGAGTACACTGGCAAAAAGCAATCTCACATTTTCCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAAc                                                                                <  2:947211/234‑1 (MQ=255)
                                                                                                                                                                  aCTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCCGCAATTGATCCCGAGTGGGGGATTGTTGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGGGGGCGAAACGTCATCCTTTTCTGGCGTTAGAGCTTTTTAACACCGAACCGGAtt                                                                         >  1:49531/1‑239 (MQ=255)
                                                                                                                                                                  aCTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGGTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGGGATCCTTTCCTGGGCTTAAAGCTTGTTAACCCTGAACTGGAtt                                                                         >  2:587867/1‑239 (MQ=255)
                                                                                                                                                                    tGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCCGAAATTGATCCCGAGTGGGGGATTGTGGCGCGGAATAATCCCCGTCTTCCAGGTGATGGTTTTAAGGCGTGGTCGAAACGGGTACCTTTGCTGGCGGTACAGTTTTTTAACCCGGAAACGGAttaa                                                                       >  1:316740/1‑237 (MQ=255)
                                                                                                                                                                    tGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGGATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGGGATCCTTTGCTGGCGtttaag                                                                                                >  2:1256906/1‑211 (MQ=255)
                                                                                                                                                                    tGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTTAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAAAGCTTGTTAACCCTGAACTGGATTTa                                                                       >  2:612090/1‑239 (MQ=255)
                                                                                                                                                                     ggCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGGATAATACCCGTTTTCCAGTTGATGTATTAAAGGCGAGGATGAAACGGGATCCTTTGCTGGCGTTAAGGCTTGGTAACCCTGAACCGGAGTTAg                                                                      >  2:857135/1‑239 (MQ=255)
                                                                                                                                                                     ggCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCCGTTGATGGATTAAAGGCGTGGCTGGAAACTGATCCTTTTCTGGGGTTAGAGCTTGTTAAACCTGAAATGGATTTaa                                                                      >  1:15368/1‑238 (MQ=255)
                                                                                                                                                                        cAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGGGGGTTGGGGCGATTAAAAATCCCCGATTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGGGTTCCTTTGCTGGCGTTTGAGTTGTTTAAACCTGAAATGGGTTTatga                                                                    >  1:307606/1‑236 (MQ=255)
                                                                                                                                                                                  cGGAGCGTATATGGTGGAGTACCCGGGCGAAAACCACTCCCCCATTTTCAAAGTGCGGCGCCCTTACCAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTg                                                           <  2:888775/239‑1 (MQ=255)
                                                                                                                                                                                       cGTATGTGGGTGACTACACTGGCGAAACGCAATCTCAGTTTTTATAAGTGCGGCGCCCTTACCAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCt                                                      <  2:1082401/239‑1 (MQ=255)
                                                                                                                                                                                        tttattggggggtgaccatGGCAAAACGCAATCCCACTTTTTGAAAGTTCGGCGCAGTTAAAAACCGCCACCCACCGCCTGGAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCt                                                      <  2:1113358/220‑1 (MQ=255)
                                                                                                                                                                                        ggatatGGTGGAGTCCCGGGGCGAAAAGCAATCCCACTTTTTTCAAGTCCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCtt                                                     <  2:49531/237‑1 (MQ=255)
                                                                                                                                                                                         tataggtgggagacaccgggCAAAAAGAAATCCACCATTTTTAAAGTCCGGCCCCCTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTa                                                    <  2:15368/221‑1 (MQ=255)
                                                                                                                                                                                               gtggAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCTGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCACTTTTCCCGTTGATGTATTAAAGGCGAGGCTGAAAAGTGATCCTTTTCTGGGGTTAGAGCTTGGTGAACCTGGACTGGGTTTTAGTccccaaccgggagctaaaggta                                               >  1:766747/1‑219 (MQ=255)
                                                                                                                                                                                               gtggAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTa                                              >  2:802820/1‑239 (MQ=255)
                                                                                                                                                                                                tggAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTaa                                             >  2:58371/1‑239 (MQ=255)
                                                                                                                                                                                                   aGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCg                                                                                                                                                                        >  1:534792/1‑114 (MQ=255)
                                                                                                                                                                                                   aGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGAc                                          >  2:732548/1‑239 (MQ=255)
                                                                                                                                                                                                   aGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGGGGTGGAGTGTGGCGGTGAATAATCCCCGTTTTCCAGTTGATGTTTTTAAGGCGGGGTTGAAAAGGGTTCCGTTGGTGGGGTTTAAATTTTTTACCCCTGAACGGTGCTTAGTCGCCCAGCTGGTGCATaaagcaaaaaa                                          >  1:67567/1‑230 (MQ=255)
                                                                                                                                                                                                   aGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCCAGTGGTGGATTGTGTCGGTGAATAATCCCCGTTTTTCCGTTGATGTATTAAAGGCAGGGTGGGAACGGGATCCTTTTCTGGCGTTATAGCTTTTTAACCCCGAACTGGGTTTAATCCCCCACACTGCGCTTAAAGGttaaaa                                          >  1:1287289/1‑234 (MQ=255)
                                                                                                                                                                                                   aGTACACTGGCAAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCg                                                                                                                                                                        <  2:534792/114‑1 (MQ=255)
                                                                                                                                                                                                       taccgggcaaaaGCCATTCTCACATTTTTCAATTTCGGCGCCCTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGAcgcg                                       <  2:1212133/231‑1 (MQ=255)
                                                                                                                                                                                                          ttggcaaaaGCAAATCTCACTTTTTTCAGTTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGACgcgcgc                                    <  2:316740/234‑1 (MQ=255)
                                                                                                                                                                                                          tGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGACGCGCGCa                                   >  2:147866/1‑239 (MQ=255)
                                                                                                                                                                                                          tGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGAGGAGCTGAATAATCCCCGCTTTCCAGTTGGTGTATTAAAGGCGGTGTGGAAACGTGATCCTTTTCTTGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTGAAAcgcgcgac                                   >  1:1223935/1‑237 (MQ=255)
                                                                                                                                                                                                          tGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGAACCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAAAGTGCTCCTTTTCTGGCGTTTGTGCGTTTTAACCCTGAACTGGATGTAATCCCCCAGCGGGGGCTTAAAggagaactcggccca                                   >  1:280508/1‑226 (MQ=255)
                                                                                                                                                                                                           cctgagaaaaaccactcctattttttAAAGGTTCGGCGCAATTAAACACCGCCCACCGCCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGAcgcgcg                                     <  1:592681/218‑1 (MQ=255)
                                                                                                                                                                                                            gCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGGTTTCCAGTTGATGTATTAAAGGCAAGGCTGAAACGTGTTCCTTTGCTGGCGTTAGAGTTTGTTAACCCTTAACCGGATTTAGTCCCCCCGCTGGCCCTTAACGTTAAGaaccccccaaa                                 >  1:1037442/1‑229 (MQ=255)
                                                                                                                                                                                                               aaaaGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTCAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGACGCCCGCAATACg                              >  2:443986/1‑239 (MQ=255)
                                                                                                                                                                                                                gaaaacgacttcaccttttttCAAGTGTGGCACCATTCACCACCCCCCACCCACCGCCTGGAGTACTGGCGGCAGAAATTGATCCCGTGTGGTGGATGGTGGCGATGAATAATCCCCGTTTTCCGGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGACGCGCGCAATa                                <  1:528075/223‑1 (MQ=255)
                                                                                                                                                                                                                  aGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGCTTGTGGCGATGAATAATCCCCGTTTTCCAGTTGGTGTATTAAAGGCGAGGCTTGAACTGGCTCCTTTTTTGGGGTTTAAGCTTGTTAACCCTGAAACTGATTTTATCCGCCAGCTGGCGCTTAACgcgaaaacccgccccaaaaggga                           >  1:206128/1‑217 (MQ=255)
                                                                                                                                                                                                                       aTCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGGGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGACGCGCGCAATAAGGGAGCaa                       >  2:1318017/1‑239 (MQ=255)
                                                                                                                                                                                                                                tttttAAAGTGCGGGCCACTAAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGACGCGCGCAATACGGGAACAAGCGATgag               <  2:1287289/239‑1 (MQ=255)
                                                                                                                                                                                                                                            tCGGCGCACTTAACAACCGCCACCACCCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGACGCGCGCAATACGGGAACAAGCGATGAGAAAACTtgat     <  2:454555/239‑1 (MQ=255)
                                                                                                                                                                                                                                                cccACGTAAAACCCGCCACACACCGGCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGACGCGCGCAATACGGGAACAAGCGATGAGAAAACTTGATGAg  <  2:529943/236‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
TTACCCGCTGTTGTTCTCAAGGAATTACAGCAAAACACGCCACGCAATGATATTTGGGCGAAGTTTTTTCTGAGGCAGGAAAACAGTTCCCGTGCGCAGGTTGACGAAGCGTTACGTGTTTATTATGCACTAGACCCCGATGCGTTAGCACAACTGGATGTACTGGCAAAACAACCGGATCGTATATGGTGGAGTACACTGGCGAAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGCGGCAGAAATTGATCCCGAGTGGTGGATTGTGGCGATGAATAATCCCCGTTTTCCAGTTGATGTATTAAAGGCGAGGCTGAAACGTGATCCTTTGCTGGCGTTAGAGCTTGTTAACCCTGAACTGGATTTAGTCCGCCAGCTGGCGCTTAACGGTAAGACGCGCGCAATACGGGAACAAGCGATGAGAAAACTTGATGAG  >  NC_000913/1476971‑1477441

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