Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	804,720	G→A	8.0%	D870D (GAC→GAT)	PP_0693 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	804,720	0	G	A	8.0%	143.5 / 2.1	50	D870D (GAC→GAT)	PP_0693	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base A (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-01

CGCCTCGTCGATGTAGATGCCGTTCGGCCCGCTGCCACCTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATCGTCATGGTGGTGCCGTTGGACGCCGCACCGTCGATGTCG  >  NC_002947/804573‑804864                                                                                                                                                    |                                                                                                                                                 cGCCTCGTCGATGTAGATGCCGTTCGGCCCGCTGCCACCTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTtcgt                                                                                                                                                 <  3:265195/149‑1 (MQ=255)  gCCTCGTCGATGTAGATGCCGTTCGGCCCGCTGCCACCTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTtcgtc                                                                                                                                                <  2:30436/149‑1 (MQ=255)  gCCTCGTCGATGTAGATGCCGTTCGGCCCGCTGCCACCTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGAACGAAGCCCAGTTGCTGGCTTGTCTCCAGGAGTGCCGCCTGGTTGCGTtcgtc                                                                                                                                                <  2:17156/149‑1 (MQ=255)           aTGTAGATGCCGTTCGGCCCGCTGCCACCTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACTAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTTGTTGCGTTCGTCGATGTCTGc                                                                                                                                       <  3:109959/149‑1 (MQ=255)           aTGTAGATGCCGTTCGGCCCGCTGCCACCTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGc                                                                                                                                       <  1:168089/149‑1 (MQ=255)            tGTAGATGCCGTTCGGCCCGCTGCCACCTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGcc                                                                                                                                      <  2:208380/149‑1 (MQ=255)             gTAGATGCCGTTCGGCCCGCTGCCACCTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGGTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGcc                                                                                                                                      <  4:253/148‑1 (MQ=255)                      gTTCGGCCCGCTGCCACCTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTc                                                                                                                              >  4:49722/1‑147 (MQ=255)                       ttCGGCCCGCTGCCACCTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAg                                                                                                                           <  4:180495/149‑1 (MQ=255)                           gCCCGCTGCCACCTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGgtag                                                                                                                       >  5:211885/1‑149 (MQ=255)                                       cTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCg                                                                                                            >  8:118481/1‑148 (MQ=255)                                                 gATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCg                                                                                                  <  6:211885/148‑1 (MQ=255)                                                   tCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACgg                                                                                                                  <  8:80569/130‑1 (MQ=255)                                                   tCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACgg                                                                                                                  >  7:80569/1‑130 (MQ=255)                                                               ccGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACgg                                                                                      <  7:245664/146‑1 (MQ=255)                                                               ccGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACgg                                                                                      >  8:245664/1‑146 (MQ=255)                                                                      cTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGc                                                                                                >  8:48711/1‑129 (MQ=255)                                                                      cTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCATCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGt                                                                            >  5:327994/1‑149 (MQ=255)                                                                      cTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCATCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGt                                                                            >  5:263720/1‑149 (MQ=255)                                                                        gtgGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTg                                                                           >  7:204667/1‑148 (MQ=255)                                                                        gtgGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTg                                                                           >  1:114783/1‑148 (MQ=255)                                                                           gTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCATCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTgcag                                                                        <  6:327994/148‑1 (MQ=255)                                                                           gTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCATCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTgcag                                                                        <  6:263720/148‑1 (MQ=255)                                                                              tCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCaa                                                                     >  4:272944/1‑148 (MQ=255)                                                                              tCACACTCGCCAGGATCGGCACGAAGCCCAGGTGCTGGCTGGTCTCCAGCAGTGCCTCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGGCACGGTCCATCAGGTGCAGCaa                                                                     <  5:146593/148‑1 (MQ=255)                                                                                acacTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACa                                                                   <  8:109298/148‑1 (MQ=255)                                                                                acacTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACa                                                                   <  7:118481/148‑1 (MQ=255)                                                                                  acTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTgcagc                                                                       <  5:101890/142‑1 (MQ=255)                                                                                  acTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTgcagc                                                                       >  6:101890/1‑142 (MQ=255)                                                                                  acTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGt                                                                <  2:225136/149‑1 (MQ=255)                                                                                    tCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTa                                                               >  1:155788/1‑148 (MQ=255)                                                                                       ccAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAgg                                                                 <  4:282875/143‑1 (MQ=255)                                                                                       ccAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAgg                                                                 >  3:282875/1‑143 (MQ=255)                                                                                        cAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGt                                                          <  2:155788/149‑1 (MQ=255)                                                                                          ggATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCa                                                                                                                            <  2:193763/81‑1 (MQ=255)                                                                                          ggATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCa                                                                                                                            >  1:193763/1‑81 (MQ=255)                                                                                             tCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGt                                                      >  4:50685/1‑148 (MQ=255)                                                                                                  aCGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCAGCAGGTGCAGCACCAGGTACATGTTGGTCAgcgc                                                >  4:193520/1‑149 (MQ=255)                                                                                                                 tggctggTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATc                                         >  7:35654/1‑141 (MQ=255)                                                                                                                 tggctggTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATc                                         <  8:35654/141‑1 (MQ=255)                                                                                                                        tctcCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATCGTCATGGTGGTGCCg                          >  2:168328/1‑149 (MQ=255)                                                                                                                          tcCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATCGTCATGGTGGTGCCGtt                        >  5:150182/1‑149 (MQ=255)                                                                                                                          tcCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATCGTCATGGTGGTGCCGtt                        >  8:161594/1‑149 (MQ=255)                                                                                                                            cagcagTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGa                                           >  6:191939/1‑128 (MQ=255)                                                                                                                            cagcagTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGa                                           <  5:191939/128‑1 (MQ=255)                                                                                                                            cagcagTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATCGTCATGGTGGTGCCGtt                        <  3:272944/147‑1 (MQ=255)                                                                                                                            cagcagTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATCGTCATGGTGGTGCCGTTgg                      <  2:114783/149‑1 (MQ=255)                                                                                                                                     cgccTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGACCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATCGTCATCGTGGTGCCGTTGGACGCAGCAc             >  6:38398/1‑149 (MQ=255)                                                                                                                                      gccTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGGCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATCGTCATGGTGGTGCCGTTGGACGCCGCAcc            <  3:193520/149‑1 (MQ=255)                                                                                                                                       ccTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATCGTCATGGTGGTGCCGTTGGACGCCGCACCg           >  4:78925/1‑149 (MQ=255)                                                                                                                                            ttGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATCGTCATGGTGGTGCCGTTGGACGCCGCACCGTCGAt      <  7:161594/149‑1 (MQ=255)                                                                                                                                               cGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTCCAGCAACAGGTACATGTTGGTCAGCTCCTTGATCGTCATGGAGGTGCCGTTGGACCCCGCACCGTCGATGTc   >  7:149166/1‑149 (MQ=255)                                                                                                                                                gTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATCGTCATGGTGGTGCCGTTGGACGCCGCACCGTCGATGTCg  >  7:299689/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CGCCTCGTCGATGTAGATGCCGTTCGGCCCGCTGCCACCTTCCAGGTCGATCGCAACGTGCGCCGAGACCTGTGGTTTCACACTCGCCAGGATCGGCACGAAGCCCAGTTGCTGGCTGGTCTCCAGCAGTGCCGCCTGGTTGCGTTCGTCGATGTCTGCCGCTTCGTCCAGGTAGTACGGCAGGCGAATGCGCCCGGCCAGGTCACGGTCCATCAGGTGCAGCAACAGGTACATGTTGGTCAGCGCCTTGATCGTCATGGTGGTGCCGTTGGACGCCGCACCGTCGATGTCG  >  NC_002947/804573‑804864

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