Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,333,218	A→G	R78R (CGT→CGC)	rplU ←	50S ribosomal subunit protein L21

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,333,218	0	A	G	100.0%	127.1 / NA	37	R78R (CGT→CGC)	rplU	50S ribosomal subunit protein L21
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/9);  total (28/9)

TGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGACGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGATCATCAGCACTTCAGCGAACT  >  NC_000913/3333092‑3333351                                                                                                                               |                                                                                                                                      tgGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGa                                                                                                                          >  1:529297/1‑140 (MQ=255)        gccAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTaa                                                                                                                             >  2:134052/1‑130 (MQ=255)        gccAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTaa                                                                                                                             <  1:134052/130‑1 (MQ=255)           aGCCTTTTTATGTGCCATTTGAGATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACttt                                                                                                                <  1:196024/140‑1 (MQ=255)                       tgCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGAcc                                                                                                    >  2:403760/1‑140 (MQ=255)                         ccATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGt                                                                                                  >  1:359637/1‑140 (MQ=255)                            tttGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAg                                                                                               >  1:522075/1‑140 (MQ=255)                                    ctctcCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAAGGATTTTACCTTACTAGCCACGCCCGTGAGCAACAACt                                                                                       >  1:273907/1‑140 (MQ=255)                                              gTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGa                                                                             >  1:27203/1‑140 (MQ=255)                                              gTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGa                                                                             >  1:473288/1‑140 (MQ=255)                                                 ttAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTa                                                                          >  2:396202/1‑140 (MQ=255)                                                                          acaTCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTAAGCAACAACTTCCGCTGTGATTACGACGCCGTCAACGAAAGGAACGcc                                                 >  2:277422/1‑140 (MQ=255)                                                                            aTCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGa                                               >  2:118593/1‑140 (MQ=255)                                                                            aTCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGa                                               >  2:114720/1‑140 (MQ=255)                                                                                gtcaaCCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGAtttt                                           <  2:404130/137‑1 (MQ=255)                                                                                  gAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGa                                         <  1:118593/140‑1 (MQ=255)                                                                                        cTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAACTTTGATTGCGCGGCCTTCAACGAATGGAATGCCGATTTTGActtctt                                   >  2:184472/1‑140 (MQ=255)                                                                                        cTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCAAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACTCCGCCATGAACCAAAGGAACGCCGATTTTGActtctt                                   >  1:233773/1‑140 (MQ=255)                                                                                         tGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACttcttc                                  >  1:537023/1‑140 (MQ=255)                                                                                         tGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACttcttc                                  >  2:504131/1‑140 (MQ=255)                                                                                           aCGATGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCAc                                <  2:473288/140‑1 (MQ=255)                                                                                               tGGCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGACCGCCGATTTAGACTTCTTCACCGtt                            >  1:288040/1‑140 (MQ=255)                                                                                                ggCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGttt                           >  2:383947/1‑140 (MQ=255)                                                                                                 gCCCTGCTGCTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTg                          <  2:522075/140‑1 (MQ=255)                                                                                                     tgctgcTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGat                      >  1:545/1‑140 (MQ=255)                                                                                                     tgctgcTTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAAAGGAAGGAACGCCGATTTTGACTTGTTCACCGTTTGCGat                      >  1:257336/1‑140 (MQ=255)                                                                                                          ctTACGATAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGAtcatca                 >  2:647361/1‑140 (MQ=255)                                                                                                                 tAGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGATCATCAGCACTTc          >  2:473736/1‑140 (MQ=255)                                                                                                                  aGTGTTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGATCATCAGCACTTCa         >  2:569505/1‑140 (MQ=255)                                                                                                                    tgtTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGTAACGCCGACTTTGACTTCTTCACCGATTGCAATCATCAGCACTTCAGc       >  2:275561/1‑140 (MQ=255)                                                                                                                    tgtTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGATCATCAGCACTTCAGc       >  1:406655/1‑140 (MQ=255)                                                                                                                    tgtTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGATCATCAGCACTTCAGc       >  1:13504/1‑140 (MQ=255)                                                                                                                    tgtTTACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGATCATCAGCACTTCAGc       <  2:529297/140‑1 (MQ=255)                                                                                                                      tttACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGATCATCAGCACTTCAGCGa     <  1:504131/140‑1 (MQ=255)                                                                                                                      tttACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGATCATCAGCACTTCAGCGa     >  1:126868/1‑140 (MQ=255)                                                                                                                      tttACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTCGCGATCATCAGCACTTCAGCGa     >  1:511018/1‑140 (MQ=255)                                                                                                                       ttACGACGGCGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGATCATCAGCACTTCAGCGaa    <  1:403760/140‑1 (MQ=255)                                                                                                                         acgacggcgAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGATCATCAGCACTTCAGCGAACt  >  1:1787/1‑140 (MQ=255)                                                                                                                         acgacggcgAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGATCATCAGCACTTCAGCGAACt  >  1:121464/1‑140 (MQ=255)                                                                                                                         acgacggcgAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTATTCGCCGTTTGCGATCAGGAGCACTTCAGCGAACt  >  2:290306/1‑140 (MQ=255)                                                                                                                               |                                                                                                                                      TGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAGTGTTTACGACGACGAAACTTAACGATTTTAACTTTCTCGCCACGACCGTGAGCAACAACTTCAGCTTTGATTACGCCGCCATCAACGAAAGGAACGCCGATTTTGACTTCTTCACCGTTTGCGATCATCAGCACTTCAGCGAACT  >  NC_000913/3333092‑3333351

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