Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,171,417	G→A	L850L (CTG→TTG)	mfd ←	transcription‑repair coupling factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,171,417	0	G	A	96.0%	73.9 / ‑4.7	25	L850L (CTG→TTG)	mfd	transcription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (8/16);  minor base C (1/0);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

CAAGACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAGTTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAG  >  NC_000913/1171174‑1171648                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        gaaGACGTTTTTTTGCATCTGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATAGGCCTTATGAGGCTAACGTCGGACGCGACCGCGTAACTGGTGCAGCTGCGCCATACCGAAGTGATCCGCGCGTTGAATGATAATAGGGTTGGCTGTCGGGATGTCGTTCCGGGTTTCGATAATGGTTGTACAATCCAGCACGTTTAAACGTTGATGATTGAAATCAGTCGTCACTCGTTCCAATtcgcg                                                                                                                                                                                                                                   <  2:166913/249‑1 (MQ=255)             gtgCATCGGGAGTCATCTCTTTTGGTTGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGAGGTCGAGCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCTTCTGc                                                                                                                                                                                                                      <  2:82322/251‑1 (MQ=255)               gCATCGGTAGTCCTCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCAAATTCGCGCTCGCGCATCTGccc                                                                                                                                                                                                                    >  1:191274/1‑251 (MQ=255)                cATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCg                                                                                                                                                                                                                   <  1:264471/251‑1 (MQ=255)                    ggTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGAc                                                                                                                                                                                                               >  1:83226/1‑251 (MQ=255)                                                  cAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCAGTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGtt                                                                                                                                                                                 <  2:199786/251‑1 (MQ=255)                                                               gCCTGATGTTGCGTACGTCCGACGCGACCGCGTACCTGGTGCAGCTGCGCAAGAGCGTCGTGATCCGCGCGTTCAATGATACTAGTGTTGGCTGTCGGGATGACGATCCCGGTTTCCATAATGGTTGTCCCAACCAGCACGTTGCAACGTAGATGATGGAATTCATTCATCACACGTACCACTTCGAGCTCGCGCCTCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCCGCCGTTc                                                                                                                                                                    >  2:255003/1‑251 (MQ=255)                                                                   gatgatGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTcggcg                                                                                                                                                                <  2:267319/251‑1 (MQ=255)                                                                       atgCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCtt                                                                                                                                                            <  1:102681/251‑1 (MQ=255)                                                                           gAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTg                                                                                                                                                        <  1:69912/251‑1 (MQ=255)                                                                            aaCGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGa                                                                                                                                                       <  2:187177/251‑1 (MQ=255)                                                                                                               gcCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTTGCTGCCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGCTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGt                                                                                                                                                                      <  2:32240/201‑1 (MQ=255)                                                                                                               gcCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGt                                                                                                                                                                      >  1:32240/1‑201 (MQ=255)                                                                                                                                   gcgTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCAc                                                                                                <  1:230207/251‑1 (MQ=255)                                                                                                                                     gTTCACTGATAATAGCGTTGGCTGTCGAACTGTCGATCCCGGTTTCGACAAAGGTTGTACAAACCAGCACGTTGAAACGTCGATGATGGAAATCATACATCACCCGTTCCAATTCGCGCCCGCGCATATGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTACTGAATGTTTTCCACATCATTGCAGAGCTAATAAACCAGTACTCCGCGCAACATTTCACGc                                                                                              >  2:23701/1‑251 (MQ=255)                                                                                                                                     gTTCAATGATCATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGc                                                                                              >  2:102485/1‑251 (MQ=255)                                                                                                                                            gATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATc                                                                                        <  1:209742/250‑1 (MQ=255)                                                                                                                                                        ggCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGacca                                                                           <  2:194544/251‑1 (MQ=255)                                                                                                                                                                        aTCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTTAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCCGCTTTCTGAATGGTGTCCACATCATTGTAGAGAGAATAAACCTGTCTTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGTACCATCATGCTGTCATAc                                                             <  1:255003/249‑1 (MQ=255)                                                                                                                                                                        aTCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTc                                                           >  2:102658/1‑251 (MQ=255)                                                                                                                                                                          cccGGTTGCGATAATGGTTGTACAAACCAGCCCGTTGAACCGGTGATGATGGAAATCATTCATCACCCGTTCCAATGCGCGCTCGCGCGTCTGCCCGTGAGCGATGGCGATCCGCGCTTCTGGCGCGATTTCTGCCAGCCGTTCGGCGGCTGTCTGAATGTTTTCCAGATCATTGTAGAGATAAGAAAGCTGTCCTCCGCTCAAAAATTCAGGCAGGATCGCCTCCCGGAGCACGATGCTTTCATACTCa                                                          <  2:46054/250‑1 (MQ=255)                                                                                                                                                                                    aTAATGGTTGTACAAACCAGCACGTTGAAACGGTGATGTTGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGAGCAACTGCCCGTGACCGATGGAGATAACCGCATATGGCACCAGTTCTGCCAGCCGTTAGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCACAATTTCCCGCAGGATCGCCTCCCGGACCCCCATGATGTCATACTCACGGACACAGGt                                               >  1:109392/1‑251 (MQ=255)                                                                                                                                                                                                    ccagcaagtagcaaCGTTCATGCTGGGAATCATTCATAACCAGATCCAATTCGCGCTCGAGCATCTGCCCGTGACAGATGGCGCTCCGCGCTTCCGGCCCCAGCTCTGCCAGCCGTTCCGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAACTTTCACGCAGGATCGCCTCCCGGACCACCCTGCTATCATACTCACGGACAAAGGTTTTACCTGCCAgacg                                >  2:182717/13‑250 (MQ=255)                                                                                                                                                                                                               aaaCGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGGGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGgcgggcgg                    <  2:83226/251‑1 (MQ=255)                                                                                                                                                                                                                                 ccattcatCACCCGTTCCAATTCGCGCTCGCGCGTCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTGCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAg  <  1:5555/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        CAAGACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAGTTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAG  >  NC_000913/1171174‑1171648

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