Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,734,543	(C)8→9	coding (492/732 nt)	yfiH ←	UPF0124 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,734,535	1	.	C	95.6%	155.3 / 0.1	45	coding (500/732 nt)	yfiH	UPF0124 family protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (25/18);  total (26/19)
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.98e-01

GTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑C‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATGAGCGGCGGCG  >  NC_000913/2734385‑2734675                                                                                                                                                        |                                                                                                                                              gTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCt                                                                                                                                                <  1:201034/151‑1 (MQ=255)  tACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑c                                                                                                                                              >  2:432786/1‑151 (MQ=255)    cAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATAACCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑c‑cc                                                                                                                                           >  1:379467/1‑151 (MQ=255)      aCGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑c‑cccc                                                                                                                                         >  1:805389/1‑151 (MQ=255)               gccgAAAATTTGCTCAACCCCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTc                                                                                                                                 <  2:249955/151‑1 (MQ=255)                    aaaTTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAAc                                                                                                                             >  1:555212/1‑150 (MQ=255)                     aaTTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAAcgc                                                                                                                           >  1:544699/1‑151 (MQ=255)                      aTTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAAcgc                                                                                                                           >  2:171923/1‑150 (MQ=255)                                      cGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGcc                                                                                                           <  2:293524/150‑1 (MQ=255)                                      cGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCg                                                                                                          >  2:225983/1‑151 (MQ=255)                                      cGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCGGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCg                                                                                                          <  1:503825/151‑1 (MQ=255)                                         tCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                                       >  1:400819/1‑151 (MQ=255)                                           gCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCt                                                                                                     >  2:139508/1‑151 (MQ=255)                                            ccAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTACCCCCATGCTAAATAAAAGCTGCCCTTGCTTTGGCGTCTACGGCCATAACCGCCCCGCGACCCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTa                                                                                                    <  1:475526/151‑1 (MQ=255)                                              aGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                                                                  <  2:64680/151‑1 (MQ=255)                                                  gCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                                                              <  2:416602/151‑1 (MQ=255)                                                    tGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCg                                                                                            >  2:449148/1‑151 (MQ=255)                                                       cGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgaga                                                                                         <  2:221889/151‑1 (MQ=255)                                                             aGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAttt                                                                                   <  2:729719/151‑1 (MQ=255)                                                             aGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAttt                                                                                   >  1:285324/1‑151 (MQ=255)                                                                tGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTcc                                                                                >  1:704238/1‑151 (MQ=255)                                                                 gATAAATATCCGCCAGATACTTATCACCATGCGTAATGAAAGCTGCACTTGCTTTAGCGTCTACAGCCATAAACGCCTCGCGAGCCTCC‑CCCCCCCACGTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTc                                                                                 <  1:550704/149‑1 (MQ=255)                                                                    aaaTATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGAt                                                                            <  1:738627/151‑1 (MQ=255)                                                                        atCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATc                                                                        >  1:363503/1‑151 (MQ=255)                                                                                   aCTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                             >  1:638289/1‑151 (MQ=255)                                                                                      tATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa                                                           <  2:668085/150‑1 (MQ=255)                                                                                       aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa                                                          <  2:103073/150‑1 (MQ=255)                                                                                       aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAAc                                                         <  2:375176/151‑1 (MQ=255)                                                                                               gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CTCTCCCCCACTTCGAACGCGCGTGGTCCAATTGCCGGCCCTAACCAGGCGCGAATATTTTCCCGATTATCAGCAAAACAGGAAGCCGtctc                                                   <  1:563243/149‑1 (MQ=255)                                                                                                   aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTcc                                               >  2:418279/1‑149 (MQ=255)                                                                                                   aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAg                                             <  2:156946/151‑1 (MQ=255)                                                                                                     tgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGc                                            >  1:59202/1‑150 (MQ=255)                                                                                                      gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCAc                                          >  1:417845/1‑151 (MQ=255)                                                                                                                   gCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAg                              >  2:803975/1‑150 (MQ=255)                                                                                                                      ttAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTcc                           <  1:323978/150‑1 (MQ=255)                                                                                                                       tAGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCAc                         >  2:95774/1‑151 (MQ=255)                                                                                                                        aGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCAc                         >  2:649795/1‑150 (MQ=255)                                                                                                                        aGCGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACg                        >  2:402132/1‑151 (MQ=255)                                                                                                                          cGTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgcc                      >  2:332171/1‑151 (MQ=255)                                                                                                                           gTCTACTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgcc                      >  2:799188/1‑150 (MQ=255)                                                                                                                              tACTGCCATAAACGCCTCGCGAACCT‑C‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgccagcc                  >  1:387694/1‑150 (MQ=255)                                                                                                                               aCTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgccagcc                  >  2:339052/1‑150 (MQ=255)                                                                                                                               aCTGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgccagcc                  >  2:823032/1‑150 (MQ=255)                                                                                                                                 tGCCATAAACGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCCCGCcagccagc               >  1:436779/1‑151 (MQ=255)                                                                                                                                        aaCGCCTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATGAgcg        <  2:119478/151‑1 (MQ=255)                                                                                                                                            ccTCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATGAgcggcgg    <  2:29351/151‑1 (MQ=255)                                                                                                                                              tCGCGAACCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATGAgcggcggcg  >  1:350495/1‑151 (MQ=255)                                                                                                                                                   aaCCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATGAgcggc      >  1:785466/1‑142 (MQ=255)                                                                                                                                                   aaCCTCC‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATGAgcggc      <  2:785522/142‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                              GTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑C‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATGAGCGGCGGCG  >  NC_000913/2734385‑2734675

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