Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	199,508	A→C	18.6%	N336T (AAC→ACC)	U712_00905 →	Phosphoglucosamine mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	199,508	0	A	C	18.6%	86.2 / 2.0	48	N336T (AAC→ACC)	U712_00905	Phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base A (26/12);  new base C (2/7);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 2.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATAAGGCGCTCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAGTCAATGTGAGAGTGACTGAT  >  CP006881/199372‑199641                                                                                                                                         |                                                                                                                                      tataAGGCGCTCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTacaaca                                                                                                                                     >  7:16333/1‑139 (MQ=255) tataAGGCGCTCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTaaaaca                                                                                                                                     >  1:242062/1‑139 (MQ=255)          cTCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACACCGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGCTTACAACACGACAgggg                                                                                                                            >  6:9835/1‑139 (MQ=255)          cTCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACACCACGACAgggg                                                                                                                            >  8:181364/1‑139 (MQ=255)          cTCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACACCACGACAggg                                                                                                                             >  5:301991/1‑138 (MQ=255)          cTCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAgggg                                                                                                                            >  2:65557/1‑139 (MQ=255)          cTCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAggg                                                                                                                             >  8:191660/1‑138 (MQ=255)          cTCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGATGTCGGCGACCGCTACGTAGTAGAAGAAATGAAAAACGACGGCTACAACGTCGGCGCAGCGCTGTCAGGCCATCGTTTTTTCCTTGATTACAACACGCAAgggg                                                                                                                            >  5:47876/1‑139 (MQ=255)           tCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACCACACGACAGGGGa                                                                                                                           >  1:133502/1‑139 (MQ=255)           tCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGa                                                                                                                           >  5:180004/1‑139 (MQ=255)           tCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGa                                                                                                                           >  6:325362/1‑139 (MQ=255)           tCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGa                                                                                                                           >  5:327433/1‑139 (MQ=255)           tCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGa                                                                                                                           >  1:57326/1‑139 (MQ=255)                                   tGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCt                                                                                                   >  5:73792/1‑139 (MQ=255)                                   tGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCt                                                                                                   >  5:37866/1‑139 (MQ=255)                                       gACCGCTGGCGGCGACCGCAAGGTAGTAGAAGAAATCAAAAACCACGGCTACAACGTCGGCGGAGAGCAGGCAGGACATCTTATTTTCCTTGATTACACCACGACAGGGGACGGATTATTGTCTGCTATTATGCtgatg                                                                                               <  7:157531/139‑1 (MQ=255)                                           gCTGTCGGCGCCCGCTACGTAGTAGAAGCAATGAATAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAGCACGACAGGGGACGGATTATTGTCTGAGATTAGGCTGATGAaca                                                                                           >  5:51267/1‑139 (MQ=255)                                            cTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAacac                                                                                          >  8:233807/1‑139 (MQ=255)                                                      ccGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGc                                                                                 >  2:265178/1‑138 (MQ=255)                                                                      gCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGt                                                                >  3:180194/1‑139 (MQ=255)                                                                        aaacaaaaacCACGGCAACACCGGCGGCGGAGAGCAGTCAGGACATATTATTTTCATTGATAACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCa                                                              <  5:325362/135‑1 (MQ=255)                                                                         aTGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGt                                                                >  7:251475/1‑136 (MQ=255)                                                                         aTGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGt                                                                >  5:117215/1‑136 (MQ=255)                                                                         aTGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGt                                                                >  3:165816/1‑136 (MQ=255)                                                                          taaataaaGACGTCTCCATCGGCGGCGGAGAGCTGTCTGGACATCATCTTTTACTTCATTAAAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCaga                                                            <  6:123769/134‑1 (MQ=255)                                                                             aatcccaCGGCTACAAGGTGGGCGGATAGCATTCAGGACATATTTTTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCt                                                         <  4:305763/133‑1 (MQ=255)                                                                             aaaaccaCGTCTACAACTTCGGCGTAGAGCAGTCAGGACATCTTATTTTCCTTTATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCt                                                         <  4:180194/133‑1 (MQ=255)                                                                             aaaaccaCGGCAACAACGGCGGCGGAGACCAGTCAGGACTTCTTATTGCCCTTCATAACACCACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCt                                                         <  4:208459/133‑1 (MQ=255)                                                                                    cGGCTACAACGTCGGCGGAGAGCAGTCAGGACATGTTATTTTCCTTGATTACAACACGACAGGGGACGCATTATTGTCTGCTATTATGCTGATGAAAACATTAAAAGCAACAGGCAAGCCGATGTCAGAGATTGCAGCt                                                  >  6:211803/1‑139 (MQ=255)                                                                                     ggCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTg                                                 <  7:233807/139‑1 (MQ=255)                                                                                       cTACACGGTGGGCGGAGAGCAGTCAGGACATCTTATTTTCCTGGATTCCACCACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGaa                                               <  8:134517/139‑1 (MQ=255)                                                                                              gggggcggAGAGCAGGCAGGACATGTTATTTTCCTTGATAACACCACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGa                                        <  5:9835/136‑1 (MQ=255)                                                                                              gTCGGCGGAGAGCAGCCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGa                                        <  7:181364/139‑1 (MQ=255)                                                                                                 ggcggAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGt                                     >  5:163542/1‑139 (MQ=255)                                                                                                 ggcggAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACCACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGt                                     >  3:10128/1‑139 (MQ=255)                                                                                                 ggcggAGAGCAGTCAGGACATCTTATTTTCCTTGAGTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGAAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGt                                     >  1:17366/1‑139 (MQ=255)                                                                                                 ggcggAGAGCAGTCAGGACATCTTATTTTCATTGATTACAACAAGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGAAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGt                                     >  3:286443/1‑139 (MQ=255)                                                                                                      agagCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCg                                <  7:191660/139‑1 (MQ=255)                                                                                                       gagCAGTCAGGACATCTTTTTTTGCTTGATTACACCACGACTGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGc                               <  4:286443/139‑1 (MQ=255)                                                                                                                ggACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAg                      >  3:227543/1‑139 (MQ=255)                                                                                                                 gACATCTTATTTTCCTTGATTACACCACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAGt                     <  4:310655/139‑1 (MQ=255)                                                                                                                 gACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAGt                     <  6:73792/139‑1 (MQ=255)                                                                                                                 gACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAGt                     <  2:242062/139‑1 (MQ=255)                                                                                                                      cTTATTTTCCTTGATTACAACCCGCCAGGGGACGGATTATTTTCTGCTATTATTCTGATGAACACTTTAAAAGCAACAGGCAAGCAGCTTTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAGTAAAtg                >  7:288456/1‑139 (MQ=255)                                                                                                                      cTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATTCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAGTCAAtg                >  8:68555/1‑139 (MQ=255)                                                                                                                         aTTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAGTCAAt                 <  4:165816/135‑1 (MQ=255)                                                                                                                         aTTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAGTCAAt                 <  8:251475/135‑1 (MQ=255)                                                                                                                           tttCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAGTCAATGTga             >  4:303983/1‑137 (MQ=255)                                                                                                                               cTTGATAACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAGTCAATGTGAGAGTGa       <  3:128843/139‑1 (MQ=255)                                                                                                                                    ttACACCACGACAGGGGACGGATTATTGTCCGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAGTCAATGTGAGAGTGACTGAt  <  1:65557/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      TATAAGGCGCTCGAAAAAGAAGGCATCAAAAGCGTGCAGACAGCTGTCGGCGACCGCTACGTAGTAGAAGCAATGAAAAAAGACGGCTACAACGTCGGCGGAGAGCAGTCAGGACATCTTATTTTCCTTGATTACAACACGACAGGGGACGGATTATTGTCTGCTATTATGCTGATGAACACTTTAAAAGCAACAGGCAAGCCGCTGTCAGAGCTTGCAGCTGAAATGCAGAAGTTCCCGCAGCTGTTAGTCAATGTGAGAGTGACTGAT  >  CP006881/199372‑199641

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