Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,568,442	G→A	L128F (CTT→TTT)	digH ←	divisome localized glycosyl hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,568,442	0	G	A	96.2%	81.6 / NA	26	L128F (CTT→TTT)	digH	divisome localized glycosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/16);  minor base C (1/0);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTACCGGGCTTCGTATTAACCGATACGCGATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGA  >  NC_000913/1568202‑1568667                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  ccGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTACCGGGCTTCGTATTAACCGATACGCGATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCAc                                                                                                                                                                                                                         <  2:139683/251‑1 (MQ=255)     gTCGAGGACAAAGCGCTCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGAGAGAGTGCTATTCAGTTCCCTGATAGCAACGCGCTTAGTCTTAACCGATACGCGATAGGGGTTAACCCAGGCGTGTACTTTCATACCACGCTTGTAGGCTTCGGCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTAATAAAATCGGACCACGGc                                                                                                                                                                                                                      >  1:235407/1‑250 (MQ=255)            aCAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTACCGGGCTTCGTATTAACCGATACGCGATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAAtttt                                                                                                                                                                                                              <  1:289400/251‑1 (MQ=255)                                            gcgGTGTTGCACATAGACGCTCGCTGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTACCGGGCTTCGTATTAACCGATACGCGATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGa                                                                                                                                                                              <  2:249449/251‑1 (MQ=255)                                                  ttGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTACCGGGCTTCGTATTAACCGATACGCGATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTgg                                                                                                                                                                         >  1:183921/1‑250 (MQ=255)                                                      acaTAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTACCGGGCTTCGTATTAACCGATACGCGATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAACATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGaaaa                                                                                                                                                                     >  1:53981/1‑250 (MQ=255)                                                                     ggTTGCTGAGACAGAGTGCTACTCAGTTTCCTGATAGTACCGGGCTTCGAATTAACCGATACGCGATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTAGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATAGGACCACGGCAAACTTTTCGATGGCCACAGGGCGGTACCGTCCGGATTGACCTGGAAAAAGACCGTGTTTATGcc                                                                                                                                                     >  2:175077/1‑251 (MQ=255)                                                                              gaCAGAGTGCTATTCAGTTCCCTGATAGTACCGGGCTTCGTATTAACCGATACGCGATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTg                                                                                                                                            <  2:128132/251‑1 (MQ=255)                                                                                         aTTCAGTTCCCTGATAGTACCGAGCTTCGTATTAACCGATACGCGATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATc                                                                                                                                   <  2:390390/249‑1 (MQ=255)                                                                                                          tACCGGGCTTCGTATTAACCGATACGCGATAGGGGTTAAACCAGGCGTGGACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAACGGTTCGTCACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGTACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTTCACTGGAAAAAAACCGTGTTTATTCCGAGAAGTTGCCGATGATCCAGTTTGTCGGTCATTGcc                                                                                                                >  2:123972/1‑251 (MQ=255)                                                                                                            cggggCTTCGTATTAAACGATACGCGATAGGGGGTAAAACATGCGGGTACTTTCATCCCTCGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCCCCAATCTTACCGGGCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGTGCGGTGCCGTCCGGTTTGACCTGAAAAAAGACCGGGTTTATGCCGATACGTTGCAGATGATCCAGTTTGTTGACCATCGCCtg                                                                                                              <  2:392035/249‑1 (MQ=255)                                                                                                                                  cgcgATATGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtg                                                                                        <  1:332332/251‑1 (MQ=255)                                                                                                                                    cgATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtggt                                                                                      >  1:317804/1‑251 (MQ=255)                                                                                                                                    cgATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtggt                                                                                      <  2:119857/251‑1 (MQ=255)                                                                                                                                     gATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtggtg                                                                                     <  2:339649/251‑1 (MQ=255)                                                                                                                                                      gCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCTAGCATAAATTGCAGCGGATCGTAGCCCGGATTATCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAAATTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGTGACGTTGCAGATGATCTATTTTGTCGATCATCGCCGATTGTTGTACACGGGGCCGGGTGGTCGTGTTTCTAATTTTAAc                                                                    <  2:255887/251‑1 (MQ=255)                                                                                                                                                      gCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAAc                                                                    >  2:10911/1‑251 (MQ=255)                                                                                                                                                                    cccACGCTTGTGGGCGTCGTCGGGCATAAATTGCAGCGGATCGTAACCCGGATTTGCGCCAATCTTACCGGTCATAAAATCGGACCACGGCAAAAGTTTCGATGGCCACAGGGCGGTACCGGCCGGCTTGACCTGGAAAAATACCGTGTTTATGCCGAGACGTTGCAGATGTTCCAGTTTGTCGATCATCGTCTGTGGTGGTACACGGGCCCGGCTGGTGGGGTTACTATTTTTGACCGAgcagtccggtg                                                      <  2:396189/251‑11 (MQ=255)                                                                                                                                                                        cGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGAGGcc                                                   >  1:274889/1‑250 (MQ=255)                                                                                                                                                                         gCTTGTGGGCTGCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGACCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAg                                                 <  2:53981/251‑1 (MQ=255)                                                                                                                                                                            tgtgGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTc                                               <  1:201398/250‑1 (MQ=255)                                                                                                                                                                                        cgAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAATACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCTGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAAcc                                  <  1:362670/251‑1 (MQ=255)                                                                                                                                                                                               aaaTTGCAGCAAATCGTAACCCGGTTTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGCGCAACCGTGgcc                            >  2:234776/1‑250 (MQ=255)                                                                                                                                                                                                           aTCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCAAAGGGCGGTACCGTCCGGCTTGACTTGGAAAAAGACCGTGGTTATGCCGAGACGTTGCAGATGATCAAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGGCCGGCTGGTGGGGTTACTAATGTCAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCAc               <  2:317804/251‑1 (MQ=255)                                                                                                                                                                                                                     ggATTTTCACCAATCTTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTTGCCAGCCAGATGCCACGCATCGGTtg     >  2:238752/1‑251 (MQ=255)                                                                                                                                                                                                                         tttCACCAATCGTACCGGTCATAAAATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGTTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGa  <  1:372458/250‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  CCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTACCGGGCTTCGTATTAACCGATACGCGATAGGGGTTAAACCAGGCGTGTACTTTCATCCCACGCTTGTGGGCTTCGTCGAGCATAAATTGCAGCGGATCGTAACCCGGATTTTCACCAATCTTACCGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGA  >  NC_000913/1568202‑1568667

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