Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,378	C→T	14.7%	L58L (TTG→TTA)	tpn ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,378	0	C	T	14.7%	110.7 / 12.8	54	L58L (TTG→TTA)	tpn	IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base C (25/20);  new base T (4/4);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01

CGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAC  >  CP000730/36238‑36515                                                                                                                                             |                                                                                                                                          cGCAAACAAACATGCTATGTATGTCCCTCTGCATCCATGGCACGATATAAATCGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTACGCTTTTTCTTCCAAATTTGATAc                                                                                                                                           <  2:66623/141‑1 (MQ=14)   cAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAATTTGGATACaa                                                                                                                                         >  2:359524/1‑141 (MQ=255)     aCCAAATATCTAATGGATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaaaa                                                                                                                                       <  2:460494/141‑1 (MQ=255)          atatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCGCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGATTTTTCTTCCAAATTTGATACAAAATTggg                                                                                                                                  >  2:152274/1‑141 (MQ=255)           tatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTgggg                                                                                                                                 <  2:413688/141‑1 (MQ=255)           tatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTgggg                                                                                                                                 <  2:348516/141‑1 (MQ=255)           tatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTATTCTTCCAAATTTGATACAACATTgggg                                                                                                                                 >  2:33162/1‑141 (MQ=255)                     atgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAACATTggggacgatccgtc                                                                                                                       >  1:419953/1‑131 (MQ=255)                      tgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGa                                                                                                                      <  2:387828/141‑1 (MQ=255)                      tgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGa                                                                                                                      <  2:387833/141‑1 (MQ=255)                        tCCCTCTTCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACTTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACACAATTTGGGCGTATTTTTGGAc                                                                                                                    >  1:45360/1‑141 (MQ=14)                        tCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCACATTTGGTATCAAATTGTGGAGTATGCTTTAAc                                                                                                                    >  1:379607/1‑141 (MQ=11)                         cccTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGACTATTCTTTAAcc                                                                                                                   >  2:85203/1‑141 (MQ=255)                          ccTCTGCATCAATGGCACGATATACATAACTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAAccc                                                                                                                  <  2:301031/141‑1 (MQ=255)                          ccTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAAccc                                                                                                                  >  2:178699/1‑141 (MQ=40)                             ctGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCTAAc                                                                                                               >  2:251818/1‑141 (MQ=255)                                 aTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTa                                                                                                           <  2:394022/141‑1 (MQ=18)                                  tCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGATTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAg                                                                                                          >  1:183236/1‑141 (MQ=255)                                         cACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTg                                                                                                   >  1:181086/1‑141 (MQ=39)                                                  aaTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatgaaaac                                                                                          >  2:461892/1‑137 (MQ=255)                                                            ttttCCTTTTATTTTGATGTACGTCTGATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGACAAAAAATTGGGGAATATATTTTAACCCAACGGTAAACCGTTGAATGATGAACCATAAcacc                                                                                >  2:287151/1‑141 (MQ=11)                                                             tttCCTTTTATTTTGTTGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacca                                                                               <  2:361280/141‑1 (MQ=255)                                                              ttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcaccac                                                                              <  1:2304/141‑1 (MQ=39)                                                               tCCTCTTACTTTGATGTACGTCCCATCAATACGCTATTTGTAACAAGCTTTATTCTGCTTTTTCTTGCAAATTTGACACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACg                                                                             <  1:324436/141‑1 (MQ=11)                                                                 cTTTTATTTTGATATACGTTTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt                                                                           <  1:52349/141‑1 (MQ=255)                                                                  ttttATTTTGATGGACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTCGTCCGTTGAATGATGAACGTTTACACCACGTTc                                                                          <  2:14616/141‑1 (MQ=255)                                                                    ttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTccc                                                                        <  1:214765/141‑1 (MQ=39)                                                                        tttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTaa                                                                    <  2:348159/141‑1 (MQ=40)                                                                           gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAatct                                                                 >  2:12353/1‑139 (MQ=255)                                                                           gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAatat                                                                 >  2:55760/1‑141 (MQ=40)                                                                                     ttcgcaATGCGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGatata                                                       <  1:83435/137‑1 (MQ=255)                                                                                     tcatcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatata                                                       >  1:59382/1‑141 (MQ=255)                                                                                        tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAAAGATTAAAGCTCGCGCTCTTAATATTTCAGATATATCa                                                    >  1:406835/1‑141 (MQ=11)                                                                                           aTACGGCAGTTGTAATAAGCTTTTTTATGCTTGTTCTTCCAAATTTGGTACAAAATTGGGGCATATTCTTGAACCCGACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGa                                                 <  2:11958/141‑1 (MQ=14)                                                                                             aCGCCATTTGTAATAAGCCTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCTTATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACTTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATa                                               <  2:268284/141‑1 (MQ=11)                                                                                              cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                              >  1:345551/1‑141 (MQ=40)                                                                                                ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatatacccccagagcc                                            >  2:249490/1‑130 (MQ=11)                                                                                                ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCATGATAACt                                            >  2:368227/1‑141 (MQ=255)                                                                                                ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAcc                                            >  2:66669/1‑140 (MQ=255)                                                                                                 cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                           >  1:188629/1‑141 (MQ=40)                                                                                                   tttGTAATAAGCTTTTTTAGGCTTTTGCTTGCAAATGTGATACGAAATTGGGGCATATTCTTGAACCCAACGGGAGACCGATGAATGATGAACGGGTTCACCAAGTTCCCTTAATATTTCAGCTAAATCagggatacacaa                                         >  2:231197/1‑130 (MQ=14)                                                                                                    ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                        <  1:372056/141‑1 (MQ=40)                                                                                                     tGTAAGAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATCCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                       <  1:439835/141‑1 (MQ=255)                                                                                                         ataaGCTTTTTCATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCCATGtatt                                   >  2:332872/1‑140 (MQ=255)                                                                                                            aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCt                                <  2:345551/141‑1 (MQ=38)                                                                                                                     aTGCTTTTTCTTCCACATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                       >  1:342785/1‑141 (MQ=25)                                                                                                                      tGCTTTTTCTCCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTCTGGAACTGGTCCTTAATATTTCAGATATATCACGATAACTCAATGTATATTTTAGAtagtag                      >  1:335742/1‑141 (MQ=11)                                                                                                                             tcttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATCGTAGCCAACgg               >  1:83020/1‑141 (MQ=25)                                                                                                                               ttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCCGATATATCACGATAACTCAATGTATATCTTAGATAGTAGccgttggct             <  1:66669/141‑8 (MQ=255)                                                                                                                                  taaaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACa          >  2:180122/2‑141 (MQ=21)                                                                                                                                   aaaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAg         >  1:114869/1‑141 (MQ=34)                                                                                                                                    aaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGt        <  2:112489/141‑1 (MQ=35)                                                                                                                                     aTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGta       >  1:83901/1‑140 (MQ=25)                                                                                                                                      tttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGa      <  2:15677/141‑1 (MQ=34)                                                                                                                                          aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc  <  2:250608/141‑1 (MQ=35)                                                                                                                                             |                                                                                                                                          CGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAC  >  CP000730/36238‑36515

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑