Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,331,890	C→T	L463L (CTC→CTT)	proP →	osmolyte:H(+) symporter ProP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,331,890	0	C	T	87.1%	66.1 / 3.2	32	L463L (CTC→CTT)	proP	osmolyte:H(+) symporter ProP
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (11/16);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

CCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATG  >  NC_000913/4331653‑4332134                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     ccGTTAAAGCGCGCTGGTGGCGGCATTTAATATTTCGGTGCTGGGTGCTGGTCTGACGCCAACGCTGGCGGCCGTGCTGGTCGAAAGTGCGCAGAATCTGATGATGCCTGTCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACTGCAAATCGT‑CCGTTGAAAGTTGCGACACCGGCGGCGTCAGATATACATGAAGCGAAGTAAATTCTTGTCGAGCATTACg                                                                                                                                                                                                                                         <  2:374237/251‑1 (MQ=255)   gTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGat                                                                                                                                                                                                                                       <  1:208830/251‑1 (MQ=255)   gTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGat                                                                                                                                                                                                                                       <  2:62911/251‑1 (MQ=255)         gcgcgcTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATc                                                                                                                                                                                                                                 <  1:353953/251‑1 (MQ=255)            cggcggcgggggcATTTAAGATCTCTGTGCTGGTTGCCGGTCTGACGCCACCGAGGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTAGTACCTGATGGTAGTGGCTGTGGGTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGGGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAAGTCTGGTCGAGCAATGCGATAATATCGAg                                                                                                                                                                                                                              <  2:235544/247‑1 (MQ=255)                        gagttAATATTTCGGTTCTGGTTGCCGTTCTGTCGCCAACGCTGGCGGTCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCATATATACAGGAAGCGAAGTAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCgat                                                                                                                                                                                                                  <  2:31994/248‑1 (MQ=255)                         aTTTAATTTTTCGGTGCTGGGTGCCGGTCTGTCGCCAACGCTGGCGGGCTGGCTGGTTGAAAGCTCGCAGAATCTGATGATGCCTGTCTATTACCTGATGGTAGTGGCGGTGGTTGGTGTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCTTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAATGAAATTCTTGTCGAGTATTACGATAATAACGAGCAGAAAATCgatg                                                                                                                                                                                                                 <  2:51908/251‑1 (MQ=255)                                atTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGa                                                                                                                                                                                                          <  1:165924/251‑1 (MQ=255)                                   tCGGTTCTGGTTGCCGGTCTGACTGCAAAGCTTGCGTCCTGGCTGGTCGAAAGCTCGCAGAATTAGATGAAGCCTTCCTATTACCGGATAGTAGTGGCGATGGTTGGTTTATTCACCGGCGTAACCATGAAAGAGATGGCAAATCGT‑CCATTGAAAGGTGCGACACCTGCCGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAACTCGATGATAttttcc                                                                                                                                                                                                        <  2:99422/250‑5 (MQ=255)                                              tgctgggcggaccccAACGTTGGCGGCCTGGCTGTGCGAAATCTCGTAGACTCTTATGATGCCTGCCTAATACCTGAGGGTATTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCCAATCGT‑CCGTTGAAAGGGGCGACACCGGCGGCGGCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACgagaatggc                                                                                                                                                                                             <  2:24979/237‑6 (MQ=255)                                                                   ggcggcCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTAc                                                                                                                                                                         <  1:158525/249‑1 (MQ=255)                                                                       gcCTGGCTGGTCGAAAGCTCGAAGAATCTGATGATGACTGCCTATTCCCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCTT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAACCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAACCGTACCCGCCt                                                                                                                                                                   <  1:8191/251‑1 (MQ=255)                                                                                      atgtCGCAGAATCTGATTAGGCGTGCATATTACTTGATGGTAGTGGCGGTGTTTGGTTTAATCACCGGGGTAACCATGAAAGAGACGGCAAAGCGT‑CCGTTGAAAGGTGCGACACCGGGGGCGTCAGATATACAGGAAGCAAAGAAAATTCTTGCCGAGCATGCCGAAAATAACGAGCAGAAAATCGATGATATGGTCCCCGAGTTTCCCTATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATcc                                                                                                                                                    <  2:336578/248‑1 (MQ=255)                                                                                                                                     ggtggtTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGCACGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTAAGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGGGCAGCAACATCCGCGAATTGATGAATAAACTGAAAAGGATGGACTGATGTGACGCTGTc                                                                                                     >  1:382534/1‑251 (MQ=255)                                                                                                                                       tggttggaTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGA‑CCGATGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAAGCGAAACGTACCCGCCTGGTGCAACAACATCCGCGAATTGATGAATAAGCTGAACCGGATTGCCTGATGGGCCGCTGTCtt                                                                                                   >  1:269870/1‑251 (MQ=255)                                                                                                                                             gttTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAgg                                                                                             >  1:324695/1‑251 (MQ=255)                                                                                                                                                 aaTCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGcc                                                                                                                                                                                             >  1:369271/1‑151 (MQ=255)                                                                                                                                                 aaTCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGcc                                                                                                                                                                                             <  2:369271/151‑1 (MQ=255)                                                                                                                                                       cGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCTAATTGAAc                                                                                   <  1:195862/251‑1 (MQ=255)                                                                                                                                                        ggCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAAct                                                                                  >  1:250610/1‑251 (MQ=255)                                                                                                                                                         gCGTAACCATGAAGGAGAGGGCAAA‑CGTACCGTTGAAAGGTGCGACACCGGCGGCGTCAGTTATACAGGAGGCGAAGGAAATTCTTGTCTAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAActc                                                                                 <  1:63987/251‑1 (MQ=255)                                                                                                                                                                     aaGAGAAGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCCCGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGTAACGGATGGCCTGATGTGACGCTGTCTTATAAGGACAATTGAACTCTTAAGGTTCACt                                                                     >  1:70153/1‑251 (MQ=255)                                                                                                                                                                               aaaccct‑CCGTGGAAAGGTGGGCCCCCGGTGGAGTCAGATATACAGGAAGCGTAGGAAGTTCTTATCGAGCATTACGATAATATCGAGCATAAAATCGAGGAGATTGGCCAGGAGATTGCCGATTTGCAGGCGAAACGTTCCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGTATGGCCTGATGTGACGCTGTCTTATCAGGCAAATTGAATTCTTAAGGTTCACTTAATCTCTGa                                                           <  2:78874/245‑1 (MQ=255)                                                                                                                                                                                aaTCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGCGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGAc                                                          >  1:323305/1‑251 (MQ=255)                                                                                                                                                                                       ccGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGCTATTGACCACGAGATTGCCGCTTTGCAGGCGGAACGTACCCGCCTGGTGCAGCAAAATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCa                                                     >  2:289859/1‑250 (MQ=255)                                                                                                                                                                                          ttGAACGGAGCGACACCGGCGGCGTCAGATAGACAGGAAGCGAAGGAAATTCTTGTAGCGCATTACGATAATATCGAGCAGAAAATCGCAGATATTGACAACGATATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGCATTGATGAATAAGTTGAAACGGCTGGACTGATGTGACGCTGTCATAACAGGACAATTGACCTCTTAAGGTTCACTTCATCTCTGACGCACAtcct                                                 >  2:374752/1‑248 (MQ=255)                                                                                                                                                                                          ttGAAAGGTGCGACACCGGCGGCGTCAGATATAAAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCCACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGAGCAtcct                                                 >  2:284345/1‑248 (MQ=255)                                                                                                                                                                                                         ccGGCGGCTTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAAc                                  <  1:319066/251‑1 (MQ=255)                                                                                                                                                                                                           ggcggcGTCAGATATACCGGAAGCGAAGGAAATTCTTGTCGAGCATGACGATAATATCGAGCAGAAAATCAGTGATATTGAACAGGAGATTGCCGATTTGCAGGCGAACGGTACCCACCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCCGTATTAGCAGGCCAATTGAACTCTTAAAGTGCACTTAGTCTCTGAAGCGCATACCAACCTCCAGGTTGACgg                                >  2:260032/1‑251 (MQ=255)                                                                                                                                                                                                            gcggcgTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACgga                               <  1:364811/251‑1 (MQ=255)                                                                                                                                                                                                                        ataCAGTAAGCGTAGGAAATTCTTGTCTAGCATTACGATAATATCTAGCAGAAAACCGATTTTATGGACCTCGAGATTGCCGATTTGCAGGCTAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGTATGGCCTTATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAAt                   <  1:138407/251‑1 (MQ=255)                                                                                                                                                                                                                                         aaTTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATg  >  1:175271/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     CCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGT‑CCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATG  >  NC_000913/4331653‑4332134

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