Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,826,720	C→A	R12R (CGC→CGA)	xthA →	exodeoxyribonuclease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,826,720	0	C	A	100.0%	131.8 / NA	32	R12R (CGC→CGA)	xthA	exodeoxyribonuclease III
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (17/15);  total (17/15)

TTAACAACCATCGATCAAATCACTTAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAAC  >  NZ_CP009273/1826587‑1826843                                                                                                                                      |                                                                                                                            ttAACAACCATCGATCAAATCACTTAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAg                                                                                                                        >  2:587674/1‑139 (MQ=255)        ccATCGATCAAATCACTTAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCAc                                                                                                                 <  1:501152/139‑1 (MQ=255)         cATCGATCAAATCACTTAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCAcc                                                                                                                <  1:587674/139‑1 (MQ=255)         cATCGATCAAATCACTTAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCAcc                                                                                                                <  1:457060/139‑1 (MQ=255)                         tAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCAtcgt                                                                                                >  1:542420/1‑139 (MQ=255)                         tAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCAtcgt                                                                                                >  1:354744/1‑139 (MQ=255)                               caGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGaa                                                                                            >  1:681941/1‑137 (MQ=255)                                 ggCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAg                                                                                         >  1:144697/1‑138 (MQ=255)                                 ggCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAg                                                                                         >  1:171638/1‑138 (MQ=255)                                  gCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTgcgagc                                                                                                                           <  2:663788/103‑1 (MQ=255)                                  gCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTgcgagc                                                                                                                           >  1:663788/1‑103 (MQ=255)                                      tAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCaa                                                                                   >  2:566711/1‑139 (MQ=255)                                            aCGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGAt                                                                             >  1:143771/1‑139 (MQ=255)                                              gcgAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCAtcg                                                                                                 <  1:224940/117‑1 (MQ=255)                                              gcgAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCAtcg                                                                                                 >  2:224940/1‑117 (MQ=255)                                                      cTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGc                                                                    <  2:143771/138‑1 (MQ=255)                                                            ccATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCa                                                               <  2:354744/137‑1 (MQ=255)                                                              aTCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCAc                                                                                                                 <  1:171153/85‑1 (MQ=255)                                                              aTCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCAc                                                                                                                 >  2:171153/1‑85 (MQ=255)                                                                          tttATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGa                                                                         >  2:695475/1‑113 (MQ=255)                                                                          tttATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGa                                                                         <  1:695475/113‑1 (MQ=255)                                                                             aTCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGa                                             <  2:171638/138‑1 (MQ=255)                                                                             aTCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGa                                             <  2:681941/138‑1 (MQ=255)                                                                                  aaTAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGatat                                       >  2:364004/1‑139 (MQ=255)                                                                                  aaTAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGatat                                       >  1:410604/1‑139 (MQ=255)                                                                                  aaTAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGatat                                       >  1:251765/1‑139 (MQ=255)                                                                                  aaTAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGatat                                       >  2:672249/1‑139 (MQ=255)                                                                                    tAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGt                                     >  2:20540/1‑139 (MQ=255)                                                                                                         tttGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGa                <  1:20540/139‑1 (MQ=255)                                                                                                             tctcTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCt            <  1:655396/139‑1 (MQ=255)                                                                                                                 ttttAATATCAACGGCCTGCGAGCCAGACCTCAcc                                                                                                                >  2:563120/1‑35 (MQ=255)                                                                                                                 ttttAATATCAACGGCCTGCGAGCCAGACCTCAcc                                                                                                                <  1:563120/35‑1 (MQ=255)                                                                                                                      atatCAACGGCCTGCGAGCCAGACCTCTCCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACa    <  2:635178/138‑1 (MQ=255)                                                                                                                       tatCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAAc  <  2:683530/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                            TTAACAACCATCGATCAAATCACTTAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAAC  >  NZ_CP009273/1826587‑1826843

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