Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	401,361	C→T	7.5%	A316V (GCG→GTG)	PP_0334 →	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	401,361	0	C	T	7.5%	137.5 / 3.0	53	A316V (GCG→GTG)	PP_0334	transposase
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

CCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCT  >  NC_002947/401251‑401507                                                                                                                |                                                                                                                                                   ccTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGAGCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATc                                                                                                               >  1:309820/1‑149 (MQ=2) ccTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCt                                                                                                              <  5:190165/149‑1 (MQ=31)  cTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTg                                                                                                             <  7:89656/149‑1 (MQ=32)  cTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTg                                                                                                             <  5:355928/149‑1 (MQ=32)     atgatCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTggg                                                                                                           <  5:263496/148‑1 (MQ=32)         tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTg                                                                                                      >  1:145248/1‑149 (MQ=32)            aGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTa                                                                                                   <  3:325344/149‑1 (MQ=32)             gggCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTa                                                                                                   >  7:330458/1‑148 (MQ=32)                cgcgcATTCCTTCTTGAGCCAGCTAAGGCCTGATTCCGATCTG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACaga                                                                                               <  5:316933/149‑1 (MQ=17)                    cATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGa                                                                                           <  3:217608/149‑1 (MQ=32)                     aTTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACaga                                                                                               <  7:311720/144‑1 (MQ=32)                     aTTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACaga                                                                                               >  8:311720/1‑144 (MQ=32)                        ccGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCa                                                                                        <  1:129025/148‑1 (MQ=32)                         cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAgg                                                                                      >  5:392126/1‑149 (MQ=32)                         cGTCTTGAGACAGCTAAGGCCTGATTCAGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCATTGGTCGCAAGGAGTCTGCGATCAGGCTGGCTAATCGAAATCTGCGTATTATCTGGGCCCTGTTACAGAGCGATCAgg                                                                                      >  2:253929/1‑149 (MQ=11)                           tCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGAGCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGtt                                                                                    <  8:269901/149‑1 (MQ=17)                               gAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTa                                                                                  <  2:99965/147‑1 (MQ=32)                               gAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTAc                                                                                 <  1:203949/148‑1 (MQ=32)                                aGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTAcc                                                                                <  3:304716/148‑1 (MQ=32)                                 gCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCaa                                                                              >  5:271764/1‑149 (MQ=32)                                  ccAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAAc                                                                             >  6:323339/1‑149 (MQ=32)                                    aGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACg                                                                            >  2:141698/1‑148 (MQ=32)                                       tAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCaaa                                                                        >  3:293602/1‑149 (MQ=32)                                        aaGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGAGATAAGGTTTACAAACGCAAAc                                                                       >  2:316093/1‑149 (MQ=12)                                           gCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcg                                                                     >  1:190296/1‑148 (MQ=32)                                             cTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGTGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgacc                                                                  >  3:203182/1‑147 (MQ=17)                                             cTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac                                                                  >  5:391316/1‑149 (MQ=32)                                             atGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGTGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac                                                                  >  3:203186/2‑149 (MQ=17)                                                aTTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaacga                                                                >  2:363857/1‑148 (MQ=32)                                                    cGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGTCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAAc                                                           <  1:369686/149‑1 (MQ=17)                                                     gATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAAcc                                                          >  2:289207/1‑149 (MQ=32)                                                      aTCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg                                                         >  1:157167/1‑149 (MQ=32)                                                              aGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcgg                                                  >  4:202532/1‑149 (MQ=32)                                                                    tCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGa                                            <  2:145248/149‑1 (MQ=32)                                                                    tCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGa                                            <  2:309820/149‑1 (MQ=32)                                                                      gccgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAACCCGGAGGCGGCTATGAGc                                          >  2:302188/1‑149 (MQ=17)                                                                           ctctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCctct                                      >  2:196652/1‑148 (MQ=32)                                                                           ctctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTg                                     <  2:190296/149‑1 (MQ=32)                                                                              ttGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGa                                    <  8:302942/147‑1 (MQ=32)                                                                              ttGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGa                                    <  1:141698/147‑1 (MQ=32)                                                                                          tCCCGGCTTGGGCGCAAGGAGGTGGCGATCAGGCTGGATAATCGAACTCTGCGGATTATCTGGGCCCTGTTACAGAGCTATCAGGTTTACCAACGCAACCCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAAt                       <  4:203182/148‑1 (MQ=11)                                                                                          tCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAAt                       <  6:271764/148‑1 (MQ=32)                                                                                          tCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAAt                       <  8:330458/148‑1 (MQ=32)                                                                                            ccGCCTTGGGCGCAAGGAGGTGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAAt                       <  4:203186/146‑1 (MQ=17)                                                                                            ccGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGc                     >  5:64912/1‑148 (MQ=32)                                                                                            ccGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGc                     >  3:292619/1‑148 (MQ=32)                                                                                                cTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCCACCGGGTAATGCCACCg                >  5:384979/1‑149 (MQ=17)                                                                                                   gggCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTaa                        <  5:365390/138‑1 (MQ=32)                                                                                                   gggCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTaa                        >  6:365390/1‑138 (MQ=32)                                                                                                      cgcAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGgcgc          <  6:392126/149‑1 (MQ=32)                                                                                                      cgcAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGgcgc          >  8:123039/1‑149 (MQ=32)                                                                                                      cgcAAGGAGGCGGCAATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGgcgc          >  8:50150/1‑149 (MQ=17)                                                                                                            gaggCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGgcg    <  1:196652/149‑1 (MQ=32)                                                                                                              ggcggcGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCt  <  4:292619/149‑1 (MQ=32)                                                                                                              ggcggcGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCt  <  1:363857/149‑1 (MQ=32)                                                                                                                |                                                                                                                                                   CCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAG‑CCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCT  >  NC_002947/401251‑401507

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