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	96.3%	108.8 / ‑2.6	55	noncoding (206/2906 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (24/28);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAG  >  NC_000913/3941518‑3941620              |                                                                                          tGAAACATCTAAGCGCCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  2:3744981/93‑1 (MQ=9) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGGCCAGAGCCTGAATCAGTa            >  2:1711967/1‑93 (MQ=9) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  1:9358890/93‑1 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  2:9437709/93‑1 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            >  1:3393652/1‑93 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  1:7663341/93‑1 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            >  1:393115/1‑93 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            <  2:185508/93‑1 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTa            >  2:5554860/1‑93 (MQ=14) tGAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACCGGGAGCAGCCCAGAGCCTGAATCAGTa            >  1:9842352/1‑93 (MQ=9)  gAAACATCTAAGCACCCGGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt           <  2:6130983/93‑1 (MQ=12)  gAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGCACGGGGAGAAGCACAGAGCCTGAATCAGTAt           >  2:4236103/1‑93 (MQ=11)  gAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt           >  1:9306383/1‑93 (MQ=17)  gAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt           <  1:8189650/93‑1 (MQ=17)  gAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt           >  1:7333960/1‑93 (MQ=17)  gAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt           >  2:691338/1‑93 (MQ=17)  gAAACATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt           <  2:7145081/93‑1 (MQ=17)  aaaaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt           <  1:2571667/92‑1 (MQ=17)   aaaCATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtg          <  2:6094417/93‑1 (MQ=33)   aaaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtg          <  1:7357382/93‑1 (MQ=18)   aaaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtg          >  1:1219914/1‑93 (MQ=18)    aaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgt         >  2:563005/1‑93 (MQ=18)    aaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgt         <  1:9990081/93‑1 (MQ=18)    aaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgt         <  1:8235952/93‑1 (MQ=18)    aaCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgt         <  1:872930/93‑1 (MQ=18)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtg        <  2:4950668/93‑1 (MQ=18)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtg        <  2:9933839/93‑1 (MQ=18)     aCATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtg        <  1:2663443/93‑1 (MQ=18)      cATCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgt       <  2:8053289/93‑1 (MQ=21)      cATATAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgg       >  1:5573464/1‑92 (MQ=255)       aTCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgtg      >  1:1531550/1‑93 (MQ=18)       aTCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgtg      >  1:8121239/1‑93 (MQ=18)       aTCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgtg      >  1:6357541/1‑93 (MQ=18)       aTCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgtg      >  2:3826444/1‑93 (MQ=18)       aTCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgtg      <  1:5854471/93‑1 (MQ=21)       aTCTAAGAACCCCGAGGAAAAGAAATCAACCGAGATTCACCCAGTAGCGGCGAGCGAACCGGGAGCAGCCCAGAGCCTGAATCAGTAtgtgtg      >  2:1731803/1‑93 (MQ=255)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGGGCGAACGGGGGGCGGACCAGAGCCTGAATCAGGAtgttggt     >  2:1710834/1‑89 (MQ=11)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     <  1:5116697/93‑1 (MQ=18)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     <  1:2925487/93‑1 (MQ=21)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     >  1:2987223/1‑93 (MQ=18)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     <  1:3986643/93‑1 (MQ=18)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     >  1:9115165/1‑93 (MQ=18)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     >  1:7653535/1‑93 (MQ=18)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt     <  1:9042169/93‑1 (MQ=21)        tCTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCGGAATCAGTATgggtgt     >  1:796613/1‑93 (MQ=255)         cTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGtt    <  2:685984/93‑1 (MQ=21)         cTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGtt    >  1:2782277/1‑93 (MQ=18)         cTAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGtt    <  1:9712793/93‑1 (MQ=21)          tAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTa   <  2:5148387/93‑1 (MQ=18)          tAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTa   <  1:3411076/93‑1 (MQ=18)          tAAGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTa   <  2:810304/93‑1 (MQ=21)          caaGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTa   <  2:8636968/92‑1 (MQ=18)           aaGTACCCCGAGGAACAGACATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAAATGTGttcg  <  2:8987800/93‑3 (MQ=11)           aaGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAg  <  1:7785886/93‑1 (MQ=21)           aaGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAg  >  2:6603674/1‑93 (MQ=18)           aaGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAg  >  2:6603675/1‑93 (MQ=18)           aaGCACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAg  <  1:8113781/93‑1 (MQ=18)              |                                                                                          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: ‑