Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,735,143:1	+T	11.8%	coding (40/3093 nt)	PP_3302 →	RND family transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,735,143	1	.	T	11.8%	111.4 / 12.9	34	coding (40/3093 nt)	PP_3302	RND family transporter
Reads supporting (aligned to +/- strand):  ref base . (16/14);  new base T (2/2);  total (18/16)
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.87e-01

ACCTGCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGGCGCCACTGCGC  >  NC_002947/3734998‑3735287                                                                                                                                                   |                                                                                                                                                 aCCTGCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAc‑gc                                                                                                                                                <  3:372281/148‑1 (MQ=255) aCCTGCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAc‑gc                                                                                                                                                <  3:284671/148‑1 (MQ=255)  ccTGCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑gcga                                                                                                                                              <  1:326260/149‑1 (MQ=255)   cTGCACGCTGGCAATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑gcgag                                                                                                                                             >  8:138178/1‑149 (MQ=255)     gCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑gcgagc                                                                                                                                            >  4:164194/1‑148 (MQ=255)      cACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCATCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcgc                                                                                                                                          <  8:92971/149‑2 (MQ=255)      cACGCTGGCGATCGGATCGTCGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcgc                                                                                                                                          <  4:218199/149‑2 (MQ=255)       aCGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCTATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGat                                                                                                                                         <  3:164194/149‑1 (MQ=255)               gatcggatcgTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGATTCAACTTAACCGCGCTAGCACTAC‑CCGAGCGATCAATCAcc                                                                                                                                 >  5:308575/1‑149 (MQ=255)                         tGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTc                                                                                                                           >  8:195947/1‑145 (MQ=255)                         tGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTc                                                                                                                           <  7:195947/145‑1 (MQ=255)                                                   gATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGTCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACg                                                                                             >  1:294506/1‑149 (MQ=255)                                                         gTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGcc                                                                                       >  3:362401/1‑149 (MQ=255)                                                            gAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCttt                                                                                    >  2:298587/1‑149 (MQ=255)                                                             aGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCtttt                                                                                   <  5:141033/149‑1 (MQ=255)                                                                  aTGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAg                                                                              >  6:69822/1‑149 (MQ=255)                                                                         aGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGcc                                                                       <  7:138178/149‑1 (MQ=255)                                                                           gCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCaa                                                                               <  4:171109/139‑1 (MQ=255)                                                                           gCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCaa                                                                               >  3:171109/1‑139 (MQ=255)                                                                               aaaCAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTg                                                                    >  8:291644/1‑146 (MQ=255)                                                                                aaCAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCaga                                                                <  2:294506/149‑1 (MQ=255)                                                                                aaCAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCATGATTGTGATGATCGCCTTCGCCGGCACGCTGGACTTTTTCAACCTGGGCCGTGCag                                                                 >  6:320668/1‑148 (MQ=255)                                                                                   aGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTc                                                                                                      <  5:81809/108‑1 (MQ=255)                                                                                   aGCTCAGGTCATGCGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTc                                                                                                      >  6:81809/1‑108 (MQ=255)                                                                                              tGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCa                                                  >  8:10059/1‑149 (MQ=255)                                                                                                gggCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCAcc                                                <  4:362401/149‑1 (MQ=255)                                                                                                gggCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGCCCGTGCAGAAGATCCGCCCTTCAcc                                                >  4:196699/1‑149 (MQ=255)                                                                                                   cACCAACAATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTc                                             <  6:308575/149‑1 (MQ=255)                                                                                                      caaccaTGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATc                                                                                                            <  6:227361/83‑1 (MQ=255)                                                                                                      caaccaTGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATc                                                                                                            >  5:227361/1‑83 (MQ=255)                                                                                                                  aGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGAtt                              <  5:69822/149‑1 (MQ=255)                                                                                                                    cTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTAt                            <  7:291644/149‑1 (MQ=255)                                                                                                                              aCTTATCCTCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTgg                  <  5:320668/149‑1 (MQ=255)                                                                                                                               cTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGc                 >  5:33433/1‑149 (MQ=255)                                                                                                                                      tcgcTCGC‑GTACTGCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGgcgc          >  8:260588/2‑149 (MQ=255)                                                                                                                                      tcgcTCGC‑GTACTGCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGgcgc          >  4:119559/2‑149 (MQ=255)                                                                                                                                          tCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGGCGCCACt      >  2:148012/1‑149 (MQ=255)                                                                                                                                              aGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGGCGCCACTgcgc  <  5:320869/149‑1 (MQ=255)                                                                                                                                              aGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGGCGCCACTgcg   >  6:320869/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                 ACCTGCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGGCGCCACTGCGC  >  NC_002947/3734998‑3735287

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