Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,323,798	G→A	K193K (AAG→AAA)	rnm →	RNase AM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,323,798	0	G	A	96.8%	94.4 / NA	31	K193K (AAG→AAA)	rnm	RNase AM
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/18);  minor base C (1/0);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTGCCGAAAGACTTGAAAAAGCGCAAATACCTGGGGCGCTGGAAGGCGCACAACGACTGGCGCAGGGCGGCGCAGTGACGCGCGGTCATTTTGCGCGTTTCCTTGTTGAGTGCGGCAAAGCCAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAGGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGC  >  NC_000913/1323558‑1324041                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    gATTGCCGAAAGACTTGAAAAAGCGCAAATACCTGGGGCGCTGGAAGGCGCACAACGACTGGCGCAGGGCGGCGCAGTGACGCGCGGTCATTTTGCGCGTTTCCTTGTTGAGTGCGGCAAAGCCAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCt                                                                                                                                                                                                                                             <  1:135912/249‑1 (MQ=255)                    aaaGCGCAAATACCTGGGGCGCTGGAAGGCGCACAACGACTGGCGCAGGGCGGCGCAGTGACGCGCGGTCATTTTGCGCGTTTCCTTGTTGAGTGCGGCAAAGCCAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAACACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAACGCGGTGCTAGCTCATCCTGGGCGGTACa                                                                                                                                                                                                                         >  1:366635/1‑250 (MQ=255)                                   ggggCGCTGGAAGGCGCACAACGACTGGCGTAGGGCTGCGCAGTGACGCGCGGTCATTTTGCGCGTTTCCTTGTTGAGTGCGGCAAAGCTAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCTGTGCTAGCTCATCCTGGGCGGTACAATCTGAGTTCTAAATg                                                                                                                                                                                                         <  2:103068/251‑1 (MQ=255)                                                  gcACAACGACTGGCGCAGGGCGGCGCAGTGACGCGCGGTCATTTTGCGCGTTTCCTTGTTGAGTGCGGCAAAGCCAGTTCAATGGAGGATGTCTTTAAACAGTATATGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGAACAATCTTAGTGCTAAATGGCTGAAAAGATTggg                                                                                                                                                                                          >  1:622520/1‑250 (MQ=255)                                                             ggCGCAGGGCGGCGCAGTGACGCGCGGTCATTTTGCGCGTTTCCTTGTTGAGTGCGGCAAAGCCAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCAtttt                                                                                                                                                                                <  1:518717/250‑1 (MQ=255)                                                              gcgcAGGACGGCGGAGTGCCGCGCGGACATTTTGCGCGCTTCCTTGTTGAGTGCGGCCAAGACAGTTCAATGGCGGATGTCTTTAAAAAGAATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCCCAGTGGTGTCCAATAGAACAAGCTATTGATGTCATTCATAATTCTGGCGGAAAAGCGGTGCTAGCTCATCCTGGGCGGGACAATCTTAGTGCTAAATGGGTGAAAAGATTGGTTGCGCAttttt                                                                                                                                                                               >  2:419983/1‑249 (MQ=255)                                                               cgcAGGGCGGCGCAGTGACGCGCGGTCATTTTGCGCGTTTCCTTGTTGAGTGCGGCAAAGCCAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGcc                                                                                                                                                                             <  1:168031/251‑1 (MQ=255)                                                                  aGGGCGGCGCAGTGACGCGCGGTCATTTTGCGCGTTTCCTTGTTGCGTGCGGCAAAGCCAGTTCAATGGCGGCTGTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGGACAATAGACCCAGCTATTGATGTCATTCATCATTCTGGCGGGAAAGCGGGGCTAGCTCATCCTGGGCGGTCCAATCTTAGGGCTAAATGGCTGAAAAAATTGGTTTCGCATTTTTCCgg                                                                                                                                                                           >  2:453936/1‑249 (MQ=255)                                                                  aGGGCGGCGCAGTGACGCGCGGTCATTTTGCGCGTTTCCTTGTTGAGTGCGGCAAAGCCAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGaa                                                                                                                                                                          <  1:302211/251‑1 (MQ=255)                                                                              ggACGCGCTGTCATTTTGCGCGGTTCCTTGTTGAGTGCGGCAAAGCCAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAAATCCGGATACGTTCCGCCACAGTGGTGGACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCGCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTGGCGCATTTTGCCGAACACCACGGTGAc                                                                                                                                                              <  2:545025/250‑1 (MQ=255)                                                                                       gTCATTTTGCGCGTTTCCTTGTTGAGTGCGGCAAAGCCAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGaa                                                                                                                                                     <  1:445830/251‑1 (MQ=255)                                                                                                  ggTTTCCTTGTTGAGTGCGGCAAAGCCAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTg                                                                                                                                          <  2:521051/250‑1 (MQ=255)                                                                                                      tCCTTGTTGAGTGCGGCAAAGCCAGTTCAATGGCGGATGTCTTTAAAAAGAATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCCCATCCTGGGCGGTACAATCTTAGGGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCAAGGTGACGCGGTGGAAGTCGCGCAGTGcc                                                                                                                                        >  2:438676/1‑249 (MQ=255)                                                                                                                    ggCAAAGCCAGTTCAATGGCGGATGTCTTTAAACAGTATCTGGCGCGCGGGGAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGTAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATg                                                                                                                         >  2:61432/1‑250 (MQ=255)                                                                                                                    ggCAAAGCCAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATg                                                                                                                         >  1:16709/1‑250 (MQ=255)                                                                                                                           ccAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCTGTGCTAGCTCATCCTGGGCGGTACAATCTTAGGGCTAAATGGCTGAAAAGATTGGGTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGGGCCAGCAATCGCCCAATGAACGTAcc                                                                                                                 <  1:592138/251‑1 (MQ=255)                                                                                                                                            gTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGTGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTg                                                                                                 <  1:453936/250‑1 (MQ=255)                                                                                                                                            gTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGCAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTgc                                                                                                >  1:305144/1‑251 (MQ=255)                                                                                                                                                ttAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGt                                                                                            >  1:488946/1‑251 (MQ=255)                                                                                                                                                   aaaaGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAg                                                                                         <  1:228487/251‑1 (MQ=255)                                                                                                                                                    aaaGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGc                                                                                        <  1:71295/251‑1 (MQ=255)                                                                                                                                                                gcgcgGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCAGCATGCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATggg                                                                            <  2:509383/251‑1 (MQ=255)                                                                                                                                                                   cgGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCAt                                                                         >  2:507552/1‑251 (MQ=255)                                                                                                                                                                          aCCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAACGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAgg                                                                  <  2:294000/251‑1 (MQ=255)                                                                                                                                                                                   gTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGAttt                                                         <  2:290847/251‑1 (MQ=255)                                                                                                                                                                                   gTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGAttt                                                         <  2:305144/251‑1 (MQ=255)                                                                                                                                                                                           acaGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGTCGGGAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGGGACTCGATGGAAGTCTCGCAGTGCCAGCAATCGCCCAATGAACGTAGCCAGCTGTCGGCCCTGGCGCGTCAGCAACATTTATGGGCGTCACAAGGAACTGATTTTCTTCAGc                                                 <  2:366635/251‑1 (MQ=255)                                                                                                                                                                                                ggtgTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGc                                            >  2:438539/1‑251 (MQ=255)                                                                                                                                                                                                          aGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCg                                  >  1:522340/1‑251 (MQ=255)                                                                                                                                                                                                                         tGTCATTCATCATTCTGGCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAAc                   >  1:227538/1‑251 (MQ=255)                                                                                                                                                                                                                                          gCGGTAAAGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGc  <  1:323851/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    GATTGCCGAAAGACTTGAAAAAGCGCAAATACCTGGGGCGCTGGAAGGCGCACAACGACTGGCGCAGGGCGGCGCAGTGACGCGCGGTCATTTTGCGCGTTTCCTTGTTGAGTGCGGCAAAGCCAGTTCAATGGCGGATGTCTTTAAAAAGTATCTGGCGCGCGGGAAAACCGGATACGTTCCGCCACAGTGGTGTACAATAGAACAAGCTATTGATGTCATTCATCATTCTGGCGGTAAGGCGGTGCTAGCTCATCCTGGGCGGTACAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGC  >  NC_000913/1323558‑1324041

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