Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,947,456	(C)5→6	noncoding (70/77 nt)	metZ →	tRNA‑Met

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,451	1	.	C	91.2%	104.1 / 3.5	34	noncoding (65/77 nt)	metZ	tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (14/17);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

ATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGG‑CCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGAT  >  NC_000913/2947314‑2947585                                                                                                                                             |                                                                                                                                       aTGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGG‑cc                                                                                                                                      <  1:180070/140‑1 (MQ=255)                    ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAAt                                                                                                                    >  2:895427/1‑140 (MQ=255)                    ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAAt                                                                                                                    >  2:785779/1‑140 (MQ=255)                    ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGG‑cccc                                                                                                                                    >  1:544880/1‑123 (MQ=255)                    ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGG‑cccc                                                                                                                                    <  2:544880/123‑1 (MQ=255)                        aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTtga                                                                                                                <  2:875979/140‑1 (MQ=255)                        aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTtga                                                                                                                <  2:807544/140‑1 (MQ=255)                             aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAg                                                                                                           <  2:177230/140‑1 (MQ=255)                                gTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTaa                                                                                                        <  1:808996/140‑1 (MQ=255)                                  tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCGACCAATTAAAATTTGATGAAGTAAAg                                                                                                      <  1:428/140‑1 (MQ=255)                                  tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAg                                                                                                      <  1:987697/140‑1 (MQ=255)                                              cgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAc                                                                                          >  1:458777/1‑140 (MQ=255)                                               gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAcg                                                                                         <  1:895427/140‑1 (MQ=255)                                                    gAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCgggg                                                                                    <  1:785779/140‑1 (MQ=255)                                                          aGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGCGATCATCACCCCAAGGCCGTCGGTTCAAAAC‑CCGC‑CCCCCGCAACCAATTAAAATTTGATGAAGTAAACCCGTGCGGTGACGCGGGGTGGagc                                                                              <  2:424117/140‑1 (MQ=17)                                                           gTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTG‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagca                                                                             <  2:429802/140‑1 (MQ=255)                                                             gcgcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGG‑CCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGcc                                                                          >  2:110959/1‑140 (MQ=255)                                                              cgcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGAAAAGCAGTACGGTGACGCGGGGTGGAGCAGcc                                                                          >  2:1113835/1‑140 (MQ=255)                                                               gcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagca                                                                             >  2:444380/1‑136 (MQ=255)                                                               gcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagca                                                                             <  1:444380/136‑1 (MQ=255)                                                               gcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGTAGCCt                                                                         >  1:79812/1‑140 (MQ=255)                                                                 aTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagc                                                                              >  2:1161550/1‑133 (MQ=255)                                                                 aTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagc                                                                              <  1:1161550/133‑1 (MQ=255)                                                                 aTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGG‑CCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGt                                                                      >  1:927921/1‑140 (MQ=255)                                                                          cgcgGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGG‑CCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtcg                                                             <  1:110959/140‑1 (MQ=255)                                                                           gcgGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCACCCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtcg                                                             <  2:572581/140‑1 (MQ=255)                                                                           gcgGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtcg                                                             <  2:458777/140‑1 (MQ=255)                                                                             ggggTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCggg                                                           >  1:572581/1‑140 (MQ=255)                                                                                 tGGATCAGCCTGGTAGATTGTCGGGCAAATAATGCGAAGGTCGTCGGTGC‑AAAT‑CCGGCCCCCCGCAACCAATTAGAATTTGAAGAAGTAAAGCAATCCGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCa                                                       <  2:351131/140‑1 (MQ=14)                                                                                  ggAGCTGGCTGTTAGTTCGTCGGGCAAATAACCCCAAGGTCGTCGGCTC‑AGATACAGA‑CCCCTGCAACCAATTAAAAATCGATGTAGTAAACGAATACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCAt                                                      <  2:978006/140‑1 (MQ=9)                                                                                             gTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGGGGAGCAGCCTGGTAGCTTGTCGACCTCATAACCCgccggc                                           >  2:247271/1‑135 (MQ=14)                                                                                               aGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGG‑CCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtcgt                                        <  2:927921/140‑1 (MQ=255)                                                                                                   cgtcgGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGTGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGCCGTCgg                                     >  1:424889/1‑140 (MQ=255)                                                                                                       gggCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCa                                 >  1:381560/1‑140 (MQ=255)                                                                                                               aaCCCGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGc                         <  1:905086/140‑1 (MQ=255)                                                                                                                   cGAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGAAACCCATTAAAATATGCTCAAGTAAAGCAGTACGGGAACAGGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGccccc                     >  2:905086/1‑140 (MQ=14)                                                                                                                    gAAGGTCGTCGGTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGCCCCCCg                    >  1:1066069/1‑140 (MQ=255)                                                                                                                               gTTC‑AAAT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCGGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaa         <  2:381560/140‑1 (MQ=255)                                                                                                                                    aaaT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAt            <  2:42781/133‑1 (MQ=255)                                                                                                                                    aaaT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAt            >  1:42781/1‑133 (MQ=255)                                                                                                                                      aT‑CCGGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtga   <  1:964039/140‑1 (MQ=255)                                                                                                                                       cccggg‑CCCCCGAAGTCAATAAAAAATAGATGAATTGAACGAATACGGTGACGCGGGGTGGAGTAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat  <  2:79812/136‑1 (MQ=14)                                                                                                                                             |                                                                                                                                       ATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTC‑AAAT‑CCGG‑CCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGAT  >  NC_000913/2947314‑2947585

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