Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,437,145	T→A	noncoding (25/111 nt)	micC →	small regulatory RNA MicC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,437,145	0	T	A	90.5%	56.0 / ‑1.4	21	noncoding (25/111 nt)	micC	small regulatory RNA MicC
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (13/6);  total (14/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

CCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAGTGAAAACATCCTAAAGGAGAAGTTCAGTTAGCCGATAAATATCATTTGTTACGCCTTAAGTAAAACCTTAATCAAATAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGATTTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCTTTTTTTTT‑AGAATTGGATAATCCTTATCCAGAGCATTTAATCGGTGTTGCTTTTTTCCGCTTTTCCGGCCATTTGTGCCAGGAAGTCATAGCGTTTTTGCAAATCAGCTGCAGCATCTTTCCATAACTGTTCTGCCACTTCTGGCTGCTGCGAATTCAGCCGAC  >  NC_000913/1436903‑1437387                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    ccTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAGTGAAAACATCCTAAAGGAGAAGTTCAGTTAGCCGATAAATATCATTTGTTACGCCTTAAGTAAAACCTTAATCAAATAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATCTGCCTTTATTGTCACAGaatttatttt                                                                                                                                                                                                                                             >  2:389408/1‑251 (MQ=255)         tattttatttGCCTTTTATCCCCATTTCGGGGAGTGAAAACATCCTAAAGGAGAAGTTCAGTTAGCCGATAAATATCATTTGTTACGCCTTAAGTAAAACCTTAATCAAATAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTgg                                                                                                                                                                                                                                      <  1:28444/250‑1 (MQ=255)                  tGCCTTTTATCCCCATTTCGGGGAGTGAAAACATCCTAAAGGAGAAGTTCAGTTAGCCGATAAATATCATTTGTTACGCCTTAAGTAAAACCTTAATCAAATAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCcattgcat                                                                                                                                                                                                                            >  2:493736/1‑251 (MQ=255)                                       ggAGTGAAAACATCCTAAAGGAGAAGTTCAGTTAGCCGATAAATATCATTTGTTACGCCTTAAGTAAAACCTTAATCAAATAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTTGGCCATTGCATTGCCACTGATTTTCCAACata                                                                                                                                                                                                       <  1:461894/251‑1 (MQ=255)                                         aGTGAAAACATCCTAAAGGAGAAGTTCAGTTAGCCGATAAATATCATTTGTTACGCCTTAAGTAAAACCTTAATCAAATAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACatat                                                                                                                                                                                                      >  2:533927/1‑250 (MQ=255)                                                               aGTTCAGTTAGCCGATAAATATCATTTGTTACGCCTTAAGTAAAACCTTAATCAAATAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGATTTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAA‑GACAAGCCCGAACAgtcgt                                                                                                                                                                              >  1:167750/1‑251 (MQ=255)                                                                   cAGTTAGCCGATAAATATCATTTGTTACGCCTTAAGTAAAACCTTAATCAAATAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCg                                                                                                                                                                           <  1:446488/251‑1 (MQ=255)                                                                      ttAGCCGATAAATATCATTTGTTACGCCTTAAGTAAAACCTTAATCAAATAATTACTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGc                                                                                                                                                                        >  1:271101/1‑251 (MQ=255)                                                                        aGCCGATAAATATCATTTGTTACGCCTTAAGTAAAACCTTAATCAAATAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCtt                                                                                                                                                                      >  1:471714/1‑251 (MQ=255)                                                                         gCCGATAAATATCATTTGTTACGCCTTAAGTAAAACCTTAATCAAAGAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCttt                                                                                                                                                                     <  2:154923/251‑1 (MQ=255)                                                                                               gCCTTAAGTAAAACCTTAATCAAATAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCTTTTTTTTT‑AGAATTGGATAATCCt                                                                                                                                              >  1:177007/1‑251 (MQ=255)                                                                                                                  tCAAATAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTTTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCTTTTTTTTT‑AGAATTGGATAATCCTTATCCAGAGCATTTACTCg                                                                                                                           >  2:168914/1‑251 (MQ=255)                                                                                                                                                                 gATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGACTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCTTTTTTTTT‑AGAATTGGATAATCCTTATCCAGAGCATTTACTCGGTGTTGCTTTTTTCCGCTTTTCCGGCCATTTGTGCCAGGAAGTCAt                                                                             >  2:491900/1‑250 (MQ=255)                                                                                                                                                                            caAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCTTTTTTTTT‑AGAATTGGATAATCCTTATCCAGAGCATTTAATCGGTGTTGCTTTTTTCCGCTTTTCCGGCCATTTGTGCCAGGAAGTCATAGCGTTTTTGc                                                                  >  2:124586/1‑250 (MQ=255)                                                                                                                                                                                aTAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCTTTTTTTTT‑AGAATTGGATAATCCTTATCCAGAGCATTTAATCGGTGTTGCTTTTTTCCGCTTTTCCGGCCATTTGTGCCAGGAAGTCATAGCGTTTTTGCAAATc                                                             >  1:1923/1‑251 (MQ=255)                                                                                                                                                                                            cTTTTAATTTGATATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGTATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCAAACATATAAAAAGCCAAGCCCGACCAGTCGTCCGGGCTTTTTTTTT‑AGAATTGGATAATCCTTATCCAGAGCATTTAATCGGTGTTGCTTTTTTCCGCTTTTCCTGCCATTTGTGCCAGGAAGTCATTGCGTTTTTGCAAATCAGCTGCAGCATc                                                 >  2:299732/1‑251 (MQ=255)                                                                                                                                                                                                tAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAGCATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCTTTTTTTTT‑AGAATTGGATAATCCTTATCCAGAGCATTTAATCGGTGTTGCTTTTTTCCGCTTTTCCGGCCATTTGTGCCAGGAAGTCATAGCGTTTTTGCATATCAGCTGCAGCATCTTTc                                             >  2:474789/1‑251 (MQ=255)                                                                                                                                                                                                 aaTTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCTTTTTTTTT‑AGAATTGGATAATCCTTATCCAGAGCATTTAATCGGTGTTGCTTTTTTCCGCTTTTCCGGCCATTTGTGCCAGGAAGTCATAGCGTTTTTGCAAATCAGCTGCAGCATCTTTc                                             <  1:389408/250‑1 (MQ=255)                                                                                                                                                                                                        tataCGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCTTTTTTTTT‑AGAATTGGATAATCCTTATCCAGAGCATTTAATCGGTGTTGCTTTTTTCCGCTTTTCCGGCCATTTGTGCCAGGAAGTCATAGCGTTTTTGCAAATCAGCTGCAGCATCTTTCCATAACTg                                     >  1:134120/1‑251 (MQ=255)                                                                                                                                                                                                             gTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGATTGTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAA‑GACAAGCCCGAACAGTCGTCCGGGCTTTTTTTTTTAGAATTGGATAATCCTTATCCAGAGCAGTTAATCGGTGTTGCTTTTTTCCGCTTTGCCGGCCATTTGTGCCAGGAAGTCATAGCGTTTTTGCAAATCAGCTGCAGCATCTTTCCATAACTGTTCTg                                <  2:167750/251‑1 (MQ=255)                                                                                                                                                                                                                      cgcgGGTTATATGCCTTTATTGTCACAGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCTTTTTTTTT‑AGAATTGGATAATCCTTATCCAGAGCATTTAATCGGTGTTGCTTTTTTCCGCTTTTCCGGCCATTTGTGCCAGGAAGTCATAGCGTTTTTGCAAATCAGCTGCAGCATCTTTCCATAACTGTTCTGCCACTTCTg                       >  1:566175/1‑251 (MQ=255)                                                                                                                                                                                                                                            tctcaGAATTTATTTTCTGTTGGGCCATTGCATTGCCACTTCTTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCATTTTTTAT‑AGAATTGGATAATCCTTATCCAGAGCATTGAATCGGTGTCGCTTTTTACCGCTTTTCCGGCCATTTGTGCCATGAAGTCATAGCGTTTTTGCAAATCAGCCGCAGCATCTTTCCATAACTGTTCTGCCACTTCTTGCTGCTGCGAATTCAGCTGAc  <  2:574715/247‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    CCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAGTGAAAACATCCTAAAGGAGAAGTTCAGTTAGCCGATAAATATCATTTGTTACGCCTTAAGTAAAACCTTAATCAAATAATTCCTCTGAAATGATAATAAGTCTGGTGAATGTATCGAAGATAATACATACAAAATAAAAATTATACTTTTAATTTGCTATACGTTATTCTGCGCGGGTTATATGCCTTTATTGTCACAGATTTTATTTTCTGTTGGGCCATTGCATTGCCACTGATTTTCCAACATATAAAAAGACAAGCCCGAACAGTCGTCCGGGCTTTTTTTTT‑AGAATTGGATAATCCTTATCCAGAGCATTTAATCGGTGTTGCTTTTTTCCGCTTTTCCGGCCATTTGTGCCAGGAAGTCATAGCGTTTTTGCAAATCAGCTGCAGCATCTTTCCATAACTGTTCTGCCACTTCTGGCTGCTGCGAATTCAGCCGAC  >  NC_000913/1436903‑1437387

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