Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,790,389	G→T	13.9%	Q286K (CAA→AAA)	ydiU ←	UPF0061 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,790,389	0	G	T	13.9%	107.5 / 4.3	41	Q286K (CAA→AAA)	ydiU	UPF0061 family protein
Reads supporting (aligned to +/- strand):  ref base G (21/14);  new base T (4/2);  total (25/16)
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

GCCGCATCCGTTCTCCATAATGCGTCAACAAAACCTGCTGATAGCTGTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCA  >  NC_000913/1790254‑1790525                                                                                                                                        |                                                                                                                                         gCCGCATCCGTTCTCCATAATGCGTCAACAAAACCTGCTGATAGCTGTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTa                                                                                                                                      >  2:699224/1‑140 (MQ=255) gCCGCATCCGTTCTCCATAATGCGTCAACAAAACCTGCTGATAGCTGTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTa                                                                                                                                      >  1:422754/1‑140 (MQ=255)     cATCCGTTCTCCATAATGCGTCAACAAAACCTGCTGATAGCTGTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCaa                                                                                                                                  <  1:835345/140‑1 (MQ=255)      aTCCGTTCTCCATAATGCGTCAACAAAACCTGCTGATAGCTGTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAAGTCCACACCGCGACGGCAGGTTGATTATCaaa                                                                                                                                 >  2:491882/1‑140 (MQ=255)      aTCCGTTCTCCATAATGCGTCAACAAAACCTGCTGATAGCTGTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATCGCGACAATGTCTGCGCCAGACGCAGTAAATAACACAACTCGACGGCAGGTAGATTATCaaa                                                                                                                                 >  2:491836/1‑140 (MQ=255)            tctcCATAATGCGTCAACAAAACCTGCTGATTGCTGTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTa                                                                                                                           <  2:1088423/140‑1 (MQ=255)                  tAATGCGTCAACAAAACCTGCTGATAGCTGTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGccc                                                                                                                     >  1:907614/1‑140 (MQ=255)                      gCGTCAACAAAACCTGCTGATAGCTGTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTtga                                                                                                                 <  2:233517/140‑1 (MQ=255)                      gCGTCAACAAAACCTGCTGATAGCTGTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTtga                                                                                                                 <  2:233318/140‑1 (MQ=255)                                    tgctgATAGCTGTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGAtt                                                                                                   <  2:96590/140‑1 (MQ=255)                                               gTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAAc                                                                                        <  2:1230290/140‑1 (MQ=255)                                                     ggCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGtt                                                                                  >  2:253007/1‑140 (MQ=255)                                                          cATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAATGCGACGGCAGGTTTATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTaa                                                                             >  2:1016484/1‑140 (MQ=255)                                                          cATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAATGCGACGGCAGGTTTATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTaa                                                                             >  2:1015138/1‑140 (MQ=255)                                                          cATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAATGCGACGGCAGGTTTATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTaa                                                                             >  2:1015025/1‑140 (MQ=255)                                                          cATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAATGCGACGGCAGGTTTATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTaa                                                                             >  2:1160124/1‑140 (MQ=255)                                                          cATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTaa                                                                             >  1:1085210/1‑140 (MQ=255)                                                          cATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGAGTATCAAAGCTGTACCGCCCTTGATGATCCGAGTGAGTACATATAAAACCGGGTTCGGaa                                                                             >  2:742164/1‑140 (MQ=255)                                                                     tCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGaaa                                                                  >  1:1149857/1‑140 (MQ=255)                                                                      cTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGaaaa                                                                 >  2:166221/1‑140 (MQ=255)                                                                           gCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCaaa                                                            >  1:121622/1‑140 (MQ=255)                                                                                 aaTGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGccc                                                      <  2:1110055/140‑1 (MQ=255)                                                                                                    ccAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGcccc                                   >  2:116145/1‑140 (MQ=255)                                                                                                    ccAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGcccc                                   >  2:1178007/1‑140 (MQ=255)                                                                                                      aGACGCTGTAAATTCCACAACGCGACGGCAGTTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCcag                                 <  1:529885/140‑1 (MQ=255)                                                                                                      aGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCcag                                 <  1:688442/140‑1 (MQ=255)                                                                                                            tGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCTTCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGa                           >  2:501933/1‑140 (MQ=255)                                                                                                            tGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGa                           >  2:242095/1‑140 (MQ=255)                                                                                                            tGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGa                           >  2:444509/1‑140 (MQ=255)                                                                                                            tGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGa                           >  2:443982/1‑140 (MQ=255)                                                                                                                    ccACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGACAtgttgt                   <  1:228936/140‑1 (MQ=255)                                                                                                                    ccACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGACAtgttgt                   <  1:927316/140‑1 (MQ=255)                                                                                                                        aaCGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATTATCCGAGTTATTACAAATAAAACCGTGTTCGTAATCATCAAGAAAACAAAACGGCCCGTAATAAAGCGTCTGCCCCAGCAGCGACATGTTGTCgga               >  1:622839/1‑139 (MQ=255)                                                                                                                        aaCGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGTGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATAAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGt               <  1:124054/140‑1 (MQ=255)                                                                                                                        aaCGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCCAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGt               >  2:145354/1‑140 (MQ=255)                                                                                                                        aaCGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCCTAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGt               >  1:622944/1‑140 (MQ=255)                                                                                                                              aCGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATtcat         <  1:795925/140‑1 (MQ=255)                                                                                                                              aCGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATtcat         <  1:341083/140‑1 (MQ=255)                                                                                                                               cGGCAGCTTTATTATCAAACCTGTAACGCCCTTTATTAACACACTGATTACAAATAAAACCGCGTACGAAAACATCAAGAAAACCAAACGCCCCGTAATCAAGCGACAGCCCCAGCAGCGACATGTTGTCGGTATtcatc        <  2:1144771/140‑1 (MQ=255)                                                                                                                                    ggTTTATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCAcccc   <  1:307013/140‑1 (MQ=255)                                                                                                                                     gTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTACTCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCa  >  2:1165399/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                         GCCGCATCCGTTCTCCATAATGCGTCAACAAAACCTGCTGATAGCTGTCCAGGGCCTCATTCAGGGCATCTACGGCAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCA  >  NC_000913/1790254‑1790525

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑