Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,423,701	A→G	noncoding (87/120 nt)	rrfD ←	5S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,701	0	A	G	100.0%	59.8 / NA	20	noncoding (87/120 nt)	rrfD	5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/4);  total (16/4)

ACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTT  >  NC_000913/3423590‑3423798                                                                                                                |                                                                                                  aCCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTACCTACTCTGGCAGGCGTAGGCCCCAAACTAACATCGGCGGTACGGCGt                                                                       >  2:420031/1‑140 (MQ=25)   ccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGttt                                                                     <  2:656121/140‑1 (MQ=255)       agTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCCCACAACCATCGGCGCTACGGCGTTCCACt                                                                 >  1:393360/1‑140 (MQ=255)         tCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                                               >  2:47259/1‑140 (MQ=255)           gAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                                                             >  2:811443/1‑140 (MQ=255)            aaCTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                                            >  1:705256/1‑140 (MQ=255)                   gACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc                                                     >  2:744001/1‑140 (MQ=255)                                     ggTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGacc                                   >  2:680651/1‑140 (MQ=255)                                         cgcTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGt                                              <  2:83578/125‑1 (MQ=255)                                         cgcTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGt                                              >  1:83578/1‑125 (MQ=255)                                         cgcTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcg                               >  1:206732/1‑140 (MQ=255)                                                        gCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAg                >  2:35365/1‑140 (MQ=255)                                                        gCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAg                >  1:759242/1‑140 (MQ=255)                                                        gCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAg                >  1:745160/1‑140 (MQ=255)                                                        gCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAg                >  1:537416/1‑140 (MQ=255)                                                          cATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGc              >  2:66966/1‑140 (MQ=255)                                                              tCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAAt          <  1:744001/140‑1 (MQ=38)                                                                      cTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAAt          >  1:90534/1‑132 (MQ=17)                                                                      cTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAAt          <  2:90534/132‑1 (MQ=17)                                                                      cTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGtttt  >  1:123494/1‑140 (MQ=17)                                                                                                                |                                                                                                  ACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTT  >  NC_000913/3423590‑3423798

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