Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,531	T→C	noncoding (206/2906 nt)	rrlC →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,941,531	0	T	C	100.0%	69.7 / NA	32	noncoding (206/2906 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/17);  total (15/17)

TGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAG  >  NC_000913/3941518‑3941620              |                                                                                          tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGCACGGGGAGCAGCCCAGAGCCTGAATCAGTa            >  1:5811715/1‑93 (MQ=9) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  2:4620526/93‑1 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  2:1276940/93‑1 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            >  2:1133161/1‑93 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  1:670648/93‑1 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACCGGGAGCAGCCCAGAGCCTGACTCAGTa            >  2:4988398/1‑93 (MQ=255)  gAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt           <  2:5121845/93‑1 (MQ=17)   aaaCATCTAAGCGCCCCGAGGAAAAGAAATCAACAGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtg          <  1:5045599/93‑1 (MQ=255)   aaaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtg          >  1:6114504/1‑93 (MQ=18)   aaaCATCGAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGCGGGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtg          <  1:5996505/93‑1 (MQ=11)    gaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgt         <  2:3806525/92‑1 (MQ=18)    aaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgt         <  1:5782973/93‑1 (MQ=18)    aaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgt         >  1:1966398/1‑93 (MQ=18)    aaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgt         <  2:3169189/93‑1 (MQ=18)    aaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgg         >  1:1388026/1‑92 (MQ=16)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCCGAGCCTGAATCAGTAtgtg        >  2:2592410/1‑93 (MQ=255)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtg        <  2:1260096/93‑1 (MQ=18)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtg        >  2:5748358/1‑93 (MQ=18)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtg        >  1:6093635/1‑93 (MQ=18)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtg        <  2:2527628/93‑1 (MQ=18)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtg        <  1:3520279/93‑1 (MQ=18)      cATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgt       <  1:2377803/93‑1 (MQ=21)      cATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgt       >  1:3667831/1‑93 (MQ=18)       aTCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgtg      >  2:3605501/1‑93 (MQ=18)       aTCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgtg      <  1:2945787/93‑1 (MQ=21)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGCGCAGCCCAGAGCCTGAATCAGTATgtgtgt     >  2:2450527/1‑93 (MQ=16)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     >  2:358011/1‑93 (MQ=18)         cTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGtt    <  2:4352543/93‑1 (MQ=21)         cTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGtt    >  2:4458769/1‑93 (MQ=18)         cTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGtt    <  2:4526691/93‑1 (MQ=21)          tAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTa   <  2:6565058/93‑1 (MQ=18)          tAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTa   <  2:2268734/93‑1 (MQ=21)          tAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTa   >  2:1335620/1‑93 (MQ=18)          tAAGCACCCCGAGCAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTa   <  2:6413005/93‑1 (MQ=16)           aaGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAg  <  2:2437004/93‑1 (MQ=21)              |                                                                                          TGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAG  >  NC_000913/3941518‑3941620

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑