Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	3,545,692	C→T	noncoding (2808/2931 nt)	LY180_16870 ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	3,545,692	0	C	T	93.3%	17.1 / NA	15	noncoding (2808/2931 nt)	LY180_16870	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (4/10);  minor base G (1/0);  total (5/10)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACTCTCAAGGAGTCAGGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATCTCTTCCGCATTTAGCTAC  >  CP006584/3545400‑3545772                                                                                                                                                                                                                                                                                                     |                                                                                 aTGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACCCTTaaaggg                                                                            >  1:440556/1‑295 (MQ=14)             cacaCTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGCTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTCAGGTTAAGCCTCAAGGTTCATTAGTACCGGTTAGCTCAACGCATCCCTCCGCCTACACACCCGGCCTATCACCGCCGTCGCCTTCAACGCTCTTTCAGAACCCCTAAAGGGTCCGGGAGAACt                                                                >  2:554989/1‑299 (MQ=11)               cacttcctccGGGGCTTCGGCGGTTCACTTCTGGGTTCGGCGTGGGGTCAGGTGGGACCACCGCGGTAAGGCGGCCCGGCAAATGTTGTTTAATCATCCCGCTTCTGCGTTCTGATTTAATCTGTAGCGGGCTGAAAATCTTCTCTCAGCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTAGCAACGTCGTCGTCTTCAACGTTCCTTCAGGACCCTTAAAGGGTCAGGGAGAACTCa                                                              <  2:466412/290‑1 (MQ=11)               cacTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACCCTTAAAGGGTCAGGGAGAACTCa                                                              >  1:622107/1‑299 (MQ=14)                      aTCGGCGCTACGGCGTTTCACTTCTGAGTTCGGAATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACCCTTAAAGGGTCAGGGAGAACTCATCTCggg                                                       <  1:68349/299‑1 (MQ=9)                             ggagcggggttCACTTGTGAGTTCGGGAGGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCGGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCGGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACCCTTAAAGGGTCAGGGAGAACTCATCTCGGGGCAAGtt                                                <  2:70678/290‑1 (MQ=11)                                     tttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACCCTTAAAGGGTCAGGGAGAACTCATCTCGGGGCAAGTTTCGTGCtt                                        <  1:51765/299‑1 (MQ=14)                                          cTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACCCTTAAAGGGTCAGGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATg                                   <  1:368434/299‑1 (MQ=14)                                                   tCGGCTTGGGATCAGGTGGGGCCCCCGGGCTTAAGCCGCCGGGCATATTCTGTTTTATCAGCCCGCCGCTTCGTTCTGATTTAATTTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACGCCTTCGGCGTGGTAAGGGTAGGCCTCACGGTTCATTGGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACCCTTAAAGGGTCAGGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCa                          <  2:579506/299‑1 (MQ=11)                                                          ggggTCAGGTGGGACGACCGCGCTAAGGCCGCCAGGCAATTTCTGTTTGTTCAGACCGCTTCTCCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCGCTCAGCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACCCTTAAAGGGTCAGGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTAtct                   <  2:91413/299‑1 (MQ=11)                                                          ggggTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACCCTTAAAGGGTCAGGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTAtct                   <  1:14623/299‑1 (MQ=14)                                                             gTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACCCTTAAAGGGTCAGGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATCTCtt                <  1:637210/299‑1 (MQ=14)                                                               cAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCACAGCATCGCTGCGCTTACACACCCGGCCTACCCAAGCCGTCGTCTTCAACGTTACCACACGCCCCTGAAAGGGCCGGGGAGAACCCTTCCCGGGCAACGGTTCGCCCTTCGTTGCTTTCAGCAACTAtctccctc              >  2:472460/1‑295 (MQ=11)                                                                     gagccaccGGGCTAAGGCCGCCAGTCAAATTTTGTTTTAGCCGCCCGCTGCTGCGTTCGTATTTAATCTGTATCAGGCTGAAAGTCTGCGCCCCGCCGCCAAGACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTCCCGGTTAGCTCAACGCATCGCTGCTCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACCCTTAAAGGGTCAGGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATCTCTTCCGCAttt        <  2:354735/296‑1 (MQ=11)                                                                           ccGCGCTAAGGCCGCCAGGCAAATTATGTTTTATCATACCGCTTCTGCGTTCTGATTTAATCTGTATCCGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCCACGCATCGCTGCTCTTACACACCCGGCCTATCAACGTCGCCGTCTTCAACGGTCCTTCAGGACCATTAAAGGGCCAGGGAGAACCCACCTCGGGGCAAGTTTCGTTCTTATATGCTTCAAGCACTTAACTCTTCCGCATTTAGCTAc  >  2:612291/1‑299 (MQ=11)                                                                                                                                                                                                                                                                                                     |                                                                                 ATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGACTCTCAAGGAGTCAGGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATCTCTTCCGCATTTAGCTAC  >  CP006584/3545400‑3545772

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