Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,606,215 (G)8→10 noncoding (72/87 nt) leuP ← tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,2071.G100.0% 34.2 / NA 12noncoding (80/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (6/6);  total (6/6)
*NC_0009134,606,2072.G83.3% 28.5 / 1.1 12noncoding (80/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (5/5);  total (6/6)
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

GGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTT  >  NC_000913/4606090‑4606314
                                                                                                                        ||                                                                                                           
gggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAAccccc                                                                                         >  2:65129/1‑139 (MQ=255)
  gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑gggg                                                                                                         >  1:197853/1‑121 (MQ=255)
  gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAACTTCGCACAGTCATCTTAATTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑gggg                                                                                                         <  2:197853/121‑1 (MQ=255)
                               acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCt                                                           >  2:70795/1‑139 (MQ=255)
                                                       tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCAAGCGCGTCTACCAATTCCGCCAcc                                   >  1:356094/1‑139 (MQ=255)
                                                            aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                             <  2:423986/139‑1 (MQ=255)
                                                             ttCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                             <  1:70795/139‑1 (MQ=255)
                                                                  ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                       >  1:60131/1‑124 (MQ=255)
                                                                  ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                       <  2:60131/124‑1 (MQ=255)
                                                                      cTTCGCACAGTCATGTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTTAATCCCCACGTCCGTAAGAACACTTACACCTGAATCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                    <  2:356094/139‑1 (MQ=21)
                                                                                tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                             >  2:497792/1‑120 (MQ=255)
                                                                                tCATCTTACTTTTTTTGATATCGCCTCCTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                             <  1:497792/120‑1 (MQ=18)
                                                                                                              tttGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt  <  1:164856/119‑1 (MQ=255)
                                                                                                              tttGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt  >  2:164856/1‑119 (MQ=255)
                                                                                                                        ||                                                                                                           
GGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTT  >  NC_000913/4606090‑4606314

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