Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,046,782	G→T	10.0%	intergenic (‑180/‑184)	yihG ← / → polA	inner membrane protein, inner membrane acyltransferase/5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,046,782	0	G	T	10.0%	113.6 / 3.4	37	intergenic (‑180/‑184)	yihG/polA	inner membrane protein, inner membrane acyltransferase/5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease
Reads supporting (aligned to +/- strand):  ref base G (20/13);  new base T (2/2);  total (22/15)
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.45e-01

AGCAAAAGGGTTATCGCAGCGAGTATTCTCGTCATAATGAATTTATTCAAAAGGTTAGCCATAGACAGTGCGCTAAAGAAGCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATCCTTTACCCATAATCTGATAAACAGGCACGGACATTATGGTTCAGATCCCCCAAAATCCACTTATCCTTGTAGATGGTTCATCTTATCTTTA  >  NC_000913/4046541‑4047021                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                aGCAAAAGGGTTATCGCAGCGAGTATTCTCGTCATAATGAATTTATTCAAAAGGTTAGCCATAGACAGTGCGCTAAAAAAGCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATACAGGGAGAACCTTAAGGAGAAAAATAATTCATAGCTATCCACTTAAGCAAAAATCCCATTATCACAATTATTAGGGATGGATTTATTTTAAACTGCATGAAAACAAAAGACAAACATCATGCTTTAAAAAGCATGATAATAACATAAAAGGCATGTa                                                                                                                                                                                                                >  4:437245/1‑275 (MQ=255)    aaaaGGGTTATCGCAGCTAGTATTCTCGTCATAATGAATTTATTCAAAAGGTTAGCCATAGACAGTGCGCTAAAGAAGCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGCTCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCCGAAAAAAGCATGATAATAAATTAAAAGCATTTTAaat                                                                                                                                                                                                             >  4:340866/1‑275 (MQ=255)                                          tttATTCAAAAGGTTAGCCATAGACAGTGCGCTAAAGAAGCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTgcgag                                                                                                                                                                       <  2:340870/275‑1 (MQ=255)                                                                     tGCGCTAAAGAAGCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAGGATCTACAAAAGAtttc                                                                                                                                            >  2:398771/1‑274 (MQ=255)                                                                                                            cGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTTTTTACAGCATGatg                                                                                                                                                                                                                                   <  6:356925/149‑2 (MQ=255)                                                                                                                ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTTTTTACAGCATGat                                                                                                                                                                                                                                    >  8:168846/1‑144 (MQ=255)                                                                                                                ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTTTTTACAGCATGat                                                                                                                                                                                                                                    <  7:168846/144‑1 (MQ=255)                                                                                                                     tAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTaa                                                                                                                                                                                                                          >  8:125058/1‑149 (MQ=255)                                                                                                                         ggATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGc                                                                                                                                                                                                                      >  8:210568/1‑149 (MQ=255)                                                                                                                          gatgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCg                                                                                                                                                                                                                     <  7:210568/149‑1 (MQ=255)                                                                                                                           atgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGa                                                                                                                                                                                                                    <  7:125058/149‑1 (MQ=255)                                                                                                                           atgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGa                                                                                                                                                                                                                    <  5:361652/149‑1 (MQ=255)                                                                                                                             gatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATCCTTTACCCa                                                                                    >  1:18492/1‑275 (MQ=255)                                                                                                                              atCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAg                                                                                                                                                                                                                                          <  5:286088/124‑1 (MQ=255)                                                                                                                              atCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAg                                                                                                                                                                                                                                          >  6:286088/1‑124 (MQ=255)                                                                                                                                        agaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTa                                                                                                                                                                                                       <  6:147819/149‑1 (MQ=255)                                                                                                                                                      tatctatcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATc                                                                                                                                                                                         <  5:139019/149‑1 (MQ=255)                                                                                                                                                          tatcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTaaaa                                                                                                                                                                                                                        <  5:168803/114‑1 (MQ=255)                                                                                                                                                          tatcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTaaaa                                                                                                                                                                                                                        >  6:168803/1‑114 (MQ=255)                                                                                                                                                          tatcCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCcaca                                                                                                                                                                                     >  6:174916/1‑149 (MQ=255)                                                                                                                                                              cacaTTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGAc                                                                                                                                                                                 >  7:325962/1‑149 (MQ=255)                                                                                                                                                                  ttAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGAtt                                                                                                                                                                             >  6:143661/1‑149 (MQ=255)                                                                                                                                                                      aaaaaaTCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAGGTTATCCACATGACGATTTgcg                                                                                                                                                                         >  7:410178/1‑149 (MQ=255)                                                                                                                                                                                 ttATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCaga                                                                                                                                                              >  5:125710/1‑149 (MQ=255)                                                                                                                                                                                     ctcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCTAGAAGAt                                                                                                                                                          >  7:291514/1‑149 (MQ=255)                                                                                                                                                                                      tcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATc                                                                                                                                                         <  8:410178/149‑1 (MQ=255)                                                                                                                                                                                      tcAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATc                                                                                                                                                         >  5:182092/1‑149 (MQ=255)                                                                                                                                                                                                        tttatttTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGa                                                                                                                                       >  6:29485/1‑149 (MQ=255)                                                                                                                                                                                                         ttatttTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGaa                                                                                                                                      >  5:236714/1‑149 (MQ=255)                                                                                                                                                                                                            tttttAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAg                                                                                                                                   >  8:369844/1‑149 (MQ=255)                                                                                                                                                                                                                aaaaaaaataaaaaaCAAAAANNAACATNATGGTGTAAAAAGCACGATAAAAAATTAAAAGCTAGGTAAATAAATTAGACACCAAGTTATCCACATGACGATTTGCGAGCGATCCAGGAGATCTACAAAAGATTTTCACGAAAAGCAGTGAAAAACTCATCTGTTCATCCTGTCTGTGGCATCCTTTACCCATAATCTGATAAACAGGCACGGACATTATGGTTCAGATCCCCCAAAATCCACTTATCCTTGTAGATGGTTCATCTTATCTTTa  <  3:18492/264‑1 (MQ=255)                                                                                                                                                                                                                      cATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCa                                                                                                                     >  1:83403/1‑153 (MQ=255)                                                                                                                                                                                                                          aaaaaCAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGtt                                                                                                                 <  3:83403/153‑1 (MQ=255)                                                                                                                                                                                                                                  aGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCa                                                                                                             <  5:29485/149‑1 (MQ=255)                                                                                                                                                                                                                                   gACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCAt                                                                                                            >  5:263119/1‑149 (MQ=255)                                                                                                                                                                                                                                      aaaCATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCct                                                                                                         >  8:371919/1‑149 (MQ=255)                                                                                                                                                                                                                                            catGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCctgtctgt                                                                                                   <  5:174916/149‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                AGCAAAAGGGTTATCGCAGCGAGTATTCTCGTCATAATGAATTTATTCAAAAGGTTAGCCATAGACAGTGCGCTAAAGAAGCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGATGTAAATAATTTATGCACAAAGTTATCCACATGACGATTTGCGAGCGATCCAGAAGATCTACAAAAGATTTTCACGAAAAGCGGTGAAAAACTCATGTTTTCATCCTGTCTGTGGCATCCTTTACCCATAATCTGATAAACAGGCACGGACATTATGGTTCAGATCCCCCAAAATCCACTTATCCTTGTAGATGGTTCATCTTATCTTTA  >  NC_000913/4046541‑4047021

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑