Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,535,527	Δ1 bp	17.6%	coding (332/894 nt)	yddE ←	PhzC‑PhzF family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,535,527	0	T	.	17.6%	88.0 / 15.8	34	coding (332/894 nt)	yddE	PhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base T (15/13);  new base . (3/3);  total (18/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 = 9.51e-01

CTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAA  >  NC_000913/1535381‑1535672                                                                                                                                                   |                                                                                                                                                   cTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGGTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGAACCTTGTTCCAGCGAAATACGAGAATCATCGTtg                                                                                                                                                  >  2:834706/1‑149 (MQ=255)  tgttgtTGCCACCTGAATTGGCAAGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTtgt                                                                                                                                                 <  1:246463/149‑1 (MQ=255)            aCCTGAATTGGCAAGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCG‑TGTGCTTTTCGATa                                                                                                                                      >  1:472270/1‑149 (MQ=255)                     gtccaGCCTGGCAGAATATCGTCCTCAGTGTGATGAAGCGCGGTGATAATCGCCGCCCGTGTTTCACCTTCCAGCGGTGGCTCACAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCAGCGTTGTGCTTTTCGATAGTCAcgcg                                                                                                                              <  2:730723/145‑1 (MQ=255)                      gCAAGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCACACGTGTGTCACCATCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGCTAATCATTG‑TGTGCTTTTAGATAGTCACGCGAt                                                                                                                            >  1:523021/1‑149 (MQ=255)                       cAAGCCGGGCAGAATATCGCCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTTCAGCGAAATACGATAATCCTCGTTGTGCTTTTCGATAGTCACGCGAt                                                                                                                            <  1:480738/149‑1 (MQ=255)                              ggCAGAATATCGTCCTCAGGGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCG‑TGTGCTTTTCGATAGTCACGCGATGTTTTCCt                                                                                                                    <  2:391019/149‑1 (MQ=255)                                  gAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGcc                                                                                                                 >  1:613152/1‑149 (MQ=255)                                      atCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGACATACTATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCg                                                                                                             >  1:199191/1‑149 (MQ=255)                                          tCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCCTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTAACGCGATGGTTTCCTGCAGCCGATGt                                                                                                         >  2:517162/1‑149 (MQ=255)                                            cTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGttt                                                                                                       <  1:573320/149‑1 (MQ=255)                                                          gcgcGTTGAGACTCGCCGCACGTGTTTCACCTTCCAGCTGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCaa                                                                                        <  2:320424/149‑1 (MQ=255)                                                              gTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGAGATGTTTGCCAG‑ATCGTGCAATTTc                                                                                    >  1:820651/1‑149 (MQ=255)                                                               ttGATAATCGCCGAACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCAGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTcc                                                                                   >  1:583835/1‑149 (MQ=255)                                                                tGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAC‑ATCGTGCAATTTCCt                                                                                  >  1:830090/1‑149 (MQ=255)                                                                     aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCG‑TGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCGG‑ATCGTGCAATTTCCTAAACCt                                                                            >  2:697242/1‑149 (MQ=255)                                                                                      aCCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCCACCTCCTTCCCACGTa                                                            >  2:105465/1‑149 (MQ=255)                                                                                       ccTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACGCCGTCGCACGTAc                                                           >  2:668972/1‑149 (MQ=255)                                                                                            cAGCGGTGGCTCAAATCCCGGCGTACCTGGTTCCAGCGAAATACGATAACCATCGTTGTGCTTTTCGATAGTCACGCGATGTGTTCCTGCCAGCGCTGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGt                                                      <  1:29336/149‑1 (MQ=255)                                                                                             aGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCAGGCGATGTTTTCCTGCCAGCGAGGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAATg                                                     >  1:58597/1‑149 (MQ=255)                                                                                             aGCGGCGGCTCAAAGCCCGGCGTACCTTGGTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTTAACCTAACACCTTCGCACGTACATAGTg                                                     <  1:517162/149‑1 (MQ=255)                                                                                              gCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCCAAATACGAGAATCATCGTTG‑GCGTTGCAATAGTCACGCTCTGTTTTCCCGCCAGCGATGTGTGCCAGAAGAGTGC‑ATTTCCTAAACCTACCACATTCGCACGTGCATAGtcgg                                                   >  2:292571/1‑146 (MQ=255)                                                                                                      tCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTAc                                            >  1:642450/1‑149 (MQ=255)                                                                                                      tCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTAc                                            >  1:756824/1‑149 (MQ=255)                                                                                                                  gTACCTTGTTCCAGCGAAATACGATAATCATCG‑TGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCg                               <  1:70890/149‑1 (MQ=255)                                                                                                                      cTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCaa                            <  2:199191/149‑1 (MQ=255)                                                                                                                           tCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTTGCTCGTGACCGCAAATCgg                       <  2:613152/149‑1 (MQ=255)                                                                                                                             cAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCa                     <  1:505582/149‑1 (MQ=255)                                                                                                                              aGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCACTTTCCTCAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCAc                    <  2:779686/149‑1 (MQ=255)                                                                                                                                        gATAATCATCG‑TGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTg         <  2:472270/149‑1 (MQ=255)                                                                                                                                             gcatcGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAGCCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACTGCAAATCGGCACTTCAACCGTTGGCGt     >  1:596712/2‑149 (MQ=255)                                                                                                                                                tcGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCACATCGGCACTTCAACCGTTGGCGTaaa  >  2:16943/1‑149 (MQ=255)                                                                                                                                                tcGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTaaa  <  1:701608/149‑1 (MQ=255)                                                                                                                                                acGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTaaa  <  1:402738/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG‑ATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAA  >  NC_000913/1535381‑1535672

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑