Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,734	Δ1 bp	noncoding (409/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,734	0	A	.	100.0%	97.5 / NA	29	noncoding (409/2906 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (17/12);  total (17/12)

GTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCT  >  NC_000913/3941663‑3941753                                                                        |                    gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAc                <  2:4581970/76‑1 (MQ=255) gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAc                >  1:4541191/1‑76 (MQ=255)  tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAcc               >  1:1300975/1‑76 (MQ=255)  tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAcc               <  1:975731/76‑1 (MQ=255)  tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAcc               >  2:1571340/1‑76 (MQ=255)   acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCa              <  2:4647221/76‑1 (MQ=9)     acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtc            >  2:963298/1‑76 (MQ=255)     acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtc            >  2:3756445/1‑76 (MQ=255)      caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcc           <  2:1646943/76‑1 (MQ=12)       agaaaTGCACGTATGGTGGGCTGGATGAGTAGGGCGGGACACGTGGTAGCCTGTCTGAATATGGG‑GGGACCAtcct          <  2:745128/74‑1 (MQ=11)       aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcct          >  1:1390670/1‑76 (MQ=255)       aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcct          <  1:1563482/76‑1 (MQ=12)       aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcct          <  1:29753/76‑1 (MQ=255)        aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGCCCAtcctc         >  2:36051/1‑76 (MQ=11)        aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctc         <  1:1367425/76‑1 (MQ=255)        aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctc         <  1:3492460/76‑1 (MQ=255)        aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctc         >  1:2614085/1‑76 (MQ=255)         aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctcc        <  2:1118584/76‑1 (MQ=9)         aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctcc        <  1:3700565/76‑1 (MQ=255)          aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCa       >  2:212294/1‑76 (MQ=255)          aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCa       >  2:2839809/1‑76 (MQ=255)          aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCa       >  1:2616230/1‑76 (MQ=255)            tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCacg     >  2:1751889/1‑74 (MQ=255)            tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg     >  2:2921897/1‑76 (MQ=255)             gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgc    >  2:2205898/1‑75 (MQ=255)              cccaTATTGTGAGCTCGATGAGTAGGGGGGGCCACGTGGTATCCTGTCTGAATATGGG‑GGGCCTGTCCTCCAAGGc   <  2:949923/74‑1 (MQ=1)              cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGc   >  1:1691022/1‑76 (MQ=12)               acaTATTGTGCGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATTGG‑GGGACCATCCTCCAAGGCt  >  1:2706115/1‑76 (MQ=11)               acaTATTGTGAGATCGATGAGAAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCt  >  1:2664888/1‑76 (MQ=11)                                                                        |                    GTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCT  >  NC_000913/3941663‑3941753

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