Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	243,084	G→A	R74C (CGC→TGC)	fadE ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	243,084	0	G	A	93.3%	90.7 / NA	30	R74C (CGC→TGC)	fadE	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (7/21);  minor base C (2/0);  total (9/21)
Fisher's exact test for biased strand distribution p-value = 8.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCGGCAGATCCGCCAGCTCATGGGTGATCTGGAAATCATTCGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCGAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTAGGA  >  NC_000913/242845‑243328                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     cggcggCAGATCCGCCAGCTCATGGGTGATCTGGAAATCATTCGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGgcg                                                                                                                                                                                                                                           <  1:100286/251‑1 (MQ=255)     ggcAGATGCGACACCTCGTGGGTGATCTGGAAATAATTCGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAAAAACGATTGCTCTTCTGAGGTCAGGCGCGGCTGCGGCTAGGTATGCAGCTTTTTACAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCACGTGGTGCCCGCATCAATCGCTTCTTTCTCAGAGAGCGACACCGGCGGCCTCACCTTACAGAAAACCCCAAACACCGgcgaga                                                                                                                                                                                                                                        >  2:263061/1‑247 (MQ=255)       cAGATCCGCCAGCTCATGGGTGATCTGGAAATCATTCGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTGTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGaaa                                                                                                                                                                                                                                     <  2:375602/251‑1 (MQ=255)        aGATCCGCCAGCTCATGGGTGATCTGGAAATCATTCGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCCTCACCTTACGGAAACCGCCAAAAACCGGCGCGGaaa                                                                                                                                                                                                                                     >  2:227939/1‑250 (MQ=255)          aTCCGCCAGCTCATGGGTGTTCTGGAAAGCATTCGCCATCCGGCGGGCTTCTGCTACCGGGCCGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGAGAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCGCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAatca                                                                                                                                                                                                                                  <  2:357014/251‑1 (MQ=255)          aTCCGCCAGCTCATGGGTGATCTGGAAATCATTCGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAatca                                                                                                                                                                                                                                  <  2:357033/251‑1 (MQ=255)            ccGCCAGCTCATGGGTGATCTGGAAATCATTCGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAatcatc                                                                                                                                                                                                                                <  1:89532/251‑1 (MQ=255)                       tGGGTGATCTGGAACTCATTCGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGGAACCGCAAAATACCGGCGGGGAAATCATCGACTTACGcg                                                                                                                                                                                                                      >  1:52505/1‑249 (MQ=255)                             aTCTGGAAATCATTCGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGgcgc                                                                                                                                                                                                               >  1:413466/1‑251 (MQ=255)                                 ggAAATCATTCGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCa                                                                                                                                                                                                              <  2:404497/248‑1 (MQ=255)                                   aaaTCATTCGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGtt                                                                                                                                                                                                         <  1:171999/251‑1 (MQ=255)                                        atgcGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAGAAACGCTTGCTCGTCGGCGGTCAGGCGCGGCGGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGATCAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAAGCGCTTCTGTCTCAGTGCGCGACATCGGCGGCATCTACGTACGGAAAGCGCAAAATACCGGCACGGAAATCAGCGACTTACGCCTAGGTGCAGAGTGAAATg                                                                                                                                                                                                    <  2:165527/248‑1 (MQ=255)                                                                    gggCCGTCGAGAAACGCTTGCTCGTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCTGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAg                                                                                                                                                                        <  2:396581/251‑1 (MQ=255)                                                                       ccGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTAc                                                                                                                                                                     >  1:325441/1‑251 (MQ=255)                                                                           cGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCAc                                                                                                                                                                 <  1:5188/251‑1 (MQ=255)                                                                                aaCGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCCATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATCGGCTCAAAGTTAAATGGCACTAGGATAATGTCCAGAGGCACCAGTACCCACGCCgc                                                                                                                                                            >  1:12517/1‑250 (MQ=255)                                                                                                 ggTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGc                                                                                                                                           >  1:338188/1‑251 (MQ=255)                                                                                                     aGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAg                                                                                                                                       <  1:375771/251‑1 (MQ=255)                                                                                                                                  tttCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCCCCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGACGAGGGCGGCTGTCCAGGCGAGCAAAATCAGa                                                                                                           <  2:90198/250‑1 (MQ=255)                                                                                                                                  tttCCAGTCCGGCTTGCCCTGGAAAAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGAc                                                                                                          >  2:397008/1‑251 (MQ=255)                                                                                                                                                        aaCAAGTCTCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCAcgcg                                                                                    <  2:121488/251‑1 (MQ=255)                                                                                                                                                        aaCAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGTGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCAcgcg                                                                                    <  2:110351/251‑1 (MQ=255)                                                                                                                                                                           aGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCg                                                                 <  2:312643/251‑1 (MQ=255)                                                                                                                                                                                        aTCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCCCCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCg                                                    <  1:226815/251‑1 (MQ=255)                                                                                                                                                                                           aaTCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAAGACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAg                                                 <  2:358538/251‑1 (MQ=255)                                                                                                                                                                                                  tctttcgCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGaa                                           <  1:247699/250‑1 (MQ=255)                                                                                                                                                                                                             gcgcGACATCGGCGGCATCACCTTACGGAAACCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAatc                                <  1:227939/250‑1 (MQ=255)                                                                                                                                                                                                                                         aaaCCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTAccc                                                                                                                                                                   >  2:157685/1‑91 (MQ=255)                                                                                                                                                                                                                                         aaaCCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTAccc                                                                                                                                                                   <  1:157685/91‑1 (MQ=255)                                                                                                                                                                                                                                          aaCCGCAAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTagga  <  1:362460/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     CGGCGGCAGATCCGCCAGCTCATGGGTGATCTGGAAATCATTCGCCATCCGGCAGGCTTCTTCTACCGGGCCGTCGAGAAACGCTTGCTCTTCGGCGGTCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACCTTACGGAAACCGCGAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTAGGA  >  NC_000913/242845‑243328

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