Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,171,788 G→T intergenic (+32/‑269) rrfB → / → murB 5S ribosomal RNA/UDP‑N‑acetylenolpyruvoylglucosamine reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,171,7880GT92.9% 37.3 / ‑0.8 29intergenic (+32/‑269)rrfB/murB5S ribosomal RNA/UDP‑N‑acetylenolpyruvoylglucosamine reductase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (11/15);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

TTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCAT  >  NC_000913/4171543‑4172005
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                         
ttACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCt                                                                                                                                                                                                                      >  1:277574/1‑251 (MQ=37)
 tACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCTGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCAGGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGcc                                                                                                                                                                                                                       <  1:523448/249‑1 (MQ=255)
     aCGCCGAAGCTGTTTGGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAAGCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGTGCCCACCTGCCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGTGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGt                                                                                                                                                                                                                 <  2:26391/251‑1 (MQ=11)
      cGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTGGGCCTTTCGtt                                                                                                                                                                                                                <  1:193163/251‑1 (MQ=37)
                 ttgtGGCTGATGAGAGAATAGTTGCAGCCTGATACAGATTAAATCAGAACGCAGTAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCTGTGGTCCCACCTGACCCCAGTCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCTAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGTCCTTTCGTTTTATCtgttg                                                                                                                                                                                                      <  2:502060/247‑1 (MQ=11)
                                gaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGct                                                                                                                                                                                      <  1:477761/251‑1 (MQ=18)
                                         cAGCCTGATAATGATTAAATCAGAACGCAGAAGCGGACTGATAAAACAGAATATGCCTGGCGGCAGCAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCCGAAGTGAAACGCCGTAGAGCCGGTGGTAGTGTGGGGCCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATCAAACGAAAGGCTCAG‑ACGAAAGACTCGACCTTTCGTTTTATCTTTTGGTTTTCGGTGAACGCTCTCCTGAgg                                                                                                                                                                             >  2:485246/1‑250 (MQ=11)
                                              tGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCTAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGGAGGAc                                                                                                                                                                        >  1:174303/1‑251 (MQ=18)
                                                aTACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGGTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACaa                                                                                                                                                                      <  1:181080/251‑1 (MQ=18)
                                                    agctttaaTCAGAGCGCAGAAGCGGTCTGATAAAACAGAGTTGGCCAGGCAGAAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGACCTCAGAAGGGGAACGACGTAGCGACGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGCCATCCAATAAAACGACAGGATCAG‑TCGCCAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGGGTAGGAAAAATc                                                                                                                                                                   >  2:41389/7‑250 (MQ=11)
                                                        tAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                                                                                                                                                               <  1:504815/250‑1 (MQ=18)
                                                              aGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCAGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCg                                                                                                                                                        >  1:499986/1‑251 (MQ=18)
                                                                    aaGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCGCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTg                                                                                                                                                  <  2:355167/250‑1 (MQ=18)
                                                                       aaGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGAGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                                                                                                               <  1:523483/251‑1 (MQ=18)
                                                                           ggTCTGATAAAACAGAATTTGCCTGGCGGCAGGAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTAGCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTg                                                                                                                                           <  2:345617/251‑1 (MQ=18)
                                                                                    aaaCAGAAGTGGCCCGGCGGCAGTGGCGCGGTGGTCCCACCTGACCCCAGGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAGCGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAAc                                                                                                                                  <  2:19127/251‑1 (MQ=11)
                                                                                                   ggcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTgg                                                                                                                    <  1:137221/250‑1 (MQ=18)
                                                                                                     cggcAGTAGCGCGGTGGTCCCACCTGACACCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGATGGTGGCgg                                                                                                                 >  1:160439/1‑251 (MQ=18)
                                                                                                                                       ccGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAgg                                                                                                                                                                          >  1:141158/1‑160 (MQ=25)
                                                                                                                                       ccGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAgg                                                                                                                                                                          <  2:141158/160‑1 (MQ=25)
                                                                                                                                         gAACTCAGAAGTGAAACTCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGAGGATTTGAACGTTGCGAAGCAACGGCCCAGGGGGTGGCGGGCAGGCCGCCCGCCCTAAACTGcccggactcacatt                                                                             >  2:82356/1‑239 (MQ=11)
                                                                                                                                                        ccgcggTAGCGACGTTGGTGGTGTGGGGTCTCCCCATGGGAGAGGAGGAAACTGCCAGGAATAAAATAAAACGAAAGGCTCAG‑TCGAAAGGTTTGGCTTTTGTTTTTATCTTTTTTTTGGCGGTGACCGCTCTCTTGATGAGAAAAAACCCGCCGGGGGCGGATTTGAATGTTGCGAAGCAACGGCCGGGGGGGTGGCGGGCAGGAGGGCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGTCCATc                                                              <  1:446739/246‑1 (MQ=2)
                                                                                                                                                            gttgcccctgATGGGAGTGGGGGGTTTCCCTAGGCGAGAGGGGGGAACTGCCAGGCATAAAAAAAAACGAACGGC‑CAGGTTAAAAGACTTGTCCTCTTGTTTTCTTTGTTGTTTGGGTGTAAACGCCCTCGTGCTTAGGAAAACTCGGCTGGGTGCGGATTTGCCCGTGGCTAACCAACGGCCGGGAGGGGGGGGGGAAGCACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGa                                                          <  1:275/242‑1 (MQ=2)
                                                                                                                                                            cGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCACCCAAAACTGCCAGGGCTCAAATTAAGGAGAAGGCCATccgac                                                          >  2:169746/1‑248 (MQ=11)
                                                                                                                                                                gggcCGAGGGTAGTGGGGGTTCTCCCCAGGCGAGAGTAGGGACCGGCCAGGAATAAAAAAAAACGAAAGGCTCAG‑TCGAAAGTCTTGGCTTTTCTTTTTTTCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAACTCCGCCGGGGGCGTATTTGAGCGTTGCGAAGAAAGGGCCGGGGGGGTGGGGGGAAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGa                                                      <  1:482883/249‑1 (MQ=2)
                                                                                                                                                                         gTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCtttt                                             >  1:485631/1‑251 (MQ=18)
                                                                                                                                                                                     ctcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAg                                                                          >  1:241412/1‑210 (MQ=18)
                                                                                                                                                                                                                    aacagaTAAAACGAAAGGCTCAG‑TCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                                                                                                            <  1:275056/140‑1 (MQ=12)
                                                                                                                                                                                                                    aTCAAATAAAACGAAAGGCTCAG‑TCGCAAGACTTGGCCTTTCGATTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACACCCGCCATACACTGCCAGGCCTCCAATTAAGCCGAAGGCCCCCCTGACGGATGGCCTTTTTGCGATTATACCAACTCTTCCTGTAGTCATATCTCCAAGCCAt  >  2:56731/1‑251 (MQ=14)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                         
TTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAG‑TCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCAT  >  NC_000913/4171543‑4172005

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