Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,141,865	G→A	P835S (CCG→TCG)	rne ←	ribonuclease E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,141,865	0	G	A	97.0%	89.0 / NA	33	P835S (CCG→TCG)	rne	ribonuclease E
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (11/21);  minor base C (1/0);  total (12/21)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTCACCAGTTCTGTGCGTGGAGCTGGCGCTTCCTGAACAATTG  >  NC_000913/1141623‑1142110                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      aGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACAAGGACTTCTTCCTGTTAGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGGCAGTTCCGAAGACgcgc                                                                                                                                                                                                                                               >  2:197301/1‑251 (MQ=255)   gATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGATCGATCGCGGCGGCGACAGGGACCTCAGTCCACATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATTTTGCGGACGTACAATTGGATAGTGGATCCAGACTTTCCCACAGGCCAGTTCCGAAGACGCGCAc                                                                                                                                                                                                                                             >  1:148832/1‑251 (MQ=255)                                   tCGGCAACCTGAACGTGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTg                                                                                                                                                                                                              <  2:77984/250‑1 (MQ=255)                                   tCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCCTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCTAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTg                                                                                                                                                                                                              <  1:233521/250‑1 (MQ=255)                                             ggaCGGGGGCTTCTACGACGCCGGCCTCTTCGTCAACAACGCGCGCGCTAACAACCGGTGCGATAGCTGCGGCGACAGGTACCTCAGTCACCATCGGCTGCCCGTGTACTTCTTCCTGTTCGCGCAGCTCTTCAACCTGTACATCTTGCGGACTTACAATGGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCTCACGCTACGGTCAACGGCATTGGCGACTGTTTTGGATAACGCTCg                                                                                                                                                                                                   <  2:108216/249‑1 (MQ=255)                                                                cccGGCCTCTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACGTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCGTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCgacgacgac                                                                                                                                                                                <  2:24337/250‑1 (MQ=255)                                                                cccGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGGTCGAGAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTGCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAAGTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTTGCGACTGGGTTGGATAACGCTCGTCACGATAGCgacgacgac                                                                                                                                                                                <  2:102671/250‑1 (MQ=255)                                                                    gCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCGCGTCACGATAGCgacgacgac                                                                                                                                                                                <  2:158980/247‑1 (MQ=255)                                                                    gCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGAACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCgccgacgac                                                                                                                                                                                >  1:158980/1‑247 (MQ=255)                                                                      cACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCCGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGcgacgacggacacgac                                                                                                                                                                          >  2:2873/1‑243 (MQ=255)                                                                              caacaaCTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCGTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTAc                                                                                                                                                                   <  2:9747/250‑1 (MQ=255)                                                                                   acTGGCGCGCTAACAACCTGTTCGATGGCTGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTTTTCGCGCTGCTCTTCAACCTGTACAGCGTGCGGACGTTCAAATGGATAGCGTATCCAGACTTTGCCAGAGGCCAGTTTCGAAAACGCGCACGCTTCGGTCATTGGCAAATGCGAGAGCGGTCGATAACGCTTGTCACGATACCGATGACGACGCTTACCACGTACGCGCAg                                                                                                                                                             <  2:90106/251‑1 (MQ=255)                                                                                          cgcTAACAACCGTTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTGCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGTCCGCTTACGCGCAGGTGACGa                                                                                                                                                      <  2:52280/251‑1 (MQ=255)                                                                                          cgcTAACAACCGGTTCGATAGAGGCGGCGACAGGGAACTCACTCACCATCGGCCGCACATGTACTTTTTCCTGTTCGCGCTGCCCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCCGACTTTGCCAGAGGCCAGTTCCGCAGACGCGCACGCCACGGTCAACGGCATTGGCGACTGGGGTGGATAACGCTCGTCCAGATAGCGACGACGACGCTGACCACTTAAGCGCCGGGGACgg                                                                                                                                                      >  2:246397/1‑250 (MQ=255)                                                                                                                        cAGGGACATCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGCCCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGcc                                                                                                                         >  2:104548/1‑250 (MQ=255)                                                                                                                                   gTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACg                                                                                                             <  1:2873/251‑1 (MQ=255)                                                                                                                                                      gTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCCTTGTTctct                                                                                          >  2:169760/1‑251 (MQ=255)                                                                                                                                                      gTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTctct                                                                                          <  1:169760/251‑1 (MQ=255)                                                                                                                                                        aCTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCtctt                                                                                         <  1:198922/250‑1 (MQ=255)                                                                                                                                                              tcCTGTTCGCGCTGCTCTTCAACCTGTACATCTGGCGTACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATTCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCt                                                                                   <  2:90291/250‑1 (MQ=255)                                                                                                                                                                  ggtCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTc                                                                               <  2:67334/248‑1 (MQ=255)                                                                                                                                                                     cgcgCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGt                                                                           <  1:261774/251‑1 (MQ=255)                                                                                                                                                                           gctCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGt                                                                      <  2:265732/250‑1 (MQ=255)                                                                                                                                                                                tCAACCTGTACAGCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTgc                                                                  <  2:8397/249‑1 (MQ=255)                                                                                                                                                                                  aaCCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTa                                                              <  1:157484/251‑1 (MQ=255)                                                                                                                                                                                                           aTTGGATAGCGGATCCAGACTTTGCCAGCGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCCGTTTGCGCTACGCCTGGCAGCGGGCCTTTCAcccg                                     >  2:58518/1‑249 (MQ=255)                                                                                                                                                                                                           aTTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTCAcccg                                     >  1:20277/1‑249 (MQ=255)                                                                                                                                                                                                                  aGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTCACCAGTTCtgtg                              <  1:129039/251‑1 (MQ=255)                                                                                                                                                                                                                      gATCCAGACTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCATCGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGCCTGGCAGCGGGACTTTCACCAGTTCTTTGCGTg                          >  2:42010/1‑251 (MQ=255)                                                                                                                                                                                                                             aCTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGa                                                                                                                                                         <  2:128790/117‑1 (MQ=255)                                                                                                                                                                                                                             aCTTTGCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGa                                                                                                                                                         >  1:128790/1‑117 (MQ=255)                                                                                                                                                                                                                                  gCCAGAGGCCAGTTCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGCCTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTCACCAGTTCTGTGCGTGGAGCTGGCGCtt              <  1:193048/251‑1 (MQ=255)                                                                                                                                                                                                                                              ttCCGAAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTCACCAGTTCTGTGCGTGGAGCTGGCGCTTCCTGCACAATTg  >  2:5071/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      AGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGACACCGGCCACTTCTTCAACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTCACCAGTTCTGTGCGTGGAGCTGGCGCTTCCTGAACAATTG  >  NC_000913/1141623‑1142110

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: ‑