Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,943,889	G→A	noncoding (186/2904 nt)	rrlC →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,943,889	0	G	A	96.7%	55.9 / ‑3.7	30	noncoding (186/2904 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (16/13);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATAC  >  NC_000913/3943739‑3944136                                                                                                                                                       |                                                                                                                                                                                                                                                         gCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt                                                                                                                                                      <  1:397997/251‑1 (MQ=11)      cAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAAAGGGGAGCAGCCCAGAGCCTGAATCAGGAtgtgtg                                                                                                                                                 >  2:229924/1‑251 (MQ=14)                            aaTCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAgg                                                                                                                          >  1:171204/1‑251 (MQ=14)                            aaTCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATGCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAgg                                                                                                                          <  1:201258/251‑1 (MQ=14)                                 gCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCTGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACGGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCTAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGgc                                                                                                                         <  1:369667/247‑1 (MQ=9)                                 gCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATGCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGGATGTGTGTTAGTGGAAGCGTC‑TGGAAAGgcgc                                                                                                                       <  1:28455/249‑1 (MQ=14)                                        gCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGACTGGGGAAACCCAGAGTGTTTCGAAACACTATCATTAACTGACTCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATAc                                                                                                                >  2:365732/1‑249 (MQ=9)                                                    tGATATGAACCGTTATAACCGGCGATTTCCGCATGTGGACACCCAGTGTGTTTAGCAACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTACGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAAAGGGGAGCAGCCCAAAGCCTGAATCAGTATGTGGGTTAGTGGAAGCGTC‑TGGCAAGGCGCCCGATACATGGTGACAGcccc                                                                                                  >  2:273696/1‑251 (MQ=11)                                                     gATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGcccc                                                                                                  <  1:49929/250‑1 (MQ=14)                                                     gATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCg                                                                                                 >  1:302647/1‑251 (MQ=14)                                                        atGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTac                                                                                              >  1:327464/1‑251 (MQ=14)                                                                   taACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGCATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATgac                                                                                   >  2:351018/1‑249 (MQ=14)                                                                     aCCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACATCCCCGTACACAAAAATGcaca                                                                                 <  2:286566/251‑1 (MQ=14)                                                                     aCCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCAGTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGcaca                                                                                 <  2:133984/251‑1 (MQ=14)                                                                          cGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATtg                                                                            >  2:322866/1‑251 (MQ=18)                                                                            aTTTCCGAGTGGGGAAAACCAGTGTGTTTCGACACACTCTCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATAAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGAGGAGCAGCCCAGAGCCTGAAACAGTATGGGGGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATAAAGGGAGAAAGCCCCGTAAACAAAACAGAACATATtgtg                                                                          >  2:382793/1‑251 (MQ=11)                                                                              tttcGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGa                                                                         <  2:295943/247‑1 (MQ=21)                                                                                    ctGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGAt                                                                  >  2:151323/2‑251 (MQ=21)                                                                                          aaaCCCAGTGTGTTTCGACACTCTATCATTAACTGAATCCATAGGCTAATGAGGAGAACCGAGGGAGCTGAAACATCTAAGTACCCCGGGGAAAAGAAATCAACAGAGATTCCACCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTTTGTATGATAGTGGAAGC‑TCATGGAAAAGCCCGCGATACAGGGTGCCATCCCCGTACAAAAAAATCAACAAATTGTCAGCTTGATGAga                                                              >  2:227290/1‑248 (MQ=11)                                                                                           aaCCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAgg                                                           >  1:93565/1‑251 (MQ=21)                                                                                            aCCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAggg                                                          <  2:9239/251‑1 (MQ=21)                                                                                                        tCGACACCCTATCATTAACTGAATCCATCGGTTCATAAGGCGACCCGAGGGAACTGAAACATCTAAGTAACCCGAGGTACAGAAATCAACCGAGATTCACCCAGTAGCGGCGAGCGCACGGGGAGAAGCCCAAAGCATGATTCCGTAAGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGAGCGATACAGGGTGACAGCCCCGTACCAAAAAATGCACACATTGGGAGCTCGCTGCGGCGGGCGGGacccgctg                                              >  2:218310/1‑245 (MQ=2)                                                                                                        tCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGTACACGTg                                               <  2:142702/250‑1 (MQ=25)                                                                                                               acTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGt                                             <  1:125816/245‑1 (MQ=25)                                                                                                               acTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGCGGt                                             >  2:125816/1‑245 (MQ=25)                                                                                                                 tATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAACAGAACTCAACCGAGCTTCCCCCAGTAGCGGCGCGCGAACGGGGAGCAGACCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGAATCCTGTc                                     >  1:270693/1‑251 (MQ=11)                                                                                                                          cTGAATCCATAGGTTAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAACCCCGGACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgg                            >  2:320101/1‑251 (MQ=21)                                                                                                                                        tAATGAGGCGAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCTAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGAGGAGTAGGTCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAtcctc               <  2:239223/250‑1 (MQ=255)                                                                                                                                                 gAACCGAGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCATTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTaaa     <  1:405132/251‑1 (MQ=21)                                                                                                                                                    ccGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGACCGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGCGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATAc  >  1:252263/1‑251 (MQ=21)                                                                                                                                                       |                                                                                                                                                                                                                                                         GCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTC‑TGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATAC  >  NC_000913/3943739‑3944136

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