Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,625,502	G→A	7.3%	intergenic (+4095/+469)	PP_5654 → / ← PP_4091	phosphotriesterase/transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,625,502	0	G	A	7.3%	156.4 / 2.7	55	intergenic (+4095/+469)	PP_5654/PP_4091	phosphotriesterase/transposase
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base A (2/2);  total (29/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 = 9.02e-01

CGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGG  >  NC_002947/4625364‑4625641                                                                                                                                            |                                                                                                                                            cgcgGCCGACAATCTGCTGTCCATCACATTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTTCTCAGCTAAACCAACATCGAGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAAtt                                                                                                                                   <  2:203652/149‑1 (MQ=18)  gcgGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGTCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAAttt                                                                                                                                  <  4:112287/149‑1 (MQ=255)        gACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCACCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTCCAGAc                                                                                                                            >  8:60965/1‑149 (MQ=255)           aaTCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGTGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAAttt                                                                                                                                  <  3:91073/140‑1 (MQ=255)           aaTCTGCTGTCCATCACCTTAACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCAAACGGCCAACTGCTCAGCGAAACCAACACCGCGGGCCTGCTGCAGTAACCCTA‑CGACGCGCTGGGCAAttt                                                                                                                                  >  4:91073/1‑140 (MQ=18)                gctgTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGtgc                                                                                                                    >  1:221866/1‑149 (MQ=255)                             ttCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAgcg                                                                                                       >  3:338389/1‑149 (MQ=255)                             ttCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCAAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAgcg                                                                                                       >  1:47309/1‑149 (MQ=255)                              tCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCACTgg                                                                                                                                         >  3:23117/1‑114 (MQ=255)                              tCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCACTgg                                                                                                                                         <  4:23117/114‑1 (MQ=255)                                  aGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAAc                                                                                                  <  6:43656/149‑1 (MQ=255)                                  aGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAAc                                                                                                  >  5:43656/1‑149 (MQ=255)                                   gACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTACTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACt                                                                                                 <  1:342584/149‑1 (MQ=255)                                       aCCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAAt                                                                                             >  6:188076/1‑149 (MQ=255)                                        ccAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACCGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAgc                                                                                                        <  1:70123/138‑1 (MQ=255)                                        ccAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAgc                                                                                                        >  2:70123/1‑137 (MQ=255)                                                   atcagcaACTCGATAACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCACTGGGCAATTTGCAGACCCTGGTGCTGccc                                                                                                               <  4:254374/117‑1 (MQ=255)                                                   agcagcAACTCGATTACACATAAGACACCAACGTGCAACTGCTCAGCGAAACCAACAGAGCGGGCCTGCTGCAGTACCGCTA‑CGACGCACTGGGCAATTTGCAGACCCTGGtgctgcac                                                                                                               >  3:254374/1‑117 (MQ=17)                                                    gcagcaACTCGATTACACCTACGACCCCAACGGGCAGCTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACg                                                                                >  1:240672/1‑149 (MQ=255)                                                    gcagcaACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACg                                                                                >  6:148091/1‑149 (MQ=255)                                                    gcagcaACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACg                                                                                >  5:296857/1‑149 (MQ=255)                                                     cagcaACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACgg                                                                               <  5:148091/149‑1 (MQ=255)                                                       gcaACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACgg                                                                               >  5:108434/1‑147 (MQ=255)                                                          aCTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCg                                                                          <  5:188076/149‑1 (MQ=255)                                                              gATTACACCTACGACCCCAACGGGCGACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGaca                                                                      >  4:200544/1‑149 (MQ=255)                                                                   cacCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTg                                                                  <  2:240672/148‑1 (MQ=255)                                                                      cTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCAcc                                                              <  5:238713/149‑1 (MQ=255)                                                                      cTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCAc                                                               >  6:238713/1‑148 (MQ=255)                                                                                aaCGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACtt                                                    <  7:60965/149‑1 (MQ=255)                                                                                         cTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCg                                           >  3:35461/1‑149 (MQ=255)                                                                                          tGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCtac                                                                 >  6:298578/1‑124 (MQ=255)                                                                                          tGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCtac                                                                 <  5:298578/126‑3 (MQ=255)                                                                                          tGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAAc                                                <  6:45163/143‑1 (MQ=255)                                                                                          tGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAAc                                                >  5:45163/1‑143 (MQ=255)                                                                                                gCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGAt                                     >  3:149744/1‑148 (MQ=255)                                                                                                 cGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGc                                                                           <  7:105351/109‑1 (MQ=255)                                                                                                 cGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGc                                                                           >  8:105351/1‑109 (MQ=255)                                                                                                 cGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGAt                                     >  5:196710/1‑147 (MQ=255)                                                                                                   aaaCCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGc                                 >  5:10369/1‑149 (MQ=255)                                                                                                      ccAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGAt                              <  2:221866/149‑1 (MQ=255)                                                                                                        aaCAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGAttt                            <  6:108434/149‑1 (MQ=255)                                                                                                        aaCAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGAttt                            >  5:259260/1‑149 (MQ=255)                                                                                                         aCAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTc                           <  6:196710/149‑1 (MQ=255)                                                                                                         aCAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTc                           >  3:24660/1‑149 (MQ=255)                                                                                                                gggCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCt                                                                                        <  6:224736/81‑1 (MQ=255)                                                                                                                       ctgcAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTTGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGt             <  6:259260/149‑1 (MQ=255)                                                                                                                       ctgcAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGt                                        >  1:175492/1‑122 (MQ=255)                                                                                                                       ctgcAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGt                                        <  2:175492/122‑1 (MQ=255)                                                                                                                       ctgcAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCg               <  4:149744/147‑1 (MQ=255)                                                                                                                          cAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCa          >  4:246344/1‑149 (MQ=255)                                                                                                                          cAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCa          >  8:264738/1‑149 (MQ=255)                                                                                                                           aGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCac         >  2:107911/1‑149 (MQ=255)                                                                                                                           aGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCac         <  6:10369/149‑1 (MQ=255)                                                                                                                            gTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCacg        <  7:264738/149‑1 (MQ=255)                                                                                                                               ccGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCacgac      <  6:296857/148‑1 (MQ=255)                                                                                                                                           cgcTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAgg  >  1:166744/1‑141 (MQ=255)                                                                                                                                           cgcTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAgg  <  2:166744/141‑1 (MQ=255)                                                                                                                                            |                                                                                                                                            CGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTA‑CGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGG  >  NC_002947/4625364‑4625641

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