Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,142,777	T→G	T531P (ACC→CCC)	rne ←	ribonuclease E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,142,777	0	T	G	97.4%	130.6 / ‑3.1	38	coding (1591/3186 nt)	rne	ribonuclease E
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (18/19);  minor base . (0/1);  total (18/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 = 9.99e-01

AGAAGCGGCTCAACAGCCCAGGTTGTGCAGGAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGTTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTACGCACGCGCAGC  >  NC_000913/1142645‑1142910                                                                                                                                     |                                                                                                                                      agaagCGGCTCAACAGCCCAGGTTGTGCAGGAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTg                                                                                                                                     >  1:751040/1‑135 (MQ=255)        gCTCAACAGCCCAGGTTGTGCAGGAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGgtt                                                                                                                         >  1:1161488/1‑140 (MQ=255)            aaCAGCCCAGGTTGTGCAGGAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAgg                                                                                                                           >  2:408359/1‑134 (MQ=255)            aaCAGCCCAGGTTGTGCAGGAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAgg                                                                                                                           <  1:408359/134‑1 (MQ=255)              cAGCCCAGGTTGTGCAGGAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCg                                                                                                                   <  2:1161488/140‑1 (MQ=255)                   cAGGTTGTGCAGGAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGGGCCAGCGGAGGTTGTTCCGGACGc                                                                                                              >  2:994936/1‑140 (MQ=255)                     ggTTGTGCAGGAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGTGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCtt                                                                                                            <  1:908454/140‑1 (MQ=255)                          tgCAGGAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGtt                                                                                                       >  2:1378327/1‑140 (MQ=255)                          tgCAGGAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGGGCCAGCGCAGGTTGTTCCGGACGCTTACGtt                                                                                                       >  1:889324/1‑140 (MQ=255)                               gAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCg                                                                                                  <  1:1378327/140‑1 (MQ=255)                               gAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCg                                                                                                  <  1:1360576/140‑1 (MQ=255)                                 gctgctGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGaa                                                                                                <  1:312978/140‑1 (MQ=255)                                 gctgctGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGaa                                                                                                <  1:499453/140‑1 (MQ=255)                                    gctgGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAActc                                                                                             >  1:558657/1‑140 (MQ=255)                                      tgGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACtctt                                                                                           <  1:687768/140‑1 (MQ=255)                                         tgtTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCCGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTcttctt                                                                                        >  2:843397/1‑140 (MQ=255)                                                   gCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAg                                                                              >  2:1338238/1‑140 (MQ=255)                                                      gCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAgcgc                                                                           <  2:939858/140‑1 (MQ=255)                                                       cTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGcc                                                                          >  2:847944/1‑140 (MQ=255)                                                                    aGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCAt                                                             <  2:1265691/140‑1 (MQ=255)                                                                          aacaGGCGCGGGAGGTTCAACCGGGGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGAGCCAGCGCAGGTTGTTCCGGACGCCTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCt                                                       >  2:523627/1‑140 (MQ=255)                                                                           acaGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCtt                                                      <  2:254641/140‑1 (MQ=255)                                                                                 gcgGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGgcag                                                >  1:267609/1‑140 (MQ=255)                                                                                     cAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATg                                            <  1:753285/140‑1 (MQ=255)                                                                                            aGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTa                                     <  2:382834/140‑1 (MQ=255)                                                                                                   gTTGGCGCAGGCGGCACATCCGGCATGGCAAAGG‑TGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTggg                             <  1:881070/140‑1 (MQ=255)                                                                                                    ttGGCGCAGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTggg                             <  1:847944/140‑1 (MQ=255)                                                                                                           aGGCGGCACATCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTttctt                      >  2:698934/1‑140 (MQ=255)                                                                                                                 cacaTCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCtt                >  1:620496/1‑140 (MQ=255)                                                                                                                 cacaTCCGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCtt                >  2:62995/1‑140 (MQ=255)                                                                                                                      ccGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTacgc           <  1:843397/140‑1 (MQ=255)                                                                                                                       cGGCATGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTacgca          >  1:153623/1‑140 (MQ=255)                                                                                                                           aTGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTTCGCAcgcg      >  1:1131042/1‑140 (MQ=255)                                                                                                                           aTGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTACGCAcgcg      <  1:1338238/140‑1 (MQ=255)                                                                                                                           aTGGCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTACGCAcgcg      <  1:62995/140‑1 (MQ=255)                                                                                                                           aTGGCAAAGGGTGCCAGCGAAGGTTGTTCCGGACGTTTACGTTAAGCGACCTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTCCGCAcgcg      <  2:889324/140‑1 (MQ=255)                                                                                                                             ggCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTACGCACGCgca    >  2:941966/1‑140 (MQ=255)                                                                                                                             ggCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTACGCACGCgca    >  1:1380405/1‑140 (MQ=255)                                                                                                                              gCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTACGCACGCgca    >  1:1223153/1‑139 (MQ=255)                                                                                                                              gCAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTACGCACGCgca    <  2:1223153/139‑1 (MQ=255)                                                                                                                               cAAAGGGTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTTCGCACGCgcagc  <  2:1131042/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      AGAAGCGGCTCAACAGCCCAGGTTGTGCAGGAGCTGCTGGTGTTGCCGGTGCAGCTTTCGGTGCTGGAGCTACAACAGGCGCGGCAGGTTCAGCTGGCGTTGGCGCAGGCGGCACATCCGGCATGGCAAAGGTTGCCAGCGCAGGTTGTTCCGGACGCTTACGTTCAGCGAACTCTTCTTCAGACGGCAGCGCCATCGCTTCTTCATGCAGCTTCGGCAGCATGTAGCTTAAGGTTGGGGTTTCTTCCCCTTTACGCACGCGCAGC  >  NC_000913/1142645‑1142910

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑