Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,083,303	C→T	A183V (GCG→GTG)	sbcB →	exodeoxyribonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,083,303	0	C	T	94.1%	112.9 / NA	34	A183V (GCG→GTG)	sbcB	exodeoxyribonuclease I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/17);  minor base G (0/2);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGCGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGTTTCCGGAATGTTTGGAGCATGGCGCGGCAATACCAGCTGGGTGGCACCGCTGGCGTGGCATCCTGAAAATCGCAATGCCGTAATTATGG  >  NC_000913/2083076‑2083539                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               cgacgtcgaAGTCGCACGCAACAGTTTGTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGTATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTATGCCTGGCGCCCGGAAGGAATAAACTGGCCTGTAAATGATGACGGTCTACCGAGCGGTCGCCTTGAGCATTTAACCAAAGCGAATGGGATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACt                                                                                                                                                                                                                          <  1:267898/250‑1 (MQ=255)   acgacgaAGTCACACGCAACATGTTTTAGCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACGCGCGCTGGGATTTACTTGATGTTATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTTGCCTTGAGCATTTAACCAATGCGAATGGTATCGAACATAGCAACGCCCACGATGGGATGGCTGATGGGTACGCCGCTAtg                                                                                                                                                                                                                       <  2:509672/251‑2 (MQ=255)    cgacgaAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTg                                                                                                                                                                                                                      <  2:470653/251‑1 (MQ=255)       cgaAGCCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTCGAGCTGGCAGCATGATAACTCGCGATGGGATTTACTAGATGTTATGCGTGCCTGTTATGCACTGCGCCCGGAAGGAATAAACTGGCCTGACAATGCTGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTACTGACCATAGCAACGCACACGATGTGATGGCTGCTGTGTACGCCACTATTGCg                                                                                                                                                                                                                    >  2:501334/1‑250 (MQ=255)                                            gATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATATACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGtt                                                                                                                                                                             >  2:387056/1‑251 (MQ=255)                                                               ggCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTTCCCAt                                                                                                                                                          <  1:165375/251‑1 (MQ=255)                                                                         tAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAac                                                                                                                                                >  1:404330/1‑251 (MQ=255)                                                                          aaCTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAac                                                                                                                                                >  1:548246/1‑250 (MQ=255)                                                                          aaCTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAac                                                                                                                                                <  1:72388/250‑1 (MQ=255)                                                                                 gcTGGGATTTACTGGCTGTTAAGAGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAACGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGCATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAGACCA‑GCCTGTTTGATTATCTCTTTACCACTCGTAATAAACACAACCTg                                                                                                                                        >  2:75413/1‑251 (MQ=255)                                                                                    gggATTTACTGGATGTTATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGAt                                                                                                                                      <  1:318893/250‑1 (MQ=255)                                                                                    gggATTTACTGGATGTTATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATg                                                                                                                                     <  1:304316/251‑1 (MQ=255)                                                                                           aCTGGATGTTATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGttg                                                                                                                               <  1:13083/250‑1 (MQ=255)                                                                                                    tATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGtt                                                                                                                      <  2:548246/250‑1 (MQ=255)                                                                                                    tATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTc                                                                                                                     >  2:275492/1‑251 (MQ=255)                                                                                                     aTGCGTGCCTGTTATGCCCTGCGCACGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAACCTGATGGCGTTGATTGATGTTcc                                                                                                                    <  2:196769/251‑1 (MQ=255)                                                                                                            ccTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACTCGTCAG‑CCACGCTTGTTTGATTATCGCTTTACCCAGCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATg                                                                                                             <  2:274725/251‑1 (MQ=255)                                                                                                            ccTGTTATGCCCTGCGCCCGGAAGGAAGAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCGTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATAGGATGGCTGATGTGTACGCCACTATTTCGTTTGCTAAGTTGGTAAAAACGCGTCAG‑CTACGCCTGTTTGTTTATCTCTTTACCCATCGTAATAAACACAAACTGAGGGCGTTGATTGACTTTCCGCAGATg                                                                                                             <  2:267214/251‑1 (MQ=255)                                                                                                                ttATGCCCTGCGCCCGGAAGGAATAACCTGGCCTGAACATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑c                                                                                                        >  2:81257/1‑251 (MQ=255)                                                                                                                               ggaaggaaTAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGttt                                                                                         >  1:224677/1‑251 (MQ=255)                                                                                                                               ggaaggaaCAAACTGGCATGAAAATGCTGACGGTCTACCGCGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGttt                                                                                         >  2:256814/1‑251 (MQ=255)                                                                                                                                 aaggaaTAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGTTTcc                                                                                       <  1:411037/251‑1 (MQ=255)                                                                                                                                      aTAAACTGACCTGAAAATGATGACGGTCGACCGAGCTTTCGCCTGGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTAGTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGTTTCCGGAAt                                                                                  <  1:376914/251‑1 (MQ=255)                                                                                                                                       tAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGCTTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGACA‑ACCCTGGTGCACGTTTCCGGactg                                                                                 >  2:307907/1‑248 (MQ=255)                                                                                                                                          aCTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGTTTCCGGAATGttt                                                                              <  2:416795/251‑1 (MQ=255)                                                                                                                                             ggCCTGAAAATGATGACGGTCTACCGAGCTTTCGTCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATTAAT‑CCCCTGGTGCACGTTGCCGGAATGTTTGGa                                                                           <  1:433524/251‑1 (MQ=255)                                                                                                                                                         atgaCGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGTTTCCGGAATGTTTGGAGCATGGCGCGGc                                                               >  1:166599/1‑251 (MQ=255)                                                                                                                                                              ggTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGACTGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGTTTCCGGAATGTTTGGAGCATGGAGCGGCAATAc                                                          >  2:473848/1‑251 (MQ=255)                                                                                                                                                                        gctgtCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGTTTCCGGAATGTTTGGAGCATGGCGCGGCAATACCAGCTGGGTg                                                <  2:274147/247‑1 (MQ=255)                                                                                                                                                                                  tGAGAATTTAACCAAAGCGAATGGTATTGAACCTAGCAACGCCCACGATGTGATGGCTGCTGTGTACGCCACTATTGCGATGGCAACGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGCTGGCGTTGATTGATGTTCCGCAGATGAAC‑CCCCTGGTGAACGTTTCCGGAATGTTTGGAGCATGGCGCGGCAATACAAGCTGGGGGGCACCGCTgg                                      >  2:353295/1‑251 (MQ=255)                                                                                                                                                                                           aaCCGAAGCGAATGGTATTGAACATAGCAACGCCCGCGATGTGATGTCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAATCGCGTCAG‑CCACGCCTATTTGATTAGCTCGTTACCCATCGTAATAAGCACAAGCTGAGGGCGTTGATTGATGTGCCGCAGATGAAGGCCCTTGGT‑CACGTTTCGGGAATGTTTGCAGCATGGCGCGGCAAAACGAGCTGTGTGGCGCCGCGGGCGGGGCATc                             <  2:200536/251‑1 (MQ=255)                                                                                                                                                                                                    gAATGGTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGTTTCCGGAATGTTTGGAGCATGGCGCGGCAATACCAGCTGGGTGGCACCGCTGGCGTGGCATCCTGAAAATc                    >  2:407902/1‑251 (MQ=255)                                                                                                                                                                                                        ggTATTGAACATAGCAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGTTTCCGGAATGTTTGGAGCATGGCGCGGCAATACCAGCTGGGTGGCACCGCTGGCGTGGCATCCTGAAAATCGCa                 <  1:387056/250‑1 (MQ=255)                                                                                                                                                                                                                      cAACGCCCACGATGTGATGGCTGATGTGTACGCCACTATTGCGATGGCAACGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGTTTCCGGAATGTTTGGAACATGGCGCGGCAATACCAGCTGGGTGGCACCGCTGGCGTGGCATCCTGAAAATCGCAATGCCGTAATTAtgg  >  2:405402/1‑251 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               CGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGATCCTTACGCCTGGAGCTGGCAGCATGATAACTCGCGCTGGGATTTACTGGATGTTATGCGTGCCTGTTATGCCCTGCGCCCGGAAGGAATAAACTGGCCTGAAAATGATGACGGTCTACCGAGCTTTCGCCTTGAGCATTTAACCAAAGCGAATGGTATTGAACATAGCAACGCCCACGATGCGATGGCTGATGTGTACGCCACTATTGCGATGGCAAAGCTGGTAAAAACGCGTCAG‑CCACGCCTGTTTGATTATCTCTTTACCCATCGTAATAAACACAAACTGATGGCGTTGATTGATGTTCCGCAGATGAAA‑CCCCTGGTGCACGTTTCCGGAATGTTTGGAGCATGGCGCGGCAATACCAGCTGGGTGGCACCGCTGGCGTGGCATCCTGAAAATCGCAATGCCGTAATTATGG  >  NC_000913/2083076‑2083539

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