Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,609,448	A→G	D10G (GAC→GGC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,609,448	0	A	G	97.3%	228.3 / ‑1.2	73	D10G (GAC→GGC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (44/27);  total (45/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTG  >  NZ_CP009273/4609371‑4609532                                                                              |                                                                                     gTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCAt                                                                          >  2:97772/1‑90 (MQ=255)  tataCTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATc                                                                         <  1:97772/90‑1 (MQ=255)   ataCTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCg                                                                        >  2:428440/1‑90 (MQ=255)     aCTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgc                                                                      <  2:698503/90‑1 (MQ=255)     aCTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgc                                                                      <  2:509457/90‑1 (MQ=255)        gatctgatcATTTTAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  2:752289/1‑90 (MQ=255)        gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgcgc                                                                    <  2:554117/89‑1 (MQ=255)        gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgcgc                                                                    >  1:554117/1‑89 (MQ=255)        gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgc                                                                      <  2:462517/87‑1 (MQ=255)        gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgc                                                                      >  1:462517/1‑87 (MQ=255)        gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  1:1052259/1‑90 (MQ=255)        gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  2:997572/1‑90 (MQ=255)        gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  2:804278/1‑90 (MQ=255)        gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  2:1162141/1‑90 (MQ=255)        gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  1:95411/1‑90 (MQ=255)        gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  1:426074/1‑90 (MQ=255)        gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  1:1191485/1‑90 (MQ=255)        gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCt                                                                   >  1:122633/1‑90 (MQ=255)           ctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACa                                                                >  2:474472/1‑90 (MQ=255)             gatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACaga                                                              >  2:1133219/1‑90 (MQ=255)             gatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACaga                                                              >  1:780312/1‑90 (MQ=255)             gatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACaga                                                              >  2:368711/1‑90 (MQ=255)              atcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACagaa                                                             <  2:495740/90‑1 (MQ=255)                  tttAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATg                                                         >  1:72407/1‑90 (MQ=255)                  tttAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATg                                                         >  1:1286023/1‑90 (MQ=255)                     aaaTTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAg                                                      <  1:908635/90‑1 (MQ=255)                           gAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCt                                                >  1:390620/1‑90 (MQ=255)                                         ggAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCggg                                  <  1:1278715/90‑1 (MQ=255)                                         ggAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCggg                                  <  2:1278387/90‑1 (MQ=255)                                         ggAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCggg                                  <  2:126436/90‑1 (MQ=255)                                          gagaACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACagaa                                                             >  1:897588/1‑62 (MQ=255)                                          gagaACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACagaa                                                             <  2:897588/62‑1 (MQ=255)                                          gagaACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAGGATGCAGAACGCTTTGGTGACGTCgggg                                 >  1:945991/1‑90 (MQ=255)                                            gaACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCt                               >  1:1042890/1‑90 (MQ=255)                                             aaCATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTg                              >  1:910538/1‑90 (MQ=255)                                               cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgcg                                                                     <  1:7353/49‑1 (MQ=255)                                               cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgcg                                                                     >  2:7353/1‑49 (MQ=255)                                               cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCa                                                       <  1:144685/63‑1 (MQ=255)                                               cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCa                                                       >  2:144685/1‑63 (MQ=255)                                               cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            <  1:1162141/90‑1 (MQ=255)                                               cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            >  2:166513/1‑90 (MQ=255)                                               cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            >  2:497239/1‑90 (MQ=255)                                               cATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            >  1:1253038/1‑90 (MQ=255)                                               cATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            <  2:1280198/90‑1 (MQ=255)                                               cATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            >  1:1280198/1‑90 (MQ=255)                                                 tatGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgc                                                                      >  2:1086092/1‑46 (MQ=255)                                                 tatGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGgc                                                                      <  1:1086092/46‑1 (MQ=255)                                                 tatGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGacac                          >  1:657221/1‑90 (MQ=255)                                                  atGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACAcc                         >  2:1270424/1‑90 (MQ=255)                                                   tGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACAccc                        <  2:72507/90‑1 (MQ=255)                                                   tGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACAccc                        <  1:166513/90‑1 (MQ=255)                                                   tGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACAccc                        <  2:159556/90‑1 (MQ=255)                                                     aaaCGTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTg                      <  1:497239/90‑1 (MQ=255)                                                         gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTggtca                  >  2:212256/1‑87 (MQ=255)                                                         gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTggggc                  >  2:602667/1‑88 (MQ=255)                                                         gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTc                  >  2:283989/1‑90 (MQ=255)                                                         gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTc                  >  2:209564/1‑90 (MQ=255)                                                         gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTc                  >  2:168117/1‑90 (MQ=255)                                                         gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTc                  >  1:820878/1‑90 (MQ=255)                                                         gTGCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTc                  >  1:885436/1‑90 (MQ=255)                                                           gCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCat                <  2:95411/90‑1 (MQ=255)                                                           gCATTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCat                <  2:1088454/90‑1 (MQ=255)                                                             aTTTATTATGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCatat              <  2:910652/90‑1 (MQ=255)                                                                  ttatGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAg         <  1:184144/90‑1 (MQ=255)                                                                  ttatGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAg         <  1:564264/90‑1 (MQ=255)                                                                  ttatGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAg         <  2:1191485/90‑1 (MQ=255)                                                                     tGGTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAg      >  2:347184/1‑90 (MQ=255)                                                                       gTGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCt    >  1:737807/1‑90 (MQ=255)                                                                        tGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCtt   >  1:555354/1‑90 (MQ=255)                                                                        tGCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCtt   <  1:428440/90‑1 (MQ=255)                                                                         gCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            >  2:657655/1‑64 (MQ=255)                                                                         gCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                            <  1:657655/64‑1 (MQ=255)                                                                         gCTGGGCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTtg  >  1:587755/1‑90 (MQ=255)                                                                              |                                                                                     GTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTG  >  NZ_CP009273/4609371‑4609532

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