Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,876,107 G→A 11.8% M145I (ATG→ATA nimT → 2‑nitroimidazole exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,876,1070GA11.8% 80.5 / 3.2 34M145I (ATG→ATAnimT2‑nitroimidazole exporter
Reads supporting (aligned to +/- strand):  ref base G (18/12);  new base A (2/2);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

TCCCTTCGCCTTACTTATTATTTGGCGGTACAGCGGTCATTGGCGGTGGGATTGCATTAGGCAATGTCTTACTGCCAGGATTAATTAAACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTCCTCTGCTGGCTCTTTTTTTATGGCTGCCACAGTGGCGAAGTCAAC  >  NC_000913/1875965‑1876246
                                                                                                                                              |                                                                                                                                           
tcCCTTCGCCTTACTTATTATTTGGCGGTACAGCGGTCATTGGCGGTGGGATTGCATTAGGCAATGTCTTACTGCCAGGATTAATTAAACGCGATTTCCCTCATTCCGTCGCCCGACTTACCGGCGCATATTCCCTGACAATg                                                                                                                                             >  2:235763/1‑143 (MQ=255)
ctgtgtCGCCGTCCTTGTGATTTGGCGGTCCAGCGGTCATTGGCGGTGGGATTGCATTAGGCAATCTCTCACTGCCAGGATTAATTGAACGTGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATg                                                                                                                                             <  2:487594/138‑1 (MQ=255)
           tacttaTTATTTGGCGGTACAGCGGTCATTGGCGGTGGGATTGCATTAGGCAATGTCTTACTGCCAGGATTAATTAAACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGc                                                                                                                                  <  2:500901/143‑1 (MQ=255)
                     ttGGCGGTACAGCGGTCATTGGCGGTGGGATTGCATTAGGCAATGTCTTACTGCCAGGATTAATTAAACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGCGACTgcgg                                                                                                                        >  1:54180/1‑140 (MQ=255)
                         cGGTACAGCGGTCATTGGCGGTGGGATTGCATTAGGCAATGTCTTACTGCCAGGATTAATTAAACGCGAGTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGCTCgg                                                                                                                    >  1:93309/1‑143 (MQ=255)
                                 cGGTCAGGGGCGGTGGGATTGCATTAGGCAATGTCTTACTGCCAGGATCAATTAAACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTc                                                                                                            <  1:88243/143‑1 (MQ=255)
                                         ggCGCTGGGACTGCATTCGGCTGGGTCTTACTGCCAAGATCAATTGAACGCGATTTCCCTCATTCCGTCGTCGGACTTACCGGCGCATATTTCCTGGCAATGGGAGCTGCTGCGGCACTGGGCTCGGCTATGGTCGTGCCGCt                                                                                                    <  2:344301/143‑1 (MQ=255)
                                               gggATTGCATTAGGCAATGTCTTACTGCCAGGATTAATTAAACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCttt                                                                                              <  1:656756/143‑1 (MQ=255)
                                                       aTTAGGCAATGTCTTACTGCCAGGATTAATTAAACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGtt                                                                                      >  2:543491/1‑143 (MQ=255)
                                                        ttAGGCAATGTCTTACTGCCAGGATTAATTAAACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACACTGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGttt                                                                                     >  1:130296/1‑143 (MQ=255)
                                                              aaTGTCTTACTGCCAGGATTAATTAAACGCGGTTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTtggctg                                                                               >  1:379580/1‑143 (MQ=255)
                                                                           cAGGATTAATTAAACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATAGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTg                                                                  >  2:204665/1‑143 (MQ=255)
                                                                           cAGGATTAATTAAACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATAGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTg                                                                  >  2:204762/1‑143 (MQ=255)
                                                                            aGGATTAATTAAACGCGATGTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGc                                                                 >  2:539668/1‑143 (MQ=255)
                                                                             ggATTAATTAAACGCGATTTCCCTCATGCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCt                                                                <  1:235605/143‑1 (MQ=255)
                                                                               attaattaAACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCa                                                              >  1:12455/1‑143 (MQ=255)
                                                                                 taattaaACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATg                                                            >  1:338778/1‑143 (MQ=255)
                                                                                     taaACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGa                                                        >  1:372249/1‑143 (MQ=255)
                                                                                          gcgATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGCCTTTGAACGGCTTCGGCTGGCCAGGCGCGTTGCTCCTTCTGAtcggt                                                   >  1:37859/1‑139 (MQ=255)
                                                                                            gATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGttg                                                 >  1:538546/1‑142 (MQ=255)
                                                                                             aTTTCCCTCAGGCCGTCGCCAGACGTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGtttt                                                <  1:124912/143‑1 (MQ=255)
                                                                                              tttCCCTCCTTGCGTGGCCTGACTTGCCGGCGCATATTGCCTGACGGTGGGAGCTGCAGCGGCACTGGGATCGGCTGTGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTc                                               <  2:15444/143‑1 (MQ=255)
                                                                                              tttCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTc                                               >  2:15506/1‑143 (MQ=255)
                                                                                                 cccTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTCctc                                            >  2:304767/1‑143 (MQ=255)
                                                                                                  cgtgaTTCGGTCGCCCGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCCGCACTGGGATCGGCTATGGTCGTGCCGCGGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTCctct                                           <  1:480927/139‑1 (MQ=255)
                                                                                                               ccAGACTTACCGGCGCAGATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGACTGTCAAGTCGCGTTGCTCATGCTGATGTGTTTTCCTCTGCTGGCTCttttt                              <  1:416345/143‑1 (MQ=255)
                                                                                                                    cTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTCCTCTGCTGGCTCTTTTTTTATg                         <  1:543491/143‑1 (MQ=255)
                                                                                                                     ttACCGGCGCATATTCCCTGACAATAGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTCCTCTGCTGGCTCTTTTTTTATgg                        <  1:204665/143‑1 (MQ=255)
                                                                                                                     ttACCGGCGCATATTCCCTGACAATAGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTCCTCTCCTGGCTCTTTTTTTATgg                        <  1:204762/143‑1 (MQ=255)
                                                                                                                         cGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGCACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTCCTCTGCTGGCTCTTTTTTTATGGCTGc                    >  2:174222/1‑143 (MQ=255)
                                                                                                                         cGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTGTCCCGCTGCTGGCGCTTGTGTCAGGGCTGc                    >  2:248629/1‑143 (MQ=255)
                                                                                                                                  ttCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTCAACGGATTTGGCTGGCAAGGCGCGGTGCTCATGCTGATGTGTTTTCCTCTGCTGACGCTTATTTTATGGCTGCCACAGTgtc           >  2:361061/1‑141 (MQ=255)
                                                                                                                                     ccTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCGAGGCGCGTTGCTCATGCTGATGGGTTGTCCTCTGCTGGCTATTTCTTTATGGCTGCCGCAGTGGCGaa        >  2:584798/1‑143 (MQ=255)
                                                                                                                                     ccTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTCCTCTGCTGGCTCTTTTTTTATGGCTGCCACAGTGGCGaa        >  1:487474/1‑143 (MQ=255)
                                                                                                                                     ccTGACAATGGGAGCTGCAGCGGCACTGGGAGCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTCCTCTGCTGGCTCTTTTTTTATGGCTGCCACAGTGGCGaa        <  1:289420/143‑1 (MQ=255)
                                                                                                                                           aaTGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTCCTCTGCTGGCTCTTTTTTTATGGCTGCCACAGTGGCGAAGTcaac  <  2:338778/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
TCCCTTCGCCTTACTTATTATTTGGCGGTACAGCGGTCATTGGCGGTGGGATTGCATTAGGCAATGTCTTACTGCCAGGATTAATTAAACGCGATTTCCCTCATTCCGTCGCCAGACTTACCGGCGCATATTCCCTGACAATGGGAGCTGCAGCGGCACTGGGATCGGCTATGGTCGTGCCGCTGGCTTTGAACGGTTTTGGCTGGCAAGGCGCGTTGCTCATGCTGATGTGTTTTCCTCTGCTGGCTCTTTTTTTATGGCTGCCACAGTGGCGAAGTCAAC  >  NC_000913/1875965‑1876246

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