Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 78,830 T→A 20.8% intergenic (+31/+18) setA → / ← leuD broad specificity sugar efflux system/3‑isopropylmalate dehydratase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091378,8300TA20.8% 53.7 / 8.8 24intergenic (+31/+18)setA/leuDbroad specificity sugar efflux system/3‑isopropylmalate dehydratase small subunit
Reads supporting (aligned to +/- strand):  ref base T (9/10);  new base A (3/2);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAG  >  NC_000913/78686‑78971
                                                                                                                                                |                                                                                                                                             
ggcGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTg                                                                                                                                            <  3:5041/148‑1 (MQ=255)
                   ttGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAAt                                                                                                                         >  2:170521/1‑148 (MQ=255)
                        caAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAt                                                                                                                                               >  4:755459/1‑121 (MQ=255)
                        caAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAt                                                                                                                                               <  3:755459/121‑1 (MQ=255)
                                        ttGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTg                                                                                                                                            <  1:300387/108‑1 (MQ=255)
                                        ttGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTg                                                                                                                                            >  2:300387/1‑108 (MQ=255)
                                                 aCTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGAc                                                                                                                                          >  1:922467/1‑101 (MQ=255)
                                                 aCTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGAc                                                                                                                                          <  2:922467/101‑1 (MQ=255)
                                                          ttGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATgg                                                                                  <  1:170521/148‑1 (MQ=255)
                                                                             cGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAtt                                                                                                                                              <  1:272988/69‑1 (MQ=255)
                                                                             cGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAtt                                                                                                                                              >  2:272988/1‑69 (MQ=255)
                                                                                attattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCa                                                                                                              <  1:971239/98‑1 (MQ=255)
                                                                                attattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCa                                                                                                              >  2:971239/1‑98 (MQ=255)
                                                                                  tattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCgg                                                                                                                                       <  2:618390/71‑1 (MQ=255)
                                                                                  tattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCgg                                                                                                                                       >  1:618390/1‑71 (MQ=255)
                                                                                              cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                   >  3:695662/1‑79 (MQ=255)
                                                                                              cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                   <  4:884902/79‑1 (MQ=255)
                                                                                              cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                   <  4:695662/79‑1 (MQ=255)
                                                                                              cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                   >  3:884902/1‑79 (MQ=255)
                                                                                                                 atgaCGTGGACGATAGCGGAAAGCCCGGTCAAATGACCg                                                                                                                                        >  3:1128162/1‑39 (MQ=37)
                                                                                                                 atgaCGTGGACGATAGCGGAAAGCCCGGTCAAATGACCg                                                                                                                                        <  4:1128162/39‑1 (MQ=37)
                                                                                                                                          cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCt                                                                                                  <  1:792201/52‑1 (MQ=255)
                                                                                                                                          cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCt                                                                                                  >  2:792201/1‑52 (MQ=255)
                                                                                                                                          cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAGGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGc                                                           <  3:854996/91‑1 (MQ=255)
                                                                                                                                          cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAGGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGc                                                           >  4:854996/1‑91 (MQ=255)
                                                                                                                                          cGGTCAAATGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAg  >  3:527078/1‑148 (MQ=255)
                                                                                                                                            gTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATgg       <  1:38876/141‑1 (MQ=255)
                                                                                                                                            gTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATgg       >  2:38876/1‑141 (MQ=255)
                                                                                                                                                |                                                                                                                                             
GGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAG  >  NC_000913/78686‑78971

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: