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	97.9%	99.8 / ‑4.5	48	noncoding (206/2906 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (23/24);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

TGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAG  >  NC_000913/3941518‑3941620              |                                                                                          tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  1:1851865/93‑1 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  1:7522493/93‑1 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  1:5954010/93‑1 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  2:5506873/93‑1 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  1:4367733/93‑1 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  2:6458886/93‑1 (MQ=14)  gAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt           <  1:621103/93‑1 (MQ=17)  gAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt           >  1:6101994/1‑93 (MQ=17)  gAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt           <  1:1656300/93‑1 (MQ=17)  agaaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGGGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt           <  1:7442242/91‑1 (MQ=255)   aaaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtg          >  2:798917/1‑93 (MQ=18)   aaaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtg          >  1:3095720/1‑93 (MQ=18)   aaaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtg          >  1:2344222/1‑93 (MQ=18)   aaaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtg          >  2:7583352/1‑93 (MQ=18)   aaaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATAAGTAtg          >  1:2851860/1‑93 (MQ=255)    aaCATCTAAGCACCCCGAGGAAAAGAATTCAACCGAGATACTGCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgt         <  2:3741570/93‑1 (MQ=11)    aaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgt         >  2:2270176/1‑93 (MQ=18)    aaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgt         >  2:3918761/1‑93 (MQ=18)    aaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGGAtgg         >  2:4479693/1‑92 (MQ=9)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGATCGGGGAGCAGCCCCGAGCCTGAATCAGTAtgag        >  1:7629565/1‑91 (MQ=11)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtg        <  1:6976326/93‑1 (MQ=18)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtg        <  1:5214325/93‑1 (MQ=18)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtg        <  2:7833538/93‑1 (MQ=18)      cATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgt       >  1:3164755/1‑93 (MQ=33)      cATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgt       <  1:7464898/93‑1 (MQ=18)      cATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgt       >  2:6197501/1‑93 (MQ=18)      cATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgt       <  2:631855/93‑1 (MQ=21)      cATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgt       >  2:7875497/1‑93 (MQ=18)      cATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgt       >  2:1802127/1‑93 (MQ=18)      cATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgt       <  1:6626623/93‑1 (MQ=21)      cATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgg       >  1:4118723/1‑92 (MQ=18)      cATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGGAtgtgt       >  2:6197502/1‑93 (MQ=16)       aTCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgtg      >  1:1683639/1‑93 (MQ=18)       aTCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgtg      >  2:3581473/1‑93 (MQ=18)       aTCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgtg      <  2:6806205/93‑1 (MQ=21)        tCTAAGCACCGCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     <  1:6515429/93‑1 (MQ=255)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCTAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     <  1:5412779/93‑1 (MQ=16)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     >  1:7931021/1‑93 (MQ=18)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     >  2:2179839/1‑93 (MQ=18)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     >  2:3081410/1‑93 (MQ=18)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     >  1:6239684/1‑93 (MQ=18)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     <  1:4768596/93‑1 (MQ=21)         cTAAGCCCCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGACCGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGtt    <  1:4330754/93‑1 (MQ=255)         cTAAGCCCCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAAC‑GGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTa   <  1:1842370/93‑1 (MQ=11)         cTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGtt    <  1:7898419/93‑1 (MQ=18)          tAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTa   <  1:5810601/93‑1 (MQ=21)          tAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTa   <  1:1878437/93‑1 (MQ=21)          tAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTa   >  2:2893921/1‑93 (MQ=18)           aaGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGGGCAGCCCCGAGCCTGAATCAGTATGGGTGTTAg  >  1:1964436/1‑93 (MQ=11)           aaGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAg  <  2:3600202/93‑1 (MQ=21)           aaGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAg  <  2:3103645/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 < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑