Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,873,583	A→G	S333S (AGT→AGC)	recF ←	DNA replication/repair protein RecF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,873,583	0	A	G	100.0%	221.7 / NA	62	S333S (AGT→AGC)	recF	DNA replication/repair protein RecF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/30);  total (32/30)

GATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTG  >  NZ_CP009273/3873450‑3873715                                                                                                                                      |                                                                                                                                     gATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGAt                                                                                                                                 >  1:835731/1‑139 (MQ=255)     cTGGAGGAGTCCTAAGCATTCGCCATCAACGTTTCTCGCTCATTTATATTTTGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATcgcg                                                                                                                             <  1:482184/139‑1 (MQ=255)        gaggagTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGATGATCGCGCt                                                                                                                           >  2:747085/1‑138 (MQ=255)           gagTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGAc                                                                                                                        >  2:711181/1‑138 (MQ=255)              tCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACa                                                                                                                       >  2:655645/1‑136 (MQ=255)                aTAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGAc                                                                                                                  >  1:129190/1‑139 (MQ=255)                aTAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGAc                                                                                                                  >  2:668503/1‑139 (MQ=255)                        tCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGAtt                                                                                                          >  1:673146/1‑139 (MQ=255)                          gACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGc                                                                                                        >  1:305139/1‑139 (MQ=255)                           aCATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCg                                                                                                       <  2:530289/139‑1 (MQ=255)                           aCATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCg                                                                                                       <  2:840732/139‑1 (MQ=255)                            cATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGt                                                                                                      <  2:643178/139‑1 (MQ=255)                             aTCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGt                                                                                                      <  2:403775/138‑1 (MQ=255)                               cAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGc                                                                                                   >  1:770986/1‑139 (MQ=255)                                       ctcGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAg                                                                                            <  2:410677/138‑1 (MQ=255)                                               tttATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGc                                                                                    >  2:102387/1‑138 (MQ=255)                                                  ataCTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAag                                                                                 <  2:478930/138‑1 (MQ=255)                                                  ataCTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAag                                                                                 <  1:655645/138‑1 (MQ=255)                                                   taCTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAagca                                                                               <  2:198845/139‑1 (MQ=255)                                                         ggTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCcg                                                                          <  1:191727/138‑1 (MQ=255)                                                          gTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCcgc                                                                         >  2:495344/1‑138 (MQ=255)                                                          gTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCcgc                                                                         >  2:34343/1‑138 (MQ=255)                                                           ttAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGa                                                                       <  1:164909/139‑1 (MQ=255)                                                             aaTCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACg                                                                     >  1:88099/1‑139 (MQ=255)                                                              aTCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACACAGAcc                                                                                                                 >  2:736288/1‑94 (MQ=255)                                                              aTCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGAcc                                                                                                                 <  1:736288/94‑1 (MQ=255)                                                               tCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCt                                                                   >  2:705985/1‑139 (MQ=255)                                                                    tATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAGGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCtcatca                                                              >  2:629241/1‑139 (MQ=255)                                                                             cTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCa                                                     <  1:34343/139‑1 (MQ=255)                                                                                ttCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCagag                                                  >  1:209239/1‑139 (MQ=255)                                                                                ttCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCagag                                                  >  2:76852/1‑139 (MQ=255)                                                                                          aaCATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAatcat                                        <  1:711181/139‑1 (MQ=255)                                                                                             aTCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTa                                     <  1:599173/139‑1 (MQ=255)                                                                                                ttCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                    >  2:602608/1‑137 (MQ=255)                                                                                                  cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCGTCAAGCTCAGAGGCAAAATCATCTAt                                    <  2:88099/135‑1 (MQ=255)                                                                                                  cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                    <  2:265070/135‑1 (MQ=255)                                                                                                  cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                    >  2:190108/1‑135 (MQ=255)                                                                                                  cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                    <  2:344748/135‑1 (MQ=255)                                                                                                  cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                    <  2:129190/135‑1 (MQ=255)                                                                                                  cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                    <  1:190108/135‑1 (MQ=255)                                                                                                  cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                    <  1:600505/135‑1 (MQ=255)                                                                                                  cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAgg                                >  1:117820/1‑139 (MQ=255)                                                                                                  cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAg                                 <  1:862674/138‑1 (MQ=255)                                                                                                   gAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGt                               <  2:117820/139‑1 (MQ=255)                                                                                                                           aCGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCAc       <  2:209239/139‑1 (MQ=255)                                                                                                                             gtgtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACgg     <  2:375745/139‑1 (MQ=255)                                                                                                                               gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGt                                                                                                      >  2:484869/1‑40 (MQ=255)                                                                                                                               gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGt                                                                                                      <  1:484869/40‑1 (MQ=255)                                                                                                                               gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                    <  2:520390/106‑1 (MQ=255)                                                                                                                               gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                    >  1:520390/1‑106 (MQ=255)                                                                                                                               gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGccgccc               >  2:759807/1‑127 (MQ=255)                                                                                                                               gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGccgccc               <  1:759807/127‑1 (MQ=255)                                                                                                                               gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGt   >  1:737662/1‑139 (MQ=255)                                                                                                                               gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGt   >  1:4673/1‑139 (MQ=255)                                                                                                                               gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGt   >  2:493934/1‑139 (MQ=255)                                                                                                                               gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGt   >  2:410281/1‑139 (MQ=255)                                                                                                                               gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGt   >  1:862109/1‑139 (MQ=255)                                                                                                                               gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGAAGAGACACCGCCGCCCGCTTTCACGGGt   >  2:573285/1‑139 (MQ=255)                                                                                                                                ttCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                    >  1:448439/1‑105 (MQ=255)                                                                                                                                ttCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                    <  2:448439/105‑1 (MQ=255)                                                                                                                                ttCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGt   >  1:691463/1‑138 (MQ=255)                                                                                                                                ttCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTg  <  2:862109/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     GATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTG  >  NZ_CP009273/3873450‑3873715

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