Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,019,751	C→T	I197I (ATC→ATT)	ssnA →	putative aminohydrolase SsnA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,019,751	0	C	T	91.3%	56.0 / NA	23	I197I (ATC→ATT)	ssnA	putative aminohydrolase SsnA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/9);  minor base G (0/2);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCGTGAAAGCCACAGGCCGTGGTTTG‑‑CATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTT  >  NC_000913/3019549‑3019963                                                                                                                                                                                                           |                                                                                                                                                                                                                       cGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGAGGACCTGTTGTGAAACTACTGTCCGTAACAACGGCATCAAAGAGTGGCAGGAAGGTGTAGAAGAAAACAGCCGTTTCGCCGGTTTGATTGATGAGGCGAAGAAAGTGACAAGCGTGCCGTATCTGGTGGAAGTTCATATTGGTGCTCATGCGCCGTTTACCGTGGCGGATGCCGGTCTGGAGAtgctg                                                                                                                                                                        <  1:190145/251‑1 (MQ=255)   tataTCGTCGTGTTGCTCTCCCCATTGCGCCATGCATTTTTAAAAGTTGGCCATCGCGCGATGACCTGTTTTGAAACTACATACCGTAACAACGGCTTCAAAGAGTTTCAGTAAGGTGTATAAGAAAACGTCGTTTTCGGCGGTTTGTTTGTTGAGGCGAAGAGATCGACATGCGAACCGTATCTGGGGGAAGGACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCCGTCTGGAGAttctt                                                                                                                                                                        <  1:293140/249‑5 (MQ=255)        cggcggGTCGCTCTCCACATTGCTCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTCCCGTGCCGGATTCCGTTCTGTAGATGCTGCGTGAAg                                                                                                                                                                 <  1:345707/251‑1 (MQ=255)             ggTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGCCATGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTAGCCAGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGAGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCGTg                                                                                                                                                            >  1:235391/1‑251 (MQ=255)                 gccctcCACATTGCGCGACGAATTTTTAAAATGTGGCCTGCGCGCTATGACCTGTTTTGAAACTACTGACCGTAACAACTGCATCAAAGAGTTGAAGGAAGGAGTAGAAGAAAACATCCGTTTCGCCCGGTTGAGGGATGATGCGAATAAAGCTACAAGTGAGCCGTATCTGGTGTAAGTACGGATTTGTTCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCGTGAAAg                                                                                                                                                        <  1:446360/248‑1 (MQ=255)                 gCTCTCCTCATTGCGCGACGCAATTTTAAAAGTTGGCCTGCGCGCGTTGACCTTTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTACAGGAAGTAGTAGAAGAAAACATCCGTTTCGCCCGTTTGAGTGATGAGGCGAAGAAAGCGACAGGCGATCCGTATCTTGTTGAAGCACTTATGGGTTCTCACGCGCCTTTTACCGTGGCGGATGCCGGGCTGGAGATGCTGCGTGAAGCCGTGAAAg                                                                                                                                                        <  1:446359/251‑1 (MQ=255)                      ccACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCGTGAAAGCcac                                                                                                                                                    >  1:58286/1‑250 (MQ=255)                      ccACATTGCGCGACGCATTTTTAAAAGTTGGCCGGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTGGCAGGAAGGTGTAGAAGAAAACAGCCGCTTCGCCCGTTTGAGGGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACAGATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCGTGAAAGCca                                                                                                                                                     <  1:394085/249‑1 (MQ=255)                                            aaaaGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCTTGAAGCCGTGAAAGCCACAGGCCGTGGTTTG‑‑CATATTCAc                                                                                                                            <  1:288395/250‑1 (MQ=255)                                                                  aCCTGTTTTGTAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGATGCGAAGAAAGCGACAAGCGAGCCGTATCTGGGGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGCGAAGCCGTGAAAGCCACAGGCCGTGGGTTG‑‑CATATTCACGCTGCGGAAGACCTTTACGACGt                                                                                                     >  1:409433/1‑251 (MQ=255)                                                                                                               gtatgcctttcgattctaaaCAGCCGTTTCCCCGGTTTGAGTGATGAGGCGAAAAAAGCAACAAGCGCGCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGTTGCGGCGTGAACCCGTGAAAGCCACAGGGCGTGGTTTG‑‑CATATTCACTCGGCGGAAGACCGGTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGTCACGAc                                                         <  1:420049/233‑1 (MQ=255)                                                                                                                 ggAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGACGAAAGCGACAAGCGCGCCGTATCTGGTGGACGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCGTGAAAGCCACAGGCCGTGGTTTG‑‑CATATTCACGCTGCGGACGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTgg                                                      >  1:388505/1‑251 (MQ=255)                                                                                                                      gtgtAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCgg                                                                                                                                                                                     >  1:211217/1‑121 (MQ=255)                                                                                                                          agaagaAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGCAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGcc                                                                                                                                                                                              >  1:430295/1‑108 (MQ=255)                                                                                                                                      ncGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCg                                                                                                                                                                                             >  1:118229/2‑97 (MQ=255)                                                                                                                                         tttCGCCCGTTTGATTGAGGAGGCGAAGAAAGCGACAAGCGACCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGGGCCGGATGCCGGTCTGGAGATGCTTCGGGAAGCCGGTAAAGCCACAGGTCGTGGGTTG‑‑CATATTCACGCTGCGGAAGACCTTGACGACGTTTCCTACAGTCACCACTGTTACGGCAAAGACCGGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGc                               <  1:294003/250‑1 (MQ=255)                                                                                                                                                     gnttgatGAGGAGAAGCAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCGTGAAAGCCACAGGCCGTGGTTTG‑‑CATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAACGCTGGTCg                  >  1:348838/1‑251 (MQ=255)                                                                                                                                                          atgaGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCGTGAAAGCCACAGGCCGTGGTTTG‑‑CATATTCACGCTGCGGAAg                                                                                                                  >  1:234147/1‑150 (MQ=255)                                                                                                                                                                   agaaAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCGTGAAAGCCACAGGCCGTGGTTTG‑‑CATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAAc                          >  1:243634/1‑229 (MQ=255)                                                                                                                                                                     aaaGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGCAGCCGTGAACGCCACAGGCCGTGGTTTG‑‑CATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACtt  >  1:101134/1‑251 (MQ=255)                                                                                                                                                                            cAAGCGAGCCGTATCTAGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCGTGAAAGCCACAGGCCGTTGTTTG‑‑CATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGt                                                                                 <  1:450477/165‑1 (MQ=255)                                                                                                                                                                                   gCCGTATCTGGTGGAAGCACATATTGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTgg                                                                                                                                                                                >  1:345179/1‑65 (MQ=255)                                                                                                                                                                                           tggtggAAGCACAGATGGGTGCTCACGCGCCGTTTACCGTGTCGGATGCCGGTCTGGAGATTCTGCGGGAAGCCGGTAAATCCACAGGCC‑‑GGTTTGTTCATATTCACGCTGGGGAAGACCTTTACGACTTTTCATACAGTCACCACTGGTACGGCAAAGACctg                                                                   <  1:309569/164‑1 (MQ=255)                                                                                                                                                                                                           |                                                                                                                                                                                                                       CGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCGTGAAAGCCACAGGCCGTGGTTTG‑‑CATATTCACGCTGCGGAAGACCTTTACGACGTTTCCTACAGTCACCACTGGTACGGCAAAGACCTGCTGGCACGACTGGCGCAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTT  >  NC_000913/3019549‑3019963

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