Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	3,480,475	T→A	6.9%	T49S (ACA→TCA)	U712_17865 ←	Amidase enhancer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,480,475	0	T	A	6.9%	158.4 / 3.6	58	T49S (ACA→TCA)	U712_17865	Amidase enhancer
Reads supporting (aligned to +/- strand):  ref base T (29/24);  new base A (2/2);  total (32/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 = 5.51e-01

TATCTCGGTTCACCAATATCACTGTATTAGGATTTTTCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGAACACACTATCAATTGTTTGCA  >  CP006881/3480337‑3480610                                                                                                                                           |                                                                                                                                        tatCTCGGTTCACCAATATCACTGTATTAGGATTTTTCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGt                                                                                                                                         <  1:130621/139‑1 (MQ=255) tatCTCGGTTCACCAATATCACTGTATTAGGATTTTTCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGt                                                                                                                                         <  3:258874/139‑1 (MQ=255)     tcGGTTCACCAATATCACTGTATTAGGATTTTTCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGc                                                                                                                                     >  3:197328/1‑139 (MQ=255)       ggTTCACCAATATCACTGTATTAGGATTTTTCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCg                                                                                                                                    <  1:133462/138‑1 (MQ=255)            aCCAATATCACTGTATTAGGATTTTTCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGAc                                                                                                                              <  6:62832/139‑1 (MQ=255)                 tatCACTGTATTAGGATTTTTCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATg                                                                                                                         <  1:77538/139‑1 (MQ=255)                     aCTGTATTAGGATTTTTCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGc                                                                                                                        >  7:11998/1‑136 (MQ=255)                        gTATTAGGATTTTTCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAGCGATCCGTACCAGCAACATTATCTCCAGTTTCTTTATAAAATCCTGGTGGCGCTAGACTGATGCTAGatt                                                                                                                  >  4:187082/1‑139 (MQ=255)                         tATTAGGATTTTTCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGattt                                                                                                                 >  4:194732/1‑139 (MQ=255)                           ttAGGATTTTTCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCAGTAACAGCAACATTATCTCCAGTTTCTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGatttat                                                                                                               >  3:85306/1‑139 (MQ=255)                                   tttCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCGAGTTACTTTATAAAATCCTGTTGGCGATTGACTGATGCTAGATTTATTTCCAtta                                                                                                       >  6:208502/1‑139 (MQ=255)                                   tttCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTTTTGGCGATAGACTGATGCTAGATTTATTTCCAata                                                                                                       >  3:62093/1‑139 (MQ=255)                                   tttCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGc                                                                                                                        <  8:22500/122‑1 (MQ=255)                                   tttCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGc                                                                                                                        >  7:22500/1‑122 (MQ=255)                                   tttCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAata                                                                                                       >  4:102715/1‑139 (MQ=255)                                   tttCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAata                                                                                                       >  7:51327/1‑139 (MQ=255)                                   tttCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAata                                                                                                       >  2:176475/1‑139 (MQ=255)                                     tCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAatat                                                                                                      >  1:144411/1‑138 (MQ=255)                                     tCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAatat                                                                                                      >  4:97518/1‑138 (MQ=255)                                         cTGGCTGTTGGATGCAAGCGGCGCAGTTTCGTCTCGTGCAGCGCCTGCAAAACGATCCGGAACATCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATaattt                                                                                                 <  2:256762/139‑1 (MQ=255)                                         cTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATaattt                                                                                                 <  2:256761/139‑1 (MQ=255)                                            gCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATaatttaa                                                                                               <  8:88715/138‑1 (MQ=255)                                            gCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATaatttaa                                                                                               <  8:11998/138‑1 (MQ=255)                                            gCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATaattta                                                                                                <  1:176475/137‑1 (MQ=255)                                                ttGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATaatttaacaatt                                                                                          <  7:280598/139‑1 (MQ=255)                                                   gATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATaatttaacaatttaa                                                                                       <  4:197328/139‑1 (MQ=255)                                                       cAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTg                                                                                    <  3:105586/138‑1 (MQ=255)                                                       cAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTg                                                                                    <  2:94777/138‑1 (MQ=255)                                                       cAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTg                                                                                    <  3:97518/138‑1 (MQ=255)                                                          gagTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGAt                                                                                >  6:69213/1‑139 (MQ=255)                                                             tCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGt                                                                                                                                         >  3:203287/1‑79 (MQ=255)                                                             tCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGt                                                                                                                                         <  4:203287/79‑1 (MQ=255)                                                             tCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGCCTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGtt                                                                             >  6:143865/1‑139 (MQ=255)                                                                             tGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCaa                                                               >  1:228411/1‑137 (MQ=255)                                                                             tGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCaa                                                               >  7:219475/1‑137 (MQ=255)                                                                                 gcgcCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGaaa                                                         <  5:69213/139‑1 (MQ=255)                                                                                       gCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGATTCTGCTGCAAAAGAAACTGATg                                                   >  2:11511/1‑139 (MQ=255)                                                                                           aaCGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAAt                                               <  2:191404/139‑1 (MQ=255)                                                                                           aaCGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAAt                                               <  2:144411/139‑1 (MQ=255)                                                                                               aTCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATc                                              <  4:116485/136‑1 (MQ=255)                                                                                               aTCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATc                                              >  3:116485/1‑136 (MQ=255)                                                                                               aTCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCaat                                           <  5:240964/139‑1 (MQ=255)                                                                                                       cAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGAt                                   <  3:162710/139‑1 (MQ=255)                                                                                                          cAACATTATCTCCAGTTACTTTATAAAATCCTGATGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTg                                 >  2:192451/1‑138 (MQ=255)                                                                                                          cAACATTATCTCCAGTTACTTTATAAAATCCTGATGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTg                                 >  2:188416/1‑138 (MQ=255)                                                                                                            aCATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTgcag                              >  5:235775/1‑139 (MQ=255)                                                                                                              aTTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAAAATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTgcagca                            <  5:208502/139‑1 (MQ=255)                                                                                                                    tcCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGaa                      >  6:79255/1‑139 (MQ=255)                                                                                                                    tcCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGaa                      >  4:72773/1‑139 (MQ=255)                                                                                                                      cAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTg                        <  3:102715/135‑1 (MQ=255)                                                                                                                      cAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGaa                      <  4:62093/137‑1 (MQ=255)                                                                                                                       aGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTTCAGCAAGTGAaca                    >  6:255059/1‑138 (MQ=255)                                                                                                                       aGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGAacac                   <  6:235775/139‑1 (MQ=255)                                                                                                                       aGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGAaca                    >  1:249098/1‑138 (MQ=255)                                                                                                                       aGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGAaca                    >  4:198845/1‑138 (MQ=255)                                                                                                                        gTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATATCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGAacaca                  >  4:217822/1‑139 (MQ=255)                                                                                                                            cTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGAACACACTAt              >  7:179035/1‑139 (MQ=255)                                                                                                                                       ccTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGAACACACTATCAATTGTTTg    >  1:51452/1‑138 (MQ=255)                                                                                                                                       ccTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGAACACACTATCAATTGTTTg    >  3:274676/1‑138 (MQ=255)                                                                                                                                        cTGATGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGAACACACTATCAATTGTTTGCa  <  1:192451/139‑1 (MQ=255)                                                                                                                                        cTGATGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGAACACACTATCAATTGTTTGCa  <  1:188416/139‑1 (MQ=255)                                                                                                                                         tGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGAACACACTATCAATTGTTTGc   >  1:69362/1‑137 (MQ=255)                                                                                                                                           |                                                                                                                                        TATCTCGGTTCACCAATATCACTGTATTAGGATTTTTCCACTGGCTGTTGGATGCAAGAGTCGCAGTTTCGTATCGTGAAGCGCCTGCAAAACGATCCGTAACAGCAACATTATCTCCAGTTACTTTATAAAATCCTGTTGGCGATAGACTGATGCTAGATTTATTTCCAATATAATTTAACAATTTAACTGAGATGTTTGAGTCTGCTGCAAAAGAAACTGATGGAATCAATAATAAGATTGCAGCAAGTGAACACACTATCAATTGTTTGCA  >  CP006881/3480337‑3480610

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