Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,352,167	G→A	6.9%	A13V (GCG→GTG)	PP_3824 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,352,167	0	G	A	6.9%	172.3 / 2.2	58	A13V (GCG→GTG)	PP_3824	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base A (2/2);  total (29/29)
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.50e-01

TCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCAACGCGGCAGTGTAACGGCACATCGTGGTATCTTTTG  >  NC_002947/4352029‑4352310                                                                                                                                            |                                                                                                                                                tCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCa                                                                                                                                       <  5:279544/149‑1 (MQ=255)   aCGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACt                                                                                                                                     <  1:183419/149‑1 (MQ=255)     gcgTATCGGCAAGCTTCTCCTTTCCGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCAAATCCATAGTTTCACTCGATAGCTACAGCTCCGCGA‑CGCTCACGCAGCCACAGAGCGAAACCGACAACGCCGACACCATCCACTTc                                                                                                                                   <  4:221370/149‑1 (MQ=255)            ggCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTtgct                                                                                                                            <  4:57536/149‑1 (MQ=255)            ggCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTtgct                                                                                                                            <  3:284199/149‑1 (MQ=255)             gCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCAt                                                                                                                                 <  4:16389/143‑1 (MQ=255)             gCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCAt                                                                                                                                 >  3:16389/1‑143 (MQ=255)             gCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTtgctg                                                                                                                           <  6:273486/149‑1 (MQ=255)                             cGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCt                                                                                                                >  6:304476/1‑144 (MQ=255)                             cGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCt                                                                                                                <  5:304476/144‑1 (MQ=255)                               cAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGaa                                                                                                         <  6:126544/149‑1 (MQ=255)                               cAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGaa                                                                                                         >  5:126544/1‑149 (MQ=255)                                aGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAAt                                                                                                        >  2:105090/1‑149 (MQ=255)                                       aTAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGATTCACTCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCTACCGACAACGCCGACACCAAGCACTTCATTTGCTGCCCCACTGCCTTCAAGAATTCGCTTc                                                                                                 >  2:170358/1‑149 (MQ=255)                                          ggCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGAc                                                                                              <  8:207385/149‑1 (MQ=255)                                           gcgcGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCtt                                                                                                  >  3:144679/1‑144 (MQ=255)                                           gcgcGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCtt                                                                                                  <  4:144679/144‑1 (MQ=255)                                                 gggCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAAtt                                                                                                       >  3:287015/1‑133 (MQ=255)                                                 gggCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAAtt                                                                                                       <  4:287015/133‑1 (MQ=255)                                                     tctTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCa                                                                                                                         >  3:207877/1‑111 (MQ=255)                                                     tctTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCa                                                                                                                         <  4:207877/111‑1 (MQ=255)                                                         gCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGtt                                                                               >  5:286566/1‑149 (MQ=255)                                                           ggtagaTCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGtt                                                                               <  7:220770/144‑1 (MQ=255)                                                                     atccatcGTTTCACGCGATTGCTCCAGCTCCGCGAGCGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGa                                                                                               >  1:48452/1‑122 (MQ=255)                                                                     atccatcGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGa                                                                                               <  2:48452/121‑1 (MQ=255)                                                                        catcGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAg                                                                                                           >  5:95277/1‑106 (MQ=255)                                                                        catcGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAg                                                                                                           <  6:95277/106‑1 (MQ=255)                                                                         atcGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTtg                                                                                                                              >  6:138658/1‑86 (MQ=255)                                                                         atcGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTtg                                                                                                                              <  5:138658/86‑1 (MQ=255)                                                                         atcGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTgg                                                               >  4:181281/1‑149 (MQ=255)                                                                                 aCGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTccc                                                                       <  6:137639/133‑1 (MQ=255)                                                                                 aCGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTccc                                                                       >  5:137639/1‑133 (MQ=255)                                                                                     gATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGc                                                              >  4:3231/1‑138 (MQ=255)                                                                                           tCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCa                                                             >  8:188356/1‑133 (MQ=255)                                                                                           tCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCa                                                             <  7:188356/133‑1 (MQ=255)                                                                                           tCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCg                                             >  6:43008/1‑149 (MQ=255)                                                                                             cAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCATCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGATTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGc                                           >  1:295022/1‑149 (MQ=255)                                                                                             cAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCgg                                            <  6:286566/148‑1 (MQ=255)                                                                                              aGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGc                                           >  2:74489/1‑148 (MQ=255)                                                                                              aGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCACCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAAtt                                                        >  3:149922/1‑135 (MQ=255)                                                                                              aGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCACCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAAtt                                                        >  3:154745/1‑135 (MQ=255)                                                                                              aGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCACCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAAtt                                                        <  4:154745/135‑1 (MQ=255)                                                                                              aGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCACCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAAtt                                                        <  4:149922/135‑1 (MQ=255)                                                                                                   cgcgA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCCCCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCaa                                     >  5:331238/1‑149 (MQ=255)                                                                                                   cgcgA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCaa                                     >  5:330597/1‑149 (MQ=255)                                                                                                         cGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCAACGCgg                                <  2:44285/149‑1 (MQ=255)                                                                                                           cTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCAACGCGGCa                              <  1:74489/149‑1 (MQ=255)                                                                                                            tCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCt                                                                   <  5:303914/111‑1 (MQ=255)                                                                                                            tCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCt                                                                   >  6:303914/1‑111 (MQ=255)                                                                                                             cACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCAACGCGGCAgt                            >  7:234028/1‑149 (MQ=255)                                                                                                              aCGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGc                                                              >  4:24028/1‑114 (MQ=255)                                                                                                              aCGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGc                                                              <  3:24028/114‑1 (MQ=255)                                                                                                              aCGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCAACGCGGCAgtg                           >  8:305467/1‑149 (MQ=255)                                                                                                               cGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCAACGCGGCAgtgt                          <  5:43008/149‑1 (MQ=255)                                                                                                                               aaCCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCAACGCGGCAGTGTAACGGCACATCGTgg           >  2:115674/1‑148 (MQ=255)                                                                                                                               aaCCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCAACGCGGCAGTGTAACGGCACATCGTGGt          >  6:193429/1‑149 (MQ=255)                                                                                                                                      aaCGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCAACGCGGCAGTGTAACGGCACATCGTGGTATctttt   <  6:330597/149‑1 (MQ=255)                                                                                                                                       aCGGCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCAACGCGGCAGTGTAACGGCACATCGTGGTATcttttg  <  1:115674/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                TCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAGATCACATCCATCGTTTCACGCGATTGCTCCAGCTCCGCGA‑CGCTCACGCAGCCACCGAGCGCAACCGACAACGCCGCCACCATCCACTTCATTTGCTGCACCTCTGCCTTCAAGAATTCGCTTCGACCACATCCAGGCGGTTTTGCTCCCTGCTTTGGCAGAATTATCTACGCCCGGCGATCAACGCGGCAGTGTAACGGCACATCGTGGTATCTTTTG  >  NC_002947/4352029‑4352310

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