Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,883	C→T	L236M (CTA→TTG)	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,638,883	0	C	T	95.5%	62.2 / ‑6.0	23	L236M (CTA→TTG)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (9/12);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTAC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCAC  >  NC_000913/4638650‑4639127                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      cgcCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGaa                                                                                                                                                                                                                                     <  1:351465/251‑1 (MQ=255)      gAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAAACAGCCACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACggggaggt                                                                                                                                                                                                                                  >  1:243596/1‑244 (MQ=255)           tAACGGAACGGCGCTGCCCATTGACCCCTTCCCCGGGTTAGAGCACGGCGGGCAGAGCGTGCATATCCCTTTACCCGAAGGGGACTGGCTGAAGCAGAACCTGAATGTGAATATGGACCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGTCGTAATAGCGCAATGACCTTAACCAGCAACTGGCCGAATCCCAGTTTTTTAGGTGATTTTGTGC‑CAGCCAACCGGTAAGTTAGCgcg                                                                                                                                                                                                                            <  2:284779/250‑3 (MQ=255)               ggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGGGGTACCGGCGATCTTTCTGTGGATCCTGCCGGGCGTAATAGCGAAATTGCCTTAACCAGCAACTGGCCGCATCCCAGGTTTTTAGGTGATTTTTTGC‑CCGCCAAACGGGAAATTaacggggaggg                                                                                                                                                                                                                       >  2:309966/1‑241 (MQ=255)                     gcgcTGACCATTGAACCCGGCTCCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGTGACTGGCGGAAGCAGAACCGGAATCTGAATATTGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTTATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCa                                                                                                                                                                                                                 <  2:243596/251‑1 (MQ=255)                          gACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGgcgc                                                                                                                                                                                                            <  1:93614/251‑1 (MQ=255)                           cccATTGAACCCGGCACCGGGTTAGAGCAAGACGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAATCTGAATATGGCCCTTAAGTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCa                                                                                                                                                                                                           <  2:328199/250‑1 (MQ=255)                                  aaCCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCa                                                                                                                                                                                                     >  1:109061/1‑250 (MQ=255)                                            cGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTGTCCCGAAGGTGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCATGCGCAATGGCAAAGCAGCTGGt                                                                                                                                                                                          <  2:258279/251‑1 (MQ=255)                                                  agagCAAGGCGGGCAGGGCGGGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCta                                                                                                                                                                                    <  1:299413/251‑1 (MQ=255)                                                      cAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCtaataa                                                                                                                                                                                <  1:59510/251‑1 (MQ=255)                                                                                  tACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGAACCGGCGATCTTTCTGTGGGGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTTGGTGATTTTTTGCACA‑CCAAACGGGAAGGTAGCGCGTCAGGTTTTCCGGCGCAATGGCAAAGCAGCTGGTTTTATAATAATCTCGGTGAGCGTTTTGCTTCAgg                                                                                                                                                        >  2:126504/1‑247 (MQ=255)                                                                                                            aaCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTccccgc                                                                                                                          >  1:62440/1‑248 (MQ=255)                                                                                                                                                    gCGATATTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCa                                                                                   >  1:128722/1‑250 (MQ=255)                                                                                                                                                           ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTaa                                                                              <  1:126504/248‑1 (MQ=255)                                                                                                                                                                               gTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTGTTAGGTGATTTTCTAC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGcc                                                       <  1:153869/251‑1 (MQ=255)                                                                                                                                                                                  aTAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCAtt                                                    >  1:217844/1‑251 (MQ=255)                                                                                                                                                                                                  aaCCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAAcg                                                                            >  1:58947/1‑210 (MQ=255)                                                                                                                                                                                                          aCTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCTTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATg                            <  1:211446/251‑1 (MQ=255)                                                                                                                                                                                                             ggCCGCATCCCAGTTTTTTAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCg                         <  1:221287/251‑1 (MQ=255)                                                                                                                                                                                                                           tttttAGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTg            <  1:351234/250‑1 (MQ=255)                                                                                                                                                                                                                              tgaGGTGATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGaaa         <  1:205781/248‑1 (MQ=255)                                                                                                                                                                                                                                    gATTTTTTGC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGGGTCGCCGTGACGACGCCACCCGATCAATACCAATTAACGGGCCGGGGGCCTACGGCCGCCATTCTGCTGATTGCACAGACTTTTATGGCGTTCTTTGTTTTTTACACGCTCAc  >  1:105185/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      CGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTAC‑CAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCAC  >  NC_000913/4638650‑4639127

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