Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,624,670	C→T	Q40Q (CAG→CAA)	ytjB ←	putative inner membrane protein Smp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,624,670	0	C	T	93.8%	84.7 / ‑3.6	33	Q40Q (CAG→CAA)	ytjB	putative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (13/17);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGC  >  NC_000913/4624424‑4624915                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       gcctcccgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGACGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTc                                                                                                                                                                                                                                                    <  1:80294/244‑1 (MQ=255)           tttttACCGTCGGGCCCCAGCCGGTCGCGCACTTCGACGCTTTCGCCCGAACGTGCGATAAGACCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTggtggt                                                                                                                                                                                                                                           <  2:96470/250‑1 (MQ=255)                     cGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAAc                                                                                                                                                                                                                                <  1:295830/251‑1 (MQ=255)                         cgcCAGCCGGTCGCGAACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTCCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCAAGCTGTGGATTACGTTGTCGCTGGTGGTTTTAACTAAACAAT‑g                                                                                                                                                                                                                           >  1:187333/1‑251 (MQ=255)                                     gcgcACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGGTCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGc                                                                                                                                                                                                                >  1:43477/1‑250 (MQ=255)                                                    ttCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAAc                                                                                                                                                                                                 <  1:64123/250‑1 (MQ=255)                                                       gCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGCGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGAGCCAGGGAGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCAGCAACGCCCCTAACaggg                                                                                                                                                                                             >  1:55798/1‑248 (MQ=255)                                                             acgggCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGAGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACACGGCAAGAc                                                                                                                                                                                       >  1:58416/5‑251 (MQ=255)                                                              cgggCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGCTCGAGGATCGCCTGAATGCGTTTTTCATACGGCGAGTCGGTACGCATCAGCGGGGCAACGTTCAGCGTTACCTGACGCCCCAGGGGGCGGGCCCGTTCATCCAGCTGTGGATGACGTTTTCGCTGGTGGTTTTGCATAACCCATTGAC‑CTCCCAGCATCCGCGCCCCAACCAAGGCCAGAc                                                                                                                                                                                       >  2:195673/4‑250 (MQ=255)                                                              cgggCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACa                                                                                                                                                                                      <  1:332829/248‑1 (MQ=255)                                                              cgggCGATAAGATAGACTTGTTCGTCCTACACGCCCGCGACGAGGATACGGCTTTCATACGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCCGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCAAGCTGTGGATTACGTTGTAGCTGGTGGTTTTGCCTAAACCAT‑GACGCTCCCTGCATAAGCGACACTACTAAGGCAAgcc                                                                                                                                                                                       >  2:321267/4‑248 (MQ=255)                                                                    agaaGATCGCTTTTTTCGTCATACACGCCCGCGGCTATGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCACCCGGCGAGTCGGTACGCATCAGCGTTGCAACTTTCAGCGTTACCTGACGCGTCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAg                                                                                                                                                                                <  2:324571/249‑1 (MQ=255)                                                                                   tCGTCATACTCTCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACTTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCGGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACACGGCAATACAGAACAGGACAATCTCTGccc                                                                                                                                                                  <  2:307536/250‑1 (MQ=255)                                                                                    cGTCATACACACCCGCGTCGAGGATACGTCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCTTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACCAACAAGGCAAGACAGAACAGGACGATCACTGCCCgg                                                                                                                                                                >  1:41237/1‑250 (MQ=255)                                                                                                                   ttCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGcc                                                                                                                                  >  1:258307/1‑250 (MQ=255)                                                                                                                        ccGTTAACTGAGCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGTTTCAACGTTCAGCGTTACCTGACGCGGCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAGCAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCtt                                                                                                                             <  1:31330/250‑1 (MQ=255)                                                                                                                               gcggtgaaGGATCGCTTGAATGCGTTTTACATCCGGCGATTCGTGAGTTGTCAGCGGTTCCACGTTGAGCTTTACCTGACGGGCCTGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTGTTGACTAACCCAT‑GACGCTCCCTGCAGCAGCGCCACTAAAAAGGCAAGTCAGAACAGGACAATCACTGCCCGATGCATCCGGAATGTCAGTTTTGTGCGAGCTATCTTCCAACCtt                                                                                                                     <  2:245741/244‑1 (MQ=255)                                                                                                                                 gATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTg                                                                                                                   <  1:211432/251‑1 (MQ=255)                                                                                                                                       ggATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAAtt                                                                                                             <  1:142685/251‑1 (MQ=255)                                                                                                                                               tGAATGCGTTTTTCATCCGGCGAGTCGGTACGCCTCAGCGGTGCAACGTTCAGCGTTACCTGAAGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGAGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACtt                                                                                                     >  1:324834/1‑251 (MQ=255)                                                                                                                                                gAATGCGTTTTTAATCCAACAAGTCGGTACGCATCAACGGTGCAACGTGCAGCGTTACCTGGCGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTACACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACCGGACAATAACTGCCCGCTGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTgagaattt                                                                                                    >  2:257117/1‑247 (MQ=255)                                                                                                                                                          ttCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCa                                                                                           <  1:195673/250‑1 (MQ=255)                                                                                                                                                                              gCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGCTGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGt                                                                        >  1:363087/1‑249 (MQ=255)                                                                                                                                                                                       gTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCgagatcggaagagc                                                                                                                                           >  2:167567/1‑160 (MQ=255)                                                                                                                                                                                       gTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCgagatcggaagagc                                                                                                                                           <  1:167567/173‑14 (MQ=255)                                                                                                                                                                                          cAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGtt                                                           >  1:202867/1‑250 (MQ=255)                                                                                                                                                                                                        tccgtaCGCGCTAGGTTGCGTGCCAGTTCTTCCAGCTGTTGATTACGTTGTCGCTGGTGGGTTTGACTAAACCAT‑GACGCTCCCTGCTTCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCTGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTAGTTTCCCGGCTTc                                             <  2:336467/245‑1 (MQ=255)                                                                                                                                                                                                             aCGCGGCAGGGTGCGGGCCAGTTCTTCCTGCTGTTGATGCCTGTGTCGCTGGTGGTTTTGACTAAACCAT‑GGCGCTCCCTGCATCAGCTACACTAACAAGGCAAGACAGAACAGGACAATCACAGCTCGATGCCGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATTCGTTATTTTCCCTGCTTCGAACGa                                       <  2:43477/251‑1 (MQ=255)                                                                                                                                                                                                                  ccAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTa                                  <  1:106966/251‑1 (MQ=255)                                                                                                                                                                                                                                        cAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACAt            <  1:321267/251‑1 (MQ=255)                                                                                                                                                                                                                                              tgGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTg                                                                                              <  1:279263/163‑1 (MQ=255)                                                                                                                                                                                                                                              tgGATTACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTg                                                                                              >  2:279263/1‑163 (MQ=255)                                                                                                                                                                                                                                                   tACGTTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGc  <  1:171490/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       GCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCAT‑GACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGC  >  NC_000913/4624424‑4624915

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