Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,399,673 G→A G141D (GGT→GAT)  miaA → tRNA dimethylallyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,399,6730GA94.4% 37.4 / NA 18G141D (GGT→GAT) miaAtRNA dimethylallyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/5);  minor base C (1/0);  total (13/5)
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.74e-01

CGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGATCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAACTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTG  >  NC_000913/4399478‑4399917
                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
cGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGATCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGATTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAg                                                                                                                                                                                                <  1:289698/250‑1 (MQ=255)
             ggCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGATCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCCAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGCACGGGCCAGCATTGAGCAACAGGCGGCAGAGCAAGATTGGGAGTCATTGCATCGTCCACTTCAGGAGGt                                                                                                                                                                                                                     >  1:374163/1‑216 (MQ=255)
                 gctgATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCTCCGCGGCGGGGCGGATCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGAGCCAGAAGTGAGCAACAGGCGGCAGAGCAAGATTGGGAGTCATTGCTTCTCCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATcc                                                                                                                                                                              <  1:458519/251‑1 (MQ=255)
                    gATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATAACCGCGGCGGGGCGGATCCCACTGTTCGTCGGCGGTACGATGTTGAATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGATTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCCAGGATTCATCCAAATGATCcac                                                                                                                                                                            >  1:483704/1‑250 (MQ=255)
                                aaTGCGCTGGCGGAAATGGCCGATATCACCGCGGCTGGGCGGATCCCACTGTTACTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCAGGCAGTAAGGGCCAGAATTGAGCAACAGGCGGCGGAGCAAGATTGGGAGTCATTGCATCGTCGACTACAGGAGGAAGATCCGGTTGCGGCAGCAAGGATTCATACAAATGATCCAACAAGGCTTTCCCg                                                                                                                                                               >  1:30715/2‑251 (MQ=255)
                                        ngCGGAAATGGCCGATATCACCGCGGCGGGGCGGATCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGATTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATcc                                                                                                                                                                                                               >  1:100815/2‑195 (MQ=255)
                                                                   ggggcggATCCCACTGTGAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGATTGGGAGTCATTGCATCGTGAACTTCAGGAGGTAGATCCGGTTGCGgc                                                                                                                                                                                                      <  1:160071/177‑1 (MQ=255)
                                                                             ccACTGTTCGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGACGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGATTGGAAGTCATTGCATCGTCAACTTCagga                                                                                                                                                                                                                        >  1:448497/1‑149 (MQ=255)
                                                                                                                ggCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGATTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGTGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAgg                                                                               <  1:479113/251‑1 (MQ=255)
                                                                                                                                  gTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACCGGCGGCAGAGCAAGATTGGGAGTCATTGAATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACg                                                                                                                       >  1:212244/1‑193 (MQ=255)
                                                                                                                                           aCCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGATTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCaa                                                                                                                                                                                                 >  1:41553/1‑110 (MQ=255)
                                                                                                                                             cGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGATTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGc                                                  >  1:388283/1‑251 (MQ=255)
                                                                                                                                                cGGCAGACCAGGAAGTACGGGCCAGCATTGAGCAACAGGCGGCAGAGCAAGATTGGGGGTCATTGCCTCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCTCAACGGCTTTCCCGGGCACTGGACGtttttt                                                                                                                                              >  1:482355/1‑156 (MQ=255)
                                                                                                                                                       cccGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGATTGGGAGTCAttn                                                                                                                                                                                                                                         <  1:309801/58‑2 (MQ=255)
                                                                                                                                                       cccGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGATTGGGAGTCATTg                                                                                                                                                                                                                                         >  2:309801/1‑58 (MQ=255)
                                                                                                                                                                                       cggcAGAGCAAGATTGGGAGTCATTGAATCGTCAACTTCAGGAGGTAGATCCGGGTGCGGCAGCAAGGATTCATCCAAATGATCCACACAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATc                                                              >  1:120330/1‑197 (MQ=255)
                                                                                                                                                                                             agCAAGCTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGCTCCGGTTGAGGCAGCAAGGATTCATCCACATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGAGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGACCCGGCGAGCCGTGAACTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTg  >  1:64342/1‑251 (MQ=255)
                                                                                                                                                                                              gCAAGATTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGGTGCGGCAGAAAGGATTCATCCAAATGATCAACCAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAACACTTTAACGGAACTGCCGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAg                                                                       >  1:495669/1‑181 (MQ=255)
                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
CGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGATCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAACTGCTCCATCAACGCATTGAGCAGCGTTTTCATCAGATGTTG  >  NC_000913/4399478‑4399917

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