Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,261	A→T	I362L (ATA→TTA)	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,261	0	A	T	96.3%	71.9 / ‑4.3	27	I362L (ATA→TTA)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (10/16);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

AGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGAC  >  NC_000913/4639020‑4639506                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      agtgacgacgCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACg                                                                                                                                                                                                                                              >  1:343827/5‑251 (MQ=255)                cGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCa                                                                                                                                                                                                                                 >  1:93819/1‑249 (MQ=255)                              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)                                                                                                                                                                                                                                            ggggCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGTCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGTTGTGATGTGGGTACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATAACGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGa   <  1:222128/251‑1 (MQ=255)                                                                                                                                                                                                                                             gggCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGAc  <  1:74858/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      AGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGAC  >  NC_000913/4639020‑4639506

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