Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,013,205	G→A	9.3%	intergenic (+33/+98)	prfC → / ← potF‑I	peptide chain release factor 3/putrescine‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,013,205	0	G	A	9.3%	117.8 / 2.5	43	intergenic (+33/+98)	prfC/potF‑I	peptide chain release factor 3/putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (15/24);  new base A (2/2);  total (17/26)
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.81e-01

CCATCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTC  >  NC_002947/1013060‑1013350                                                                                                                                                  |                                                                                                                                                  ccATCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCgg                                                                                                                                                <  4:32612/149‑1 (MQ=255) ccATCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCgg                                                                                                                                                >  4:11512/1‑149 (MQ=255)    tCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCGGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCa                                                                                                                                             <  3:206474/149‑1 (MQ=255)      gACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                           >  4:58484/1‑149 (MQ=255)       aCGGCGGCGGGCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGc                                                                                                                                          <  8:247463/149‑1 (MQ=255)            ggcggTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCt                                                                                                                                     <  4:66478/149‑1 (MQ=255)              cggTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                   <  2:124842/149‑1 (MQ=255)                  cACCTGACCTACCGGGCGCCGCCCCGGGTCAACCTGGCGCGGATGGACGAGCGCTGGACGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                               <  4:11036/149‑1 (MQ=39)                   aCCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                              >  5:293125/1‑149 (MQ=255)                     cTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGc                                                                                                                                              >  6:116243/1‑131 (MQ=255)                     cTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGc                                                                                                                                              <  5:116243/131‑1 (MQ=255)                     cTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAgg                                                                                                                             >  1:159764/1‑148 (MQ=255)                       gACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAAt                                                                                                                          >  3:57815/1‑149 (MQ=255)                        acctacctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCTGGCACGCCCGCTCCCACAGGAATc                                                                                                                         <  4:57815/149‑1 (MQ=255)                             cctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAcc                                                                                                                      <  7:331309/147‑1 (MQ=255)                             cctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCa                                                                                                                     <  3:58484/148‑1 (MQ=255)                               tGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAc                                                                                                                       >  6:40827/1‑144 (MQ=255)                               tGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAc                                                                                                                       <  5:40827/144‑1 (MQ=255)                                       aCCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGt                                                                                                          >  6:174230/1‑149 (MQ=255)                                          cGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAt                                                                                                       <  2:159764/149‑1 (MQ=255)                                           gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGATCACCATTAAGGCCTGAATTCGGGTCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGATGTCAtg                                                                                                      <  2:120805/149‑1 (MQ=255)                                                 aCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCgg                                                                                                                                                <  2:40037/101‑1 (MQ=255)                                                 aCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCgg                                                                                                                                                >  1:40037/1‑101 (MQ=255)                                                 aCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc                                                                                                <  3:11512/149‑1 (MQ=255)                                                 aCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCACGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                      <  8:355115/143‑1 (MQ=255)                                                 aCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCACGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                      <  8:359924/143‑1 (MQ=255)                                                 aCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCACGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                      >  7:355115/1‑143 (MQ=255)                                                 aCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCACGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                      >  7:359924/1‑143 (MQ=255)                                                   cTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCt                                                                                              <  1:81084/149‑1 (MQ=255)                                                      gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCTCAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgtg                                                                                           <  3:47127/149‑1 (MQ=25)                                                       cgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCTGGCACGCCCGCTCCCACAGGAATCACCATACCTGGAGGTCATGTGATACCTGTgg                                                                                          <  4:98790/149‑1 (MQ=255)                                                         cTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc                                                                                                >  6:164335/1‑141 (MQ=255)                                                         cTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc                                                                                                <  5:164335/141‑1 (MQ=255)                                                         cTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGa                                                                                        <  6:293125/149‑1 (MQ=255)                                                               gAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGc                                                                                  <  5:174230/149‑1 (MQ=255)                                                                 agagCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGt                                                                                <  7:179754/149‑1 (MQ=255)                                                                  gagCGCTGGCCGGACATCGAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCgg                                                                                                                                                >  1:211432/1‑84 (MQ=255)                                                                  gagCGCTGGCCGGACATCGAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCgg                                                                                                                                                <  2:211432/84‑1 (MQ=255)                                                                                                  ccGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                     >  6:100472/1‑95 (MQ=255)                                                                                                  ccGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                     <  5:100472/95‑1 (MQ=255)                                                                                                                 aGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGa                                 >  8:66771/1‑148 (MQ=255)                                                                                                                                            cccTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                            >  2:247487/1‑62 (MQ=255)                                                                                                                                               ttCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCGGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGCAGGGCAGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTc  >  7:113976/1‑149 (MQ=18)                                                                                                                                                  |                                                                                                                                                  CCATCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTC  >  NC_002947/1013060‑1013350

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