Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	893,926	G→A	9.7%	D41D (GAC→GAT)	PP_0775 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	893,926	0	G	A	9.7%	140.5 / 2.9	51	D41D (GAC→GAT)	PP_0775	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base A (3/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCAGCCGCAAGCGAAAGGCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCTTTTGCAAAAGGTCTTGCCCCATGCTGATCGCCGCCAAC  >  NC_002947/893782‑894064                                                                                                                                                 |                                                                                                                                            gcACCAGCCGCAAGCGAAAGGCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCt                                                                                                                                         >  3:233423/1‑149 (MQ=255)  caccagccgccgtggcaggcacgcccctgggtgctcgacccggcacaacgtgcgcaacacatcggccagcctttcgcttgcggctggtgGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCATCCt                                                                                                                                         >  7:109180‑M2/90‑148 (MQ=255)  caccagccgccgtggcaggcacgcccctgggtgctcgacccggcacaacgtgcgcaacacatcggccagcctttcgcttgcggctggtgGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCATCCt                                                                                                                                         >  7:114512‑M2/90‑148 (MQ=255)    ccagccgccgtggcaggcacgcccctgggtgctcgacccggcacaacgtgcgcaacacatcggccagcctttcgcttgcggctggtgGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCtgat                                                                                                                                      >  1:36389‑M2/88‑149 (MQ=255)    CCAGCCGCAAGCGAAAGGCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCtgat                                                                                                                                      >  5:312953‑M1/1‑149 (MQ=255)    CCAGCCGCAAGCGAAAGGCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCtgat                                                                                                                                      >  2:331185‑M1/1‑149 (MQ=255)       GCCGCAAGCGAAAGGCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCtgatga                                                                                                                                    <  2:219516‑M1/148‑1 (MQ=255)            aaGCGAAAGGCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCg                                                                                                                              >  5:177313/1‑149 (MQ=255)            aaGCGAAAGGCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCg                                                                                                                              <  1:206599/149‑1 (MQ=255)             aGCGAAAGGCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGt                                                                                                                             <  6:230377/149‑1 (MQ=255)                 aaaGGCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGt                                                                                                                             >  3:359164/1‑145 (MQ=255)                 aaaGGCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGt                                                                                                                             <  4:359164/145‑1 (MQ=255)                  aaGGCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATc                                                                                                                        <  4:233423/149‑1 (MQ=255)                     gCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCtgat                                                                                                                                      >  1:357701/1‑132 (MQ=255)                     gCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCtgat                                                                                                                                      <  2:357701/132‑1 (MQ=255)                         gCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGa                                                                                                                 >  1:360747/1‑149 (MQ=255)                          ccGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑t                                                                                                               <  6:177313/149‑1 (MQ=255)                            gATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGTGCAGGTCCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑Tgg                                                                                                             <  1:122893/149‑1 (MQ=255)                            gATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑Tgg                                                                                                             <  4:181472/149‑1 (MQ=255)                                    gcgcACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACTGCGGCTGCTGGAGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGAGGCCTTCTTCATCCTGATGAAAGCCGATGATCCGACTGA‑TGGGCCGTTGc                                                                                                     >  5:256739/1‑149 (MQ=255)                                      gacccggcacaacgtgcgcaacacatcggccagcctttcgcttgcggctggtgGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCATCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTc                                                                                                   <  8:114512‑M2/96‑1 (MQ=255)                                      gacccggcacaacgtgcgcaacacatcggccagcctttcgcttgcggctggtgGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCATCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTc                                                                                                   <  8:109180‑M2/96‑1 (MQ=255)                                             gtgCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGca                                                                                            >  7:124096/1‑149 (MQ=255)                                                ccGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGt                                                                                                                             <  3:17353/114‑1 (MQ=255)                                                ccGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGt                                                                                                                             >  4:17353/1‑114 (MQ=255)                                                ccGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGcacc                                                                                          <  2:360747/148‑1 (MQ=255)                                                 cGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGcaccac                                                                                        >  7:151386/1‑149 (MQ=255)                                                     tCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCATGTTc                                                                                    <  5:307491/149‑1 (MQ=255)                                                      cgcaacacatcggccagcctttcgcttgcggctggtgGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCg                                                                                   <  2:36389‑M2/112‑1 (MQ=255)                                                         gCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCaa                                                                                >  3:206579/1‑149 (MQ=255)                                                           aCCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGATTGGGCCGTTGCTCGTCTGCACCACGTTCGCAAt                                                                               <  1:45258/149‑1 (MQ=255)                                                             ccAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCg                                                                             >  2:157043/1‑148 (MQ=255)                                                                          tgccACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGt                                                                                                                             <  3:238944/88‑1 (MQ=255)                                                                          tgccACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGt                                                                                                                             >  4:238944/1‑88 (MQ=255)                                                                              aCGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTaca                                                           <  8:151386/149‑1 (MQ=255)                                                                                  cggctggtgGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTg                                                       >  1:256861‑M2/10‑149 (MQ=255)                                                                                         tgGCTCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCTCCCAGTGGCATTCTTCGTCCTGATGAAAGCCGTTGATCCGAGTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGTACGATGTCAAGGTTACACTTGTCCGctc                                                <  6:256739‑M2/147‑1 (MQ=255)                                                                                             gcAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAAt                                                                               <  2:344714/114‑1 (MQ=255)                                                                                             gcAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAAt                                                                               >  1:344714/1‑114 (MQ=255)                                                                                             gcAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACa                                            >  6:361839/1‑149 (MQ=255)                                                                                                    cTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCtttt                                     <  2:256861/149‑1 (MQ=255)                                                                                                               cGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCTTTTGCAAAAGGTCt                          >  2:228678/1‑149 (MQ=255)                                                                                                               cGCACGAAAGTTCCGCCACCCAGTGGCCATCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCTTTTGCAAAAGGTCt                          >  1:211242/1‑149 (MQ=255)                                                                                                                 cACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCTTTTGCAAAAGGTCTTg                        <  1:219516/149‑1 (MQ=255)                                                                                                                  aCGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCTTTTGCAAAAGGTCTTGc                       >  1:235862/1‑149 (MQ=255)                                                                                                                  aCGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCTTTTGCAAAAGGTCTTGc                       <  1:157043/149‑1 (MQ=255)                                                                                                                       aGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCTTTTGCAAAAGGTCTTGCCCCAt                  <  2:235862/149‑1 (MQ=255)                                                                                                                           ccGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCTTTTGCAAAAGGTCTTGCCCCATGCt               >  6:228898/1‑148 (MQ=255)                                                                                                                             gCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCTTTTGCAAAAGGTCTTGCCCCATGCTGAt            <  4:206579/149‑1 (MQ=255)                                                                                                                                 cccAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCTTTTGCAAAAGGTCTTGCCCCATGCTGATcgcc        >  6:289475/1‑149 (MQ=255)                                                                                                                                        gCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCTTTTGCAAAAGGTCTTGCCCCATGCTGATCGCCGCcaac  >  4:81941/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                            GCACCAGCCGCAAGCGAAAGGCTGGCCGATGTGTTGCGCACGTTGTGCCGGGTCGAGCACCCAGGGGCGTGCCTGCCACGGCGGCTGGTGGCGCAGGTGCTGGGTATGGCCGCACGAAAGTTCCGCCACCCAGTGGCCTTCTTCGTCCTGATGAAAGCCGTTGATCCGACTGA‑TGGGCCGTTGCTCGTCTGCACCACGTTCGCAATCGGACGATGGCAAGGTTACACTTGTCCGCTCTACATTCTTTTGCAAAAGGTCTTGCCCCATGCTGATCGCCGCCAAC  >  NC_002947/893782‑894064

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: ‑