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

TTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCAT  >  NC_000913/4171543‑4172005
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                         
ttACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCt                                                                                                                                                                                                                      >  1:277574/1‑251 (MQ=37)
 tACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCTGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCAGGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGcc                                                                                                                                                                                                                       <  1:523448/249‑1 (MQ=255)
     aCGCCGAAGCTGTTTGGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAAGCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGTGCCCACCTGCCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGTGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGt                                                                                                                                                                                                                 <  2:26391/251‑1 (MQ=11)
      cGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGtt                                                                                                                                                                                                                <  1:193163/251‑1 (MQ=37)
                 ttgtGGCTGATGAGAGAATAGTTGCAGCCTGATACAGATTAAATCAGAACGCAGTAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCTGTGGTCCCACCTGACCCCAGTCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCTAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGTCCTTTCGTTTTATCtgttg                                                                                                                                                                                                      <  2:502060/247‑1 (MQ=11)
                                gaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGct                                                                                                                                                                                      <  1:477761/251‑1 (MQ=18)
                                         cAGCCTGATAATGATTAAATCAGAACGCAGAAGCGGACTGATAAAACAGAATATGCCTGGCGGCAGCAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCCGAAGTGAAACGCCGTAGAGCCGGTGGTAGTGTGGGGCCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATCAAACGAAAGGCTCAGACGAAAGACTCGACCTTTCGTTTTATCTTTTGGTTTTCGGTGAACGCTCTCCTGAgg                                                                                                                                                                             >  2:485246/1‑250 (MQ=11)
                                              tGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCTAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGGAGGAc                                                                                                                                                                        >  1:174303/1‑251 (MQ=18)
                                                aTACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGGTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACaa                                                                                                                                                                      <  1:181080/251‑1 (MQ=18)
                                                    agctttaaTCAGAGCGCAGAAGCGGTCTGATAAAACAGAGTTGGCCAGGCAGAAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGACCTCAGAAGGGGAACGACGTAGCGACGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGCCATCCAATAAAACGACAGGATCAGTCGCCAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGGGTAGGAAAAATc                                                                                                                                                                   >  2:41389/7‑250 (MQ=11)
                                                        tAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                                                                                                                                                               <  1:504815/250‑1 (MQ=18)
                                                              aGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCAGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCg                                                                                                                                                        >  1:499986/1‑251 (MQ=18)
                                                                    aaGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCGCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTg                                                                                                                                                  <  2:355167/250‑1 (MQ=18)
                                                                       aaGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGAGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                                                                                                               <  1:523483/251‑1 (MQ=18)
                                                                           ggTCTGATAAAACAGAATTTGCCTGGCGGCAGGAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTAGCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTg                                                                                                                                           <  2:345617/251‑1 (MQ=18)
                                                                                    aaaCAGAAGTGGCCCGGCGGCAGTGGCGCGGTGGTCCCACCTGACCCCAGGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAGCGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAAc                                                                                                                                  <  2:19127/251‑1 (MQ=11)
                                                                                                   ggcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTgg                                                                                                                    <  1:137221/250‑1 (MQ=18)
                                                                                                     cggcAGTAGCGCGGTGGTCCCACCTGACACCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGATGGTGGCgg                                                                                                                 >  1:160439/1‑251 (MQ=18)
                                                                                                                                       ccGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAgg                                                                                                                                                                          >  1:141158/1‑160 (MQ=25)
                                                                                                                                       ccGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAgg                                                                                                                                                                          <  2:141158/160‑1 (MQ=25)
                                                                                                                                         gAACTCAGAAGTGAAACTCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGAGGATTTGAACGTTGCGAAGCAACGGCCCAGGGGGTGGCGGGCAGGCCGCCCGCCCTAAACTGcccggactcacatt                                                                             >  2:82356/1‑239 (MQ=11)
                                                                                                                                                        ccgcggTAGCGACGTTGGTGGTGTGGGGTCTCCCCATGGGAGAGGAGGAAACTGCCAGGAATAAAATAAAACGAAAGGCTCAGTCGAAAGGTTTGGCTTTTGTTTTTATCTTTTTTTTGGCGGTGACCGCTCTCTTGATGAGAAAAAACCCGCCGGGGGCGGATTTGAATGTTGCGAAGCAACGGCCGGGGGGGTGGCGGGCAGGAGGGCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGTCCATc                                                              <  1:446739/246‑1 (MQ=2)
                                                                                                                                                            gttgcccctgATGGGAGTGGGGGGTTTCCCTAGGCGAGAGGGGGGAACTGCCAGGCATAAAAAAAAACGAACGGCCAGGTTAAAAGACTTGTCCTCTTGTTTTCTTTGTTGTTTGGGTGTAAACGCCCTCGTGCTTAGGAAAACTCGGCTGGGTGCGGATTTGCCCGTGGCTAACCAACGGCCGGGAGGGGGGGGGGAAGCACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGa                                                          <  1:275/242‑1 (MQ=2)
                                                                                                                                                            cGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCACCCAAAACTGCCAGGGCTCAAATTAAGGAGAAGGCCATccgac                                                          >  2:169746/1‑248 (MQ=11)
                                                                                                                                                                gggcCGAGGGTAGTGGGGGTTCTCCCCAGGCGAGAGTAGGGACCGGCCAGGAATAAAAAAAAACGAAAGGCTCAGTCGAAAGTCTTGGCTTTTCTTTTTTTCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAACTCCGCCGGGGGCGTATTTGAGCGTTGCGAAGAAAGGGCCGGGGGGGTGGGGGGAAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGa                                                      <  1:482883/249‑1 (MQ=2)
                                                                                                                                                                         gTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCtttt                                             >  1:485631/1‑251 (MQ=18)
                                                                                                                                                                                     ctcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAg                                                                          >  1:241412/1‑210 (MQ=18)
                                                                                                                                                                                                                    aacagaTAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                                                                                                            <  1:275056/140‑1 (MQ=12)
                                                                                                                                                                                                                    aTCAAATAAAACGAAAGGCTCAGTCGCAAGACTTGGCCTTTCGATTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACACCCGCCATACACTGCCAGGCCTCCAATTAAGCCGAAGGCCCCCCTGACGGATGGCCTTTTTGCGATTATACCAACTCTTCCTGTAGTCATATCTCCAAGCCAt  >  2:56731/1‑251 (MQ=14)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                         
TTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCAT  >  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: