Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,284	A→G	Q369Q (CAA→CAG)	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,284	0	A	G	92.6%	72.1 / ‑2.0	27	Q369Q (CAA→CAG)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (9/16);  total (10/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 = 8.62e-01

TAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTAC  >  NC_000913/4639049‑4639516                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                         tAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGGATTGCAAgac                                                                                                                                                                                                                           >  1:311555/1‑249 (MQ=255)   acggACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAggctg                                                                                                                                                                                                                         >  1:84772/4‑248 (MQ=255)    cggACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTgg                                                                                                                                                                                                                        <  1:72277/249‑1 (MQ=255)         cGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAggccggcgca                                                                                                                                                                                                                    >  2:74858/1‑242 (MQ=255)                                   gATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCtgt                                                                                                                                                                                         <  1:130521/251‑1 (MQ=255)                                                   aTGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGgtgg                                                                                                                                                                         >  1:75335/1‑250 (MQ=255)                                                     ggCGTTCTTTTTTTTGGAAACGTTCACCGCGCATCGTTTACTCCCAATGCAATATTTGCTGGTGGGGCTTTCATTTGTGATGTTTTATTTGCTCTTGATGGCGCTTTCTGAACATATCGGTGTTACCGTGGCATGTATAATCGCCAGTCTGATTGGGGCGTGAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGGTTACCCTCGCGCTGTTGTTGCTGGATGGtgtg                                                                                                                                                                        <  2:235578/250‑1 (MQ=255)                                                              tGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATTTTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGa                                                                                                                                                               <  1:54217/250‑1 (MQ=255)                                                                  tttGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCt                                                                                                                                                          >  1:168824/1‑251 (MQ=255)                                                                               ccGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGGAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGCTGGTGTGATGTGGGGACTGCTCAACTCTGCCGAt                                                                                                                                             >  1:308999/1‑251 (MQ=255)                                                                                       cGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTTCTGGATGGTGTGATGTGGGGACTGCTCAATTCTGCCGATAGCGCGCt                                                                                                                                     <  1:104411/251‑1 (MQ=255)                                                                                        gTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGGTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGTATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCtg                                                                                                                                    <  1:44838/251‑1 (MQ=255)                                                                                           tACACCCAATGCAATATTTGCTGGTGGGTCTTTCATTGGGTATGTGTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGTATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCtgtt                                                                                                                                  <  2:175734/250‑1 (MQ=255)                                                                                             cacCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTgttgtt                                                                                                                               <  1:93950/251‑1 (MQ=255)                                                                                                   aTGCAATCTTTGCTGGTGGGGCTTTCATTGGCGAAGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTTAAAGGCTGGCGCAACAGCCTGTTGTTTACCCTCGCGCTGTTGTTGCTGGGTGGTGTGATGTGGGGACTGCTCAACTCTGCCGATTACGGGCTGTTGTTgggta                                                                                                                          >  1:337906/1‑248 (MQ=255)                                                                                                      cAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAg                                                                                                                      <  1:7001/251‑1 (MQ=255)                                                                                                                                  tGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGGATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGAGGGTGTGGGGTGGGGGCTGTACAACTCTGTCGATAGCGCGCTGTTGTTGGGGACAAGTGTGCtaggggggggtctgggcggactg                                                                                          >  2:299231/1‑229 (MQ=255)                                                                                                                                                         gcgcTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGTGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTtgt                                                                                  <  1:41039/236‑1 (MQ=255)                                                                                                                                                         gcgcTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGa                                                                   <  1:218435/251‑1 (MQ=255)                                                                                                                                                                            gTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGAAGGCGGTATTGAAAGGCTTGCGCGACAGCACGTTGTTTAACCTCGCGCTGTTGTTGCTGGATGGGGTGATGTGGGGACTGCTCAACTCTGCCGATAACGCGCGGGTTTTGGGAACCAGTGTGCTGGGGGTGGCGCTGGCCGGCATGGTGTTTTGGACTCGTAATATCGACTGGTATGGGTTTTCACTg                                                >  1:103121/1‑251 (MQ=255)                                                                                                                                                                                tacgggtgCATGTATAAACGCTAGTCTGATAGGGGCGTTAATTAACGGTATTCATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGTTTGTTGTTTACCCGCGCGCGTTTGTTGCTGGATGGTTTGATGTGGGGATTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGTTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACGCGTAATATCGACTTGTATGCGTTTTTACTGCCg                                             <  2:343827/243‑1 (MQ=255)                                                                                                                                                                                                   gCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAaaga                         <  1:111018/251‑1 (MQ=255)                                                                                                                                                                                                           taTTGTGGCGTTAATGAACGGTATTTATTTGTAGGCGGTATGGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGTATGGTGTGATGTGGTGTCTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATTATGGTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTATAGAAGTTa                    <  1:42803/247‑1 (MQ=255)                                                                                                                                                                                                            agttgggcattAATGACCGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTCCAACAGCATGTTGTTTACCCTCGCGCTGTGGTTGCTGGATGGTGTGATGTGGGGATTGCTCATCTCCTCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTTGCGCTGGCCGGCATGATGTTTGTGACCCGTAATAACGACTGTTATGCGTTTTCCCTGCCGAAAATGAAAGCCAGTAAAGAAGATACaa                 <  2:93819/240‑1 (MQ=255)                                                                                                                                                                                                               ggggCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGTCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGTTGTGATGTGGGTACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATAACGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGa             <  1:222128/251‑1 (MQ=255)                                                                                                                                                                                                                gggCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGAc            <  1:74858/251‑1 (MQ=255)                                                                                                                                                                                                                          gcaCGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAAAAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGAGATGTGGGGACTGCACACCTCTGCCGATAGCGCGCTGTTGTTGGGAAACAGTGTGCTGGGGGTGGCGATGGCCGGCATGATGATTGGGACCCGTAATATCGACGGGTGTGCGTTTTAACTTCCGAACTCGAAAGCCAGAAAAGAAGGTAAAACGGaaggttctgtgc  >  2:101945/3‑240 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                         TAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTAC  >  NC_000913/4639049‑4639516

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