Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,511,466	C→T	V293V (GTG→GTA)	yhfZ ←	putative DNA‑binding transcriptional regulator YhfZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,511,466	0	C	T	96.2%	69.4 / ‑3.5	26	V293V (GTG→GTA)	yhfZ	putative DNA‑binding transcriptional regulator YhfZ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (13/12);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

TTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTGACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGA  >  NC_000913/3511239‑3511710                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     tttCTTCACCACGGCGACTGAGCATAACCGCCATCGACATCTGTGTCATCGCCATCGTTCCTTCCTCAATGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTGACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCg                                                                                                                                                                                                                                 >  2:12921/1‑249 (MQ=255)         ccacGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTGACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCgt                                                                                                                                                                                                                       >  1:193457/1‑251 (MQ=255)              gCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTGACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGttt                                                                                                                                                                                                                   >  1:77522/1‑250 (MQ=255)               cGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTGACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTAt                                                                                                                                                                                                                 >  2:173858/1‑251 (MQ=255)                      gcATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTGACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTTCCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAAc                                                                                                                                                                                                          <  1:117127/251‑1 (MQ=255)                           aGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTGACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCAc                                                                                                                                                                                                     >  2:171706/1‑251 (MQ=255)                                                                           aCCGTCAGATGGCACTCTGTTTCCAGTACGTTGCCGACCTGCATCATGCCGTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCATATGCCTTTCGTTTAATATCTCGGTTCCTGTTCCCCCCTTACTACCCGTTGTTGATGTGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGGTGTTGCATTGGGTAATCATCTACTCTCGTCAGGACCACgg                                                                                                                                                     <  2:147920/251‑1 (MQ=255)                                                                                        aCTCTGTTTCCAGTACGTTGACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTg                                                                                                                                        >  2:14020/1‑251 (MQ=255)                                                                                                              gacCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCtgtg                                                                                                                  >  1:72941/1‑251 (MQ=255)                                                                                                                        atGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCg                                                                                                        <  1:147633/251‑1 (MQ=255)                                                                                                                        atGCCTTTGAAGATGTCCTTATCACTAACGCCTGCCTCGCAAAGCAGCCTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCg                                                                                                        >  2:137184/1‑251 (MQ=255)                                                                                                                                gcagaagccctaacccaTAACGCCTGACTCGCAAAGCCGATTGAGTCTGCTTTCCATATCAGATGACTTTCGTTAAATAACTAGGCTCCTGTTACCCGCTTACTACACGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACACCGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTAGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAAt                                                                                                >  2:156765/17‑251 (MQ=255)                                                                                                                                                  aaCGCCTGCCTCGCATAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTGAATAACTCGGTTCCGGTTCCCCGCTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTCTTTATCTACAACGGCACGCAGCATTTGTTGCAGGGTGTAATCATCGATTCGCGTCAGGACCACTTCTTCGGTGGCCTGTAAAAATCGCGGATCGTTTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTc                                                                               <  1:152624/250‑1 (MQ=255)                                                                                                                                                         gcctCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTcc                                                                        <  2:81027/250‑1 (MQ=255)                                                                                                                                                         gcctCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTcc                                                                        <  1:240454/250‑1 (MQ=255)                                                                                                                                                                   gCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATAGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGATCGTTTTCCGCAACCACGtt                                                             >  1:147012/1‑251 (MQ=255)                                                                                                                                                                     cggTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTcc                                                           >  2:131270/2‑251 (MQ=255)                                                                                                                                                                                   ttCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCg                                             >  2:203609/1‑251 (MQ=255)                                                                                                                                                                                     ccATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGAc                                           <  2:205220/251‑1 (MQ=255)                                                                                                                                                                                        tatCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATc                                        >  2:173307/1‑251 (MQ=255)                                                                                                                                                                                                     tCGTTTGATAACTCGGGTCCTGTTCCCCACGTACTACCCGTTGTTGATGGGCGAGCAGCGCGTGTTTATCTACAACGGCCCGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGTTGGCCTGTAAAAATCTCGAATCGTCTGTGAGCCGCGTCGCCTCTAATTCCAGCATGGTCAGTTCTTTTTCCGCCACCACGTTCCAGATCACCGCAGCGACATCGCCttgaagaat                            <  2:210538/250‑8 (MQ=255)                                                                                                                                                                                                      cGTTTAATAACTCGGTGCCTGTTCCCCACTTACTACCCGTGGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATg                           <  2:186609/250‑1 (MQ=255)                                                                                                                                                                                                                  cGGTTCCTGTTCCCCACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCATTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCTACATCGCCTTTAACAATGCGTTGTAAAct                <  2:223073/249‑1 (MQ=255)                                                                                                                                                                                                                             ccccACTTACTACCCGTTGTTTATGGTCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTAGTGATAaga    <  1:21059/250‑1 (MQ=255)                                                                                                                                                                                                                             ccccACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAaga    <  2:178436/250‑1 (MQ=255)                                                                                                                                                                                                                              cccACTTACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACTTCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAagaga  <  1:173858/251‑1 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     TTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTGACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGA  >  NC_000913/3511239‑3511710

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: ‑