Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,053	T→G	T292T (ACT→ACG)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,639,053	0	T	G	90.3%	75.2 / ‑1.2	31	T292T (ACT→ACG)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (13/15);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

GTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTG‑GCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGT  >  NC_000913/4638824‑4639299                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        gTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTGTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGcc                                                                                                                                                                                                                                    <  1:153869/251‑1 (MQ=255)    aTAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCAtt                                                                                                                                                                                                                                 >  1:217844/1‑251 (MQ=255)                            aCTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCTTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATg                                                                                                                                                                                                         <  1:211446/251‑1 (MQ=255)                               ggCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCg                                                                                                                                                                                                      <  1:221287/251‑1 (MQ=255)                                             tttttAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTg                                                                                                                                                                                         <  1:351234/250‑1 (MQ=255)                                                tgaGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGaaa                                                                                                                                                                                      <  1:205781/248‑1 (MQ=255)                                                      gATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGGGTCGCCGTGACGACGCCACCCGATCAATACCAATTAACGGGCCGGGGGCCTACGGCCGCCATTCTGCTGATTGCACAGACTTTTATGGCGTTCTTTGTTTTTTACACGCTCAc                                                                                                                                                                               >  1:105185/1‑251 (MQ=255)                                                       gttttttgccacgcAAACGGGAAGGTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCGGGCAATGATACCGGCTGGAAAAACTTCCCCGCGTTTAGCGTCGCAGTGATGACGCCAGTCGATCAATACCAATTAACGGTCCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCAcc                                                                                                                                                                              <  2:109061/238‑1 (MQ=255)                                                       attttttgccAGCCAAACGGGAATTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGCTACCGGCTGGGAAAAATTCCCCGCGTTTAGCGTCGACGTGACGCCGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTAATCCATTCTGCTGATTGCACTCACTTTTACGGCGTTCTTTGTTTTTGAAACGCTCAcc                                                                                                                                                                              >  2:57378/9‑251 (MQ=255)                                                                     aaaCGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTaca                                                                                                                                                                <  1:309966/251‑1 (MQ=255)                                                                        cGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGCCTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACAccc                                                                                                                                                             >  1:376688/1‑251 (MQ=255)                                                                             aGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATCAACTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGACAACTTCCCCGCGTTTAGCGTCGCCGTGACGACGCCAGCCGGTCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGAGATCTTTGTTTTTGCAACGCTCACCGCGCAACGTTTACACCCAATGc                                                                                                                                                        >  1:172378/1‑251 (MQ=255)                                                                                     aGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTg                                                                                                                                                <  1:250878/251‑1 (MQ=255)                                                                                                   gcgcAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTGTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTc                                                                                                                                  <  1:131134/251‑1 (MQ=255)                                                                                                                 cagcTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGAc                                                                                                                                                                                                                                              >  1:335766/1‑129 (MQ=255)                                                                                                                 cagcTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGAGACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGAc                                                                                                                                                                                                                                              <  2:335766/129‑1 (MQ=255)                                                                                                                           gCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGctct                                                                                                           <  1:57378/250‑1 (MQ=255)                                                                                                                              aataatCTCGGTTAGCGTTTTGCTTCAGGCAATGATTCCGGCTGGGAAAACTTCCTCTCGTTTAGCGTCGCAGTGGCTACGCCAGCCGATCAATACCAATTTACGGACCGGGCGATTAAGTACGCCATTCTGCTGATTGCACTGACTTTTAGGGCTTTCTTTGTTGTTGAAACGTTCACCGCGCAACGTTTTCACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCt                                                                                                       <  1:356427/251‑1 (MQ=255)                                                                                                                              aataagCTCGGTGGGCGTTTTGCTTCGGGCAATGATACCGGCGGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCt                                                                                                       <  1:110793/251‑1 (MQ=255)                                                                                                                                     ccggggatctttttGATTCAGTCAATGAAACCGGCGGGGAAAACTTCCCCGTGTTTAGGGTCGCAGTGACGGCGTCTGCCGAACAATACCAATTAACTGACCGGGCTACTAAGGACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTGTTGAAACGCTCACCGCGCAACGTTTACCCCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTG‑GCGCtt                                                                                               <  2:217844/241‑1 (MQ=255)                                                                                                                                         tGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGGCTAAGTACGCCATTCTTCTGATTGCACTGACTTTTAAGGCGTTCTTTGTTTTTGAACCGCTCACCGCGAAACGTTTAAACCCAATGCAATCTTTGCTGGTGTGGCTTTCTTTTGTGATGTTTTATTTGCTCTTG‑TGTGCGCTTTCTg                                                                                           >  2:130521/1‑251 (MQ=255)                                                                                                                                                      tCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTG‑GCGCTTTCTGAACATATCGGttt                                                                              >  1:138792/1‑251 (MQ=255)                                                                                                                                                             aTGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTG‑GCGCTTTCTGAACATATCGGTTTTACCGTg                                                                       >  1:175734/1‑251 (MQ=255)                                                                                                                                                                  cccGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTG‑GCGCTTTCTGAACATATCGGTTTTACCGTGGCATg                                                                  <  1:377818/250‑1 (MQ=255)                                                                                                                                                                  aCCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTG‑GCGCTTTCTGAACATATCGGTTTTACCGTGGCATg                                                                  <  1:82269/251‑1 (MQ=255)                                                                                                                                                                                                     agtgacgacgCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTG‑GCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACg                               >  1:343827/5‑251 (MQ=255)                                                                                                                                                                                                         acgacgCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTG‑GCGCttta                                                                                             >  2:206642/1‑184 (MQ=255)                                                                                                                                                                                                         acgacgCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTG‑GCGCTTTc                                                                                             <  1:206642/185‑1 (MQ=255)                                                                                                                                                                                                             cgCCAGCCGAACAATACCAATTAACGGAACGGACGACTAAGTACGCCATTCTGCTGATTGCACTGAATTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTAATTGGTGATGTTTTAATTGCTCTTGCTG‑TCGCTTTCTGAACATAACGGTTTTACCGAGGCATGAAAACTCGCCAGATTGATTgggacgtttatgataggtatttt                        >  2:356645/1‑230 (MQ=255)                                                                                                                                                                                                                    cGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTG‑GCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCa                  >  1:93819/1‑249 (MQ=255)                                                                                                                                                                                                                                  tAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTG‑GCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGGATTGCAAgac  >  1:311555/1‑249 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        GTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTG‑GCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGT  >  NC_000913/4638824‑4639299

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