Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,943,846	C→T	noncoding (143/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,846	0	C	T	88.5%	28.3 / 1.7	26	noncoding (143/2904 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (12/11);  total (13/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 = 6.97e-01

GGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGG  >  NC_000913/3943748‑3944092                                                                                                   |                                                                                                                                                                                                                                                       ggCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGtt                                                                                                  >  2:219005/1‑249 (MQ=14)            aCGTGCTAATCTGCGATAAGCGTCGGTAAGGTGAGATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGGTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCt                                                                                     <  2:635926/251‑1 (MQ=14)                          gATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGcgcgcg                                                                        >  1:290706/1‑250 (MQ=14)                           aTAAGCGGCGGTAAGGTGATATGAAAAGGTATAACCGGCGATTACCGAATGGGAAAACGCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAAGGAGGCGAACCGGGGGAACTGAAACATCTAAGTACGCCGAGGAAAAGAAATCAACCGAGATTACTCCAGTAGCGGCGAGCGAACGGGGACAAGCGAAAAGACTGGATAAGTATGAGTGTTTGTGGGAGCGTCTGGAAAgcggatcggt                                                                      >  2:18502/1‑242 (MQ=2)                              aGCGTAGGTAAGGAGATACGAACCGTTATAACCGGCGATTGTCGAATGTGGAAACCCTGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACCGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGGAGCGGCGAGCGAACGGGGAGCAGCCCAAAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGAGATACa                                                                   <  2:327564/251‑1 (MQ=11)                                   cGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGCTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCACAGTAGCGGCGCGCGAACGGGGAGCAGCCCAGAGCCTGAATCGGTACGTGTGTTAGTGGAAGCgataggcaagccgcgcgctacagcg                                                                >  2:240482/1‑225 (MQ=11)                                     gTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGAc                                                            <  2:330492/251‑1 (MQ=31)                                       aaGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACa                                                           <  2:235500/250‑1 (MQ=14)                                       aaGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAg                                                          <  2:352575/251‑1 (MQ=14)                                            gATATGAACCGTTATAACCTGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTTAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAAGCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGTCGCGCGATACAGGGTGACAGcccc                                                      <  1:625647/250‑1 (MQ=9)                                               atGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGGGTGGTAGTGGAAGCGTCTGGCAAGGCGCGCGATACAGGGTGACAGCCCCGTac                                                  >  2:361576/1‑251 (MQ=9)                                                      gTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGGGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACaaaa                                            >  2:180844/1‑250 (MQ=14)                                                        tataACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGGGACAGCCCCGTACACAAAAATg                                         >  2:369783/1‑251 (MQ=14)                                                           aaCCGGCGATTTCCGAATGGGGAATCCCAGTGTGTTTCGATCCACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGCATCAGTATGTGTGTTAGTGGCAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCTTTACACAAAAATGca                                       <  1:145968/250‑1 (MQ=11)                                                           aaCCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGgc                                                                             <  1:215084/212‑1 (MQ=14)                                                           aaCCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGgc                                                                             >  2:215084/1‑212 (MQ=14)                                                                      tCCGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGTCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAg                            <  2:480238/250‑1 (MQ=25)                                                                       ccGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCt                          <  2:153184/251‑1 (MQ=34)                                                                        cGAATGGGGAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTc                         <  1:72821/251‑1 (MQ=21)                                                                             ggggAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGa                    <  1:423888/251‑1 (MQ=21)                                                                               ggAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATg                                         <  1:399607/228‑1 (MQ=14)                                                                               ggAAACCCAGTGTGTTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCACGAGGAAAAGAAATAAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGAAGCCCAGAGCCTGAATCAGAATGGGGGTTAGTGGAAGCGACTGGAAAGGAGAGAGATACAGGGTGACAGCCCCGTACACAAAAATg                                         >  2:399607/1‑228 (MQ=11)                                                                                      cAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCgg           >  2:285015/1‑251 (MQ=34)                                                                                         tgtgtTTCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCgggcc        >  2:259785/1‑249 (MQ=21)                                                                                               tCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCCGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTgg  >  2:14288/1‑251 (MQ=21)                                                                                               tCGATACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAAAGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGGGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGAACACAAAAATGCACATATTGTGAGCTCGATTAGTAGGGCGGGACACGTgg  >  2:38724/1‑251 (MQ=255)                                                                                                   |                                                                                                                                                                                                                                                       GGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGG  >  NC_000913/3943748‑3944092

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