Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,400,254 C→T intergenic (+52/‑34) miaA → / → hfq tRNA dimethylallyltransferase/RNA‑binding protein Hfq

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,400,2540CT96.7% 104.7 / NA 30intergenic (+52/‑34)miaA/hfqtRNA dimethylallyltransferase/RNA‑binding protein Hfq
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (13/16);  minor base G (0/1);  total (13/17)
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.98e-01

ATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTC  >  NC_000913/4400011‑4400493
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               
aTGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTAc                                                                                                                                                                                                                                          <  1:197543/251‑1 (MQ=255)
                aaGGCGAAATCTCATACGCTGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACCAATAAGCATATAAgg                                                                                                                                                                                                                          >  2:151202/1‑251 (MQ=255)
                  ggctaaaTCTCATCCGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGGTGGCGAAGCGGCATATAAGCTGGCTGCTTGTTTGGGAAGGGGTTCACTGGCTTGACAGTGAAATATCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGTCTGAATGTGTGCAAGTGAGACTTATCGTGCGCAATGTTTTCATGATCGAAAGGTTGAAAGTATAAATTAGCATATAAGGaa                                                                                                                                                                                                                        <  2:44260/247‑1 (MQ=255)
                     gAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACTATACAGTTGGCGAAGCGTCAGATAACCTTGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGTCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGTAAAag                                                                                                                                                                                                                     <  1:234203/251‑1 (MQ=255)
                                    gAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAgg                                                                                                                                                                                                      >  1:75068/1‑251 (MQ=255)
                                     aaaTGGTTTATCGAGGTGTTTGCGCAACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAggg                                                                                                                                                                                                     >  2:277817/1‑251 (MQ=255)
                                                   ggTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAg                                                                                                                                                                                       <  1:92964/251‑1 (MQ=255)
                                                     tgtttgtcccaacacgcAGTAGGCGAAGCGGCCGATAACTCGACTGAGTGGTCGGGCAGGGGCTCACTGGCTTGAAAGTGAAAAACCAGACCAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGGTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTTCGACATTTTTTCCGAATCGAAAGGTTTAAAGTACAAATAAACATCTCAGGAAAAGAGAGCATGGCTAATGGGCAATCtttcca                                                                                                                                                                                         >  2:36281/17‑244 (MQ=255)
                                                                   gaCAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGCAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGACAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCAGTTCCTGAACGc                                                                                                                                                                       >  2:272296/1‑251 (MQ=255)
                                                                             gAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCACTTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCgg                                                                                                                                                             >  2:27878/1‑251 (MQ=255)
                                                                                  ggCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACg                                                                                                                                                        >  1:73825/1‑251 (MQ=255)
                                                                                                  gCGTGGTTGGGAAGGGGTTCCCTGGCATGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTGCCAGTTTCtatt                                                                                                                                        <  2:73825/251‑1 (MQ=255)
                                                                                                   gggggTTGGGAAGGGGTTCACTGGCTTTACAGTGAAAAACCAGAACAGGCGCGTTACGAAGTATTACAGGGTGTTTGTGCTATTGCAGGCTGAATGTGTAAAATTGAGACGTAACGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGGAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtattt                                                                                                                                       <  2:287914/248‑1 (MQ=255)
                                                                                                                     cACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTa                                                                                                                     <  2:169389/251‑1 (MQ=255)
                                                                                                                             tGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCa                                                                                                              >  2:45518/1‑250 (MQ=255)
                                                                                                                                    gAAAAAACAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAATGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTTTATTTATTTGGGGAATGGTTTTAAGCttcatgtgcaaa                                                                                                      <  2:301268/251‑12 (MQ=255)
                                                                                                                                             gAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCtt                                                                                             <  1:153776/251‑1 (MQ=255)
                                                                                                                                               aCAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCtttt                                                                                           >  2:2192/1‑251 (MQ=255)
                                                                                                                                                   gcgcGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                                       <  1:81114/251‑1 (MQ=255)
                                                                                                                                                                   tACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGtt                                                                       >  1:227553/1‑251 (MQ=255)
                                                                                                                                                                    aCAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGtt                                                                       <  2:75068/250‑1 (MQ=255)
                                                                                                                                                                                      cGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCa                                                     <  1:45518/250‑1 (MQ=255)
                                                                                                                                                                                      cGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCa                                                     <  1:2192/250‑1 (MQ=255)
                                                                                                                                                                                              taaTGTGTACAATTGAGACGTATCGTGCGCAATATTTTCAGAATCGAAAGGGTTAAAGTACAAATAAGCATATAAGGAAAAGAGATAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCCGCCAGATgggtt                                             <  2:239378/249‑4 (MQ=255)
                                                                                                                                                                                              gAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGtt                                              <  2:262917/249‑1 (MQ=255)
                                                                                                                                                                                                         aTTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGGGAATGGGATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgcgg                                 >  1:154/1‑250 (MQ=255)
                                                                                                                                                                                                                  tATCGGGCGCAATTTTTTCAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCCAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                        >  1:161996/1‑251 (MQ=255)
                                                                                                                                                                                                                     cGTGCGCAATGTTTACAGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTGTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTGCGTGATCCTGTTGAAAAACACTGTCAGCCAGATGGTTTACAAGCAcgcgctgtcttctgtg                      <  2:200247/250‑13 (MQ=255)
                                                                                                                                                                                                                                     aGAATCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGt     >  1:306972/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aTCGAAAGGTTTAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtc  <  1:151202/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               
ATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTC  >  NC_000913/4400011‑4400493

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