Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,805,353	T→C	E189E (GAA→GAG)	radC ←	DNA repair protein RadC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,805,353	0	T	C	95.3%	219.7 / ‑0.7	64	E189E (GAA→GAG)	radC	DNA repair protein RadC
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (29/32);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACAAAGATCCCGAAGGATCGCGCATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGC  >  NZ_CP009273/3805220‑3805488                                                                                                                                      |                                                                                                                                        aCAAAGATCCCGAAGGATCGCGCATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAg                                                                                                                                     >  1:938924/1‑138 (MQ=255) aCAAAGATCCCGAAGGATCGCGCATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAg                                                                                                                                     >  2:87125/1‑138 (MQ=255) aCAAAGATCCCGAAGGATCGCGCATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAg                                                                                                                                     >  1:821307/1‑138 (MQ=255)           cGAAGGATCGCGCATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTa                                                                                                                          >  2:6863/1‑139 (MQ=255)               ggATCGCGCATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCg                                                                                                                      <  1:325903/139‑1 (MQ=255)                 aTCGCGCATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCt                                                                                                                    >  2:188255/1‑139 (MQ=255)                   cgcgCATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGCGTTTATCCGCtt                                                                                                                   >  1:916753/1‑138 (MQ=255)                    gcgcATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTa                                                                                                                 <  1:432190/139‑1 (MQ=255)                      gcATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACt                                                                                                               <  1:1017357/139‑1 (MQ=255)                      gcATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACt                                                                                                               <  2:1010688/139‑1 (MQ=255)                      gcATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACt                                                                                                               <  1:34905/139‑1 (MQ=255)                              ggTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAg                                                                                                       <  1:792650/139‑1 (MQ=255)                               gTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAg                                                                                                       >  2:1150157/1‑138 (MQ=255)                                ttAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGca                                                                                                     >  2:1108918/1‑139 (MQ=255)                                  aaaTCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGcac                                                                                                    <  2:798523/138‑1 (MQ=255)                                  aaaTCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGcac                                                                                                    <  1:981579/138‑1 (MQ=255)                                  aaaTCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGcac                                                                                                    <  2:938924/138‑1 (MQ=255)                                       cAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCg                                                                                              <  1:407340/139‑1 (MQ=255)                                       cAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCg                                                                                              <  2:1133199/139‑1 (MQ=255)                                       cAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCg                                                                                              <  2:101759/139‑1 (MQ=255)                                       cAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCGTTACTGGGTTCAGCACAACCCg                                                                                              <  2:848045/139‑1 (MQ=255)                                         gCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCg                                                                                              >  2:148326/1‑137 (MQ=255)                                         gCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGaa                                                                                            <  1:301813/139‑1 (MQ=255)                                         gCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGaa                                                                                            <  1:188255/139‑1 (MQ=255)                                          ccGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAg                                                                                           >  1:1157250/1‑139 (MQ=255)                                               ttCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGa                                                                                      >  2:393718/1‑139 (MQ=255)                                                  ggCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGAtt                                                                                    <  2:821307/138‑1 (MQ=255)                                                  ggCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGAtt                                                                                    <  1:1150157/138‑1 (MQ=255)                                                  ggCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGAtt                                                                                    <  1:367687/138‑1 (MQ=255)                                                  ggCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGAtt                                                                                    <  2:816133/138‑1 (MQ=255)                                                                   ctcCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGa                                                                                                                               >  1:325203/1‑78 (MQ=255)                                                                   ctcCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGa                                                                                                                               <  2:325203/78‑1 (MQ=255)                                                                   ctcCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTAt                                                                                  >  2:631543/1‑123 (MQ=255)                                                                   ctcCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTAt                                                                                  <  1:631543/123‑1 (MQ=255)                                                                    tcCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGcc                                                                 >  2:902439/1‑139 (MQ=255)                                                                        cGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAgg                                                             >  2:477918/1‑139 (MQ=255)                                                                              aTCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTa                                                       >  2:572436/1‑139 (MQ=255)                                                                              aTCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTa                                                       >  2:963252/1‑139 (MQ=255)                                                                                         gTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGAtt                                                                                    >  2:996625/1‑99 (MQ=255)                                                                                         gTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGAtt                                                                                    <  1:996625/99‑1 (MQ=255)                                                                                             aGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTcgc                                        >  2:267973/1‑139 (MQ=255)                                                                                             aGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTcgc                                        >  2:859366/1‑139 (MQ=255)                                                                                               cACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAAt                                                                                                                                 <  1:556022/48‑1 (MQ=255)                                                                                               cACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAAt                                                                                                                                 >  2:556022/1‑48 (MQ=255)                                                                                               cACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTcgcg                                       >  2:918998/1‑138 (MQ=255)                                                                                               cACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTcgcg                                       >  1:761077/1‑138 (MQ=255)                                                                                                 cgcgTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCa                                                                             >  1:101019/1‑98 (MQ=255)                                                                                                 cgcgTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCa                                                                             <  2:101019/98‑1 (MQ=255)                                                                                                 cgcgTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGa                                     >  2:803428/1‑138 (MQ=255)                                                                                                 cgcgTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGAt                                    <  1:477918/139‑1 (MQ=255)                                                                                                    gTAAATCCATGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAAt                                 >  1:1136592/1‑139 (MQ=255)                                                                                                             tGAACTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAg                         >  2:613547/1‑138 (MQ=255)                                                                                                                aCTGACAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATg                     <  1:749709/139‑1 (MQ=255)                                                                                                                    aCAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGAc                  <  1:286397/138‑1 (MQ=255)                                                                                                                    aCAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGAc                  <  1:963252/138‑1 (MQ=255)                                                                                                                    aCAACTCTTTATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACt                 >  1:218756/1‑139 (MQ=255)                                                                                                                           tttATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAAc            <  1:902439/137‑1 (MQ=255)                                                                                                                           tttATTATCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAAc            <  1:1013260/137‑1 (MQ=255)                                                                                                                             tattatCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACa           <  1:393718/136‑1 (MQ=255)                                                                                                                             tattatCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGCTGATTATGTGCAAGGATCATCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACa           <  2:916753/136‑1 (MQ=255)                                                                                                                                  tCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGAt                      <  2:1033565/120‑1 (MQ=255)                                                                                                                                  tCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGAt                      >  1:1033565/1‑120 (MQ=255)                                                                                                                                  tCCGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAg   >  2:361245/1‑139 (MQ=255)                                                                                                                                   ccGCTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGc  >  2:1114766/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                        ACAAAGATCCCGAAGGATCGCGCATAGCGGGTTAAATCCAGCCGCGTTCGGCAAAAGAAACATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGC  >  NZ_CP009273/3805220‑3805488

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