Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,916,850	Δ1 bp	8.0%	coding (103/369 nt)	PP_1717 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,916,850	0	G	.	8.0%	156.4 / 7.0	50	coding (103/369 nt)	PP_1717	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/26);  new base . (2/2);  total (22/28)
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

TGGCCTGGCATCTGGGGGTGGCTGGCTGGAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGCCGAGCAGCCATTGCGCAGTTCGCTGGTGATCAGCCAAGGCGCCAGCCGGTTGCCGCGTACTGGC  >  NC_002947/1916708‑1916993                                                                                                                                                |                                                                                                                                                tggcctggcATCTGGGGGTGGCTGGCTGGAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCAGTCGTCAACGTCACTTGCTGAGATACCGAGCGCT‑GATGGAGCGCt                                                                                                                                           >  6:295523/1‑149 (MQ=255)   gcgtggcATCTGGGGGTGGCTGGCTGGAGCGAGCCATCATGAGGCAGGTCAGGAGCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCTCGATC‑AGCGCtcg                                                                                                                                         <  3:277302/146‑3 (MQ=255)     ctggcATCTGGGGGTGGCTGGCTGGAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCTCGATC‑AGCGCtcgg                                                                                                                                        >  4:238386/1‑145 (MQ=255)     ctggcATCTGGGGGTGGCTGGCTGGAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCTCGATC‑AGCGCtcgg                                                                                                                                        >  3:69854/1‑145 (MQ=255)     ctggcATCTGGGGGTGGCTGGCTGGAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCa                                                                                                                                       <  5:8958/149‑1 (MQ=255)        gcATCTGGGGGTGGCTGGCTGGAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCTCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGc                                                                                                                                     <  5:290710/148‑1 (MQ=255)        gcATCTGGGGGTGGCTGGCTGGAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCTCGATC‑AGCGCtcggta                                                                                                                                      <  3:238386/147‑6 (MQ=255)              gggggTGGCTGGCTGGAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCtt                                                                                                                                          <  8:99437/137‑1 (MQ=255)              gggggTGGCTGGCTGGAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCtt                                                                                                                                          >  7:99437/1‑137 (MQ=255)                gggTGGCTGGCTGGAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGcc                                                                                                                            >  7:103151/1‑149 (MQ=255)                   tggctggcTGGAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCgg                                                                                                                         <  8:103151/149‑1 (MQ=255)                         gcttgAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTGGTCAACGTCACTTTCTGACATACCGAGCGCT‑GATCGTGGGCTTGCAGCTACCCAGCCCGGCCATg                                                                                                                    <  1:138621/148‑1 (MQ=255)                         gctggAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCTTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCTGCCATGc                                                                                                                   <  3:28598/149‑1 (MQ=255)                             gagcgagcCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGTGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGAc                                                                                                               <  3:52952/149‑1 (MQ=255)                                  agcCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGa                                                                                                          >  8:51076/1‑149 (MQ=255)                                     catcatGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTg                                                                                                       >  7:243796/1‑149 (MQ=255)                                      atcatGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGc                                                                                                      >  2:70407/1‑149 (MQ=255)                                         atGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCcgc                                                                                                   <  1:70407/149‑1 (MQ=255)                                             gCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCt                                                                                                >  4:84919/1‑148 (MQ=255)                                             gCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCt                                                                                                <  3:84919/148‑1 (MQ=255)                                             gCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTg                                                                                               <  2:7090/149‑1 (MQ=255)                                             gCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTg                                                                                               <  1:283182/149‑1 (MQ=255)                                                 ggTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAAt                                                                                           <  7:51076/149‑1 (MQ=255)                                                   tCAGGATCTGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCAGCCATGCAGACGCTGACTGCCGCGCTGGAATAc                                                                                         <  5:74044/149‑1 (MQ=255)                                                      ggATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTg                                                                                      >  6:122379/1‑149 (MQ=255)                                                            gCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTg                                                                                >  8:247383/1‑149 (MQ=255)                                                            gCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTg                                                                                <  5:122379/149‑1 (MQ=255)                                                                  aGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGAt                                                                          <  2:132161/149‑1 (MQ=255)                                                                          gATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTg                                                                             <  5:118352/138‑1 (MQ=255)                                                                          gATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTg                                                                             >  6:118352/1‑138 (MQ=255)                                                                                         gAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGCCGAGCAGCCATTgcg                                                   >  4:333054/1‑149 (MQ=255)                                                                                           gCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGCCGAGCAGCCATTGCGCa                                                 <  2:18471/149‑1 (MQ=255)                                                                                               gACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGAt                                                                          <  3:154173/120‑1 (MQ=255)                                                                                               gACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGAt                                                                          >  4:154173/1‑120 (MQ=255)                                                                                                     gCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGCCGAGCAGCCATTGCGCAGTTCGCTgg                                        <  3:137936/148‑1 (MQ=255)                                                                                                          tcgtcAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGcc                                                                 >  1:109502/1‑118 (MQ=255)                                                                                                          tcgtcAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGcc                                                                 <  2:109502/118‑1 (MQ=255)                                                                                                           cgtcAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGCCGAGCAGCCATTGCGCAGTTCGCTGGTGATCAg                                 >  5:107946/1‑149 (MQ=255)                                                                                                               aaCGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGc                                                                  <  4:125216/112‑1 (MQ=255)                                                                                                               aaCGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGc                                                                  >  3:125216/1‑112 (MQ=255)                                                                                                                aCGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCAGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGCCGAGCAGCCATTGCGCAGTTCGCTGGTGATCAGCCAAg                            >  3:134626/1‑149 (MQ=255)                                                                                                                  gTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGCCGAGCAGCCATTGCGCAGTTCGCTGGTGATCAGCCAAGgc                          <  3:333054/149‑1 (MQ=255)                                                                                                                      cTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGAt                                                                          <  2:242940/97‑1 (MQ=255)                                                                                                                      cTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGAt                                                                          >  1:242940/1‑97 (MQ=255)                                                                                                                       ttGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCACAGCCCGGCCATGCATACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATTACCGTGCCGAGCCGCCATTGAGCAGTTCGATGGTGATAAGACAAGGAGCCAg                     >  1:97802/1‑149 (MQ=255)                                                                                                                        tGCTGACATACCGAGCGCT‑GATCGAGCGCATGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCGGGAATACCTGGCGGTGCTGGATGCCCGTGCCGAGCAGCCATTGCGCAGTTCGCTGGTGATCAGCCAAGGCGCCAGc                    >  6:263704/1‑149 (MQ=255)                                                                                                                         gCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCACGCAGACGCTGACTGCCGCGCTGGAATACCt                                                                                       >  8:190826/1‑81 (MQ=255)                                                                                                                          cTGACATAGCGAGCGCT‑GATCGAGCGCTTGCAGCTGCCCAGCCGGGGCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGCCGAGCAGCCATTGCGCAGTTCGCTGGTGATCAGCCAAGGCGCCAGCCg                  <  4:134626/149‑1 (MQ=255)                                                                                                                               aTACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGCCGAGCAGCCATTg                                                     <  2:289571/109‑1 (MQ=255)                                                                                                                                             aTCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTTGAATACCTGGCGGTGCTGGATGCCCGTGCCGAGCAGCCATTGCGCAGTTCGCTGGTGATCAGCCAAGGCGCCAGCCGGTTGCCGCGTACTGGc  <  8:243796/147‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                TGGCCTGGCATCTGGGGGTGGCTGGCTGGAGCGAGCCATCATGAGCCAGGTCAGGATCAGCCGTCAGCACTGGGATAACCTGCTGGGCGAGCTTGACCTGGCCCGTCGTCAACGTCACTTGCTGACATACCGAGCGCT‑GATCGAGCGCTTGCAGCTCCCCAGCCCGGCCATGCAGACGCTGACTGCCGCGCTGGAATACCTGGCGGTGCTGGATGCCCGTGCCGAGCAGCCATTGCGCAGTTCGCTGGTGATCAGCCAAGGCGCCAGCCGGTTGCCGCGTACTGGC  >  NC_002947/1916708‑1916993

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