Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,816,592	(A)8→9	coding (349/429 nt)	yqaA ←	DedA family protein YqaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,816,584	1	.	A	93.2%	142.5 / 2.9	44	coding (357/429 nt)	yqaA	DedA family protein YqaA
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base A (20/21);  total (21/23)
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.71e-01

AATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCAC‑‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATCCAGCTTAATAATAGCGTGACTGCACCATAGCGTTTCAGCCAGCCGGTAGCTTTCTC  >  NC_000913/2816460‑2816718                                                                                                                              |                                                                                                                                        aaTCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTAc                                                                                                                          <  2:551747/140‑1 (MQ=255)   tCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCg                                                                                                                        <  2:459556/140‑1 (MQ=255)    cGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCgg                                                                                                                       >  2:339941/1‑140 (MQ=255)       tttAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTcc                                                                                                                    >  2:61396/1‑140 (MQ=255)               gCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACAAAAAAAAAAGATTACCGGTCCCCACgag                                                                                                             >  1:597485/1‑140 (MQ=255)                        tAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGcatcc                                                                                                   >  1:37469/1‑140 (MQ=255)                         aaTGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTATCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGcatcca                                                                                                  >  2:592472/1‑140 (MQ=255)                          aTGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGcatccat                                                                                                 <  2:120500/140‑1 (MQ=255)                            ggTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCatccatcc                                                                                               <  2:515804/140‑1 (MQ=255)                            ggTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCatccatcc                                                                                               <  2:209847/140‑1 (MQ=255)                            ggTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCatccatcc                                                                                               <  1:61396/140‑1 (MQ=255)                            ggTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCatccatcc                                                                                               <  2:105697/140‑1 (MQ=255)                              tGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTATCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCg                                                                                             <  1:592472/140‑1 (MQ=255)                                aCCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCg                                                                                             >  1:270279/1‑138 (MQ=255)                                aCCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCg                                                                                             <  2:270279/138‑1 (MQ=255)                                      tgtgCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGa                                                                                     >  1:51671/1‑140 (MQ=255)                                       gtgCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGac                                                                                    >  1:393152/1‑140 (MQ=255)                                                       aGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCTGCTAACAGACACAGTAAATCGCCAAc                                                                    >  1:777962/1‑140 (MQ=255)                                                       aGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGcc                                                                       >  2:476394/1‑137 (MQ=255)                                                       aGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGcc                                                                       <  1:476394/137‑1 (MQ=255)                                                       aGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAAc                                                                    >  2:738535/1‑140 (MQ=255)                                                       aGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAAc                                                                    >  2:287046/1‑140 (MQ=255)                                                         tGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCa                                                                  >  2:611078/1‑140 (MQ=255)                                                              ccacATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGcc                                                                       >  2:675711/1‑130 (MQ=255)                                                              ccacATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGcc                                                                       <  1:675711/130‑1 (MQ=255)                                                                 catcatGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATc                                                          <  1:738535/140‑1 (MQ=255)                                                                 catcatGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATc                                                          <  2:51671/140‑1 (MQ=255)                                                                  atcatGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATcc                                                         <  1:468459/140‑1 (MQ=255)                                                                  atcatGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATcc                                                         <  2:388852/140‑1 (MQ=255)                                                                   tcatGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATcc                                                         <  1:61968/139‑1 (MQ=255)                                                                   tcatGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATcc                                                         >  2:61968/1‑139 (MQ=255)                                                                            gAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATcc                                                         >  2:307451/1‑130 (MQ=255)                                                                            gAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATcc                                                         <  1:307451/130‑1 (MQ=255)                                                                            gAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATc                                                          >  2:684436/1‑129 (MQ=255)                                                                            gAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATc                                                          <  1:684436/129‑1 (MQ=255)                                                                             aaCGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCAC‑‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATCCAGCTTAataata                                             >  2:187793/1‑140 (MQ=255)                                                                                      tgctgcAACTGCAACATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATCCAGCTTAATAATAGCGTGACt                                     <  2:37469/140‑1 (MQ=255)                                                                                      tgctgcAACTGCAACATAGCGTAACGCTTTACCAAGGCAC‑‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATCCAGCTTAATAATAGCGTGACTg                                    <  1:187793/140‑1 (MQ=255)                                                                                             actgcaacATAGCGTAACGCTTTACCAAGGCAC‑‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATCCAGCTTAATAATAGCGTGACTGCACCATa                             <  2:621201/140‑1 (MQ=255)                                                                                                  aacATAGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATCCAGCTTAATAATAGCGTGACTGCACCATAGCGt                         >  2:759313/1‑140 (MQ=255)                                                                                                       aGCGTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATCCAGCTTAATAATAGCGTGACTGCACCATAGCGTTTcag                    >  1:752575/1‑140 (MQ=255)                                                                                                          gTAACGCTTTACCAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATCCAGCTTAATAATAGCGTGACTGCACCATAGCGTTTcagcca                 <  1:611078/140‑1 (MQ=255)                                                                                                                     ccAAGGCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATCCAGCTTAATAATAGCGTGACTGCACCATAGCGTTTCAGCCAGCCGGTAGCtt      <  1:759313/140‑1 (MQ=255)                                                                                                                         ggCACA‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATCCAGCTTAATAATAGCGTGACTGCACCATAGCGTTTCAGCCAGCCGGTAGCTTtctc  >  1:366713/1‑140 (MQ=255)                                                                                                                              |                                                                                                                                        AATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCAC‑‑AAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATCGCCAACCACGGGCATCCAGCTTAATAATAGCGTGACTGCACCATAGCGTTTCAGCCAGCCGGTAGCTTTCTC  >  NC_000913/2816460‑2816718

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