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	97.1%	105.1 / NA	34	L463L (CTC→CTT)	proP	osmolyte:H(+) symporter ProP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/19);  minor base G (0/1);  total (14/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 = 1.00e+00

CGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTT  >  NC_000913/4331643‑4332131                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  cgtcGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAAGCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGTTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCttgtc                                                                                                                                                                                                                                                <  2:169534/248‑1 (MQ=255)    cGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAg                                                                                                                                                                                                                                             <  1:109432/251‑1 (MQ=255)    cGCATATCCGTTACAGCGCGCTGGCGGCGGCATGTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAAGCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAg                                                                                                                                                                                                                                             <  1:274609/251‑1 (MQ=255)           ccGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACg                                                                                                                                                                                                                                      <  1:256544/251‑1 (MQ=255)             gTTACAGCGCGCTGGCGGCGGCATTTAATATTGCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGat                                                                                                                                                                                                                                    <  1:367626/251‑1 (MQ=255)                aCAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGataat                                                                                                                                                                                                                                 <  2:95864/251‑1 (MQ=255)                      cgcTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCATCGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTGTTGCCTGATGGTAGTGGCGGTGGTTTGTTGAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGAGATACAGGAATGGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAg                                                                                                                                                                                                                           <  1:462294/251‑1 (MQ=255)                               cggcATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTTGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATc                                                                                                                                                                                                                  <  2:461142/251‑1 (MQ=255)                                                        tGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGcc                                                                                                                                                                                          <  1:342621/250‑1 (MQ=255)                                                          ccGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGAt                                                                                                                                                                                       <  1:29187/251‑1 (MQ=255)                                                                          gCTGGCGGCCTGGCTGGTCGAAAGCTTGCAGAATCTGATTATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACGGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGATCAGAAAAGCGATGATATTGACCACGAGATTGCCGATTTGCAGTCGAAACGTa                                                                                                                                                                       <  2:403291/251‑1 (MQ=255)                                                                                             gAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGTTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACa                                                                                                                                                    <  1:448748/251‑1 (MQ=255)                                                                                                                                       gTAGTGGCGGTGGTTGGTTAAATCACCGGCGTAACCATGAAAGAGACGGAAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTAAGATATACAGGAAGCGAAGGAATTTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTTCCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGAtgtg                                                                                                          <  1:107324/251‑1 (MQ=255)                                                                                                                                              cGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGt                                                                                                   >  1:4756/1‑251 (MQ=255)                                                                                                                                              cGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATAGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGt                                                                                                   >  1:120423/1‑251 (MQ=255)                                                                                                                                                      ggttTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCa                                                                                            <  1:397102/250‑1 (MQ=255)                                                                                                                                                           aaTCACTGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGACACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCa                                                                                       >  1:128042/1‑250 (MQ=255)                                                                                                                                                           aaTCACCGGCGTAACCATGAAAGAGACGGCAAAGCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATAGACAGGAAGCGAAGGCAATTCTTGTCGAGCATTACGATAATATCGGGCAGAAAATCGGTGATATTGACCACGAGAGTGCCGATTTGAAGGCGAAACGTACCCGACTGGAGCAGCAAAGTCCGCAAATTGATGAATAAGCTGGAAAGGAAGGCCTGATGTGACGCTGTCTTATCAGGCCa                                                                                       >  1:192451/1‑250 (MQ=255)                                                                                                                                                                ccGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGaa                                                                                 <  2:292142/251‑1 (MQ=255)                                                                                                                                                                ccGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTGGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGaa                                                                                 <  1:340649/251‑1 (MQ=255)                                                                                                                                                                       aaCCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGCTAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGATGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTa                                                                           >  2:428533/1‑250 (MQ=255)                                                                                                                                                                          cATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGt                                                                       >  2:321463/1‑251 (MQ=255)                                                                                                                                                                                  agaCGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTACGGTTCACTTa                                                                >  1:42490/1‑250 (MQ=255)                                                                                                                                                                                      ggCAAATCGTCCGTTGAAAGGAGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGCGCAGACAATCGATGATATTGACCACGAGATTGCCGATTTGCAGTCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAAtctc                                                           >  1:328042/1‑251 (MQ=255)                                                                                                                                                                                           aTCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAGATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCCGGCCAATTGAACTCTTAAGGTTCACTTACTCTCTGAcg                                                      >  1:37332/1‑251 (MQ=255)                                                                                                                                                                                                  gTTGAAGGGTGCGACACCGGCGGCGTCAGATATACAGGCAGCGAAGGAAATTCTTGTCGAGCCTTACGATAATATCGGGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGAGACGCTGTCTTATAAGGCTAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATAc                                               >  2:386975/1‑251 (MQ=255)                                                                                                                                                                                                  gTTGAAAGGGGCAACACCGGGGGCGTAAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATTTCGAGCAGAAAATCGATGATATTGACCACGAGTTTGCCGATTTGCCGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGGTCGCTTAATCTCTGACTCGCATAc                                               <  2:484869/251‑1 (MQ=255)                                                                                                                                                                                                                       gtcgTCAGATATACGGGAAGCGAAGTAAATTCTGGACGAGCATTGCGATAAGACCGAGAAGAAAATGGATGATTTTGACCCCGAGATTGCCGATTTGCAGGCGAACCGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACGTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACggagg                          <  2:244311/249‑1 (MQ=255)                                                                                                                                                                                                                           tCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGgaggag                        >  1:445067/1‑249 (MQ=255)                                                                                                                                                                                                                             aGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTg                    >  1:24964/1‑251 (MQ=255)                                                                                                                                                                                                                                   aCAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATg               <  1:168680/250‑1 (MQ=255)                                                                                                                                                                                                                                          gCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCg       <  1:326537/251‑1 (MQ=255)                                                                                                                                                                                                                                            gAAGGAAATTCTTGTCGAGCATTCAGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCGTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACggagg                          <  2:18961/230‑1 (MQ=255)                                                                                                                                                                                                                                            gAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACggagg                          >  1:18961/1‑230 (MQ=255)                                                                                                                                                                                                                                               ggAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGttt  >  1:8835/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  CGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTT  >  NC_000913/4331643‑4332131

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