Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,610,344	T→C	W309R (TGG→CGG)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,610,344	0	T	C	81.8%	109.0 / 24.4	55	W309R (TGG→CGG)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (4/6);  new base C (31/14);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCT  >  NZ_CP009273/4610258‑4610421                                                                                       |                                                                              cgcTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTggg                                                                            >  2:561326/1‑90 (MQ=255) cgcTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCggg                                                                            >  2:608819/1‑90 (MQ=255) cgcTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCggg                                                                            >  2:110987/1‑90 (MQ=255)  gcTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCgggg                                                                           >  2:100437/1‑90 (MQ=255)  gcTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCgggg                                                                           >  1:712429/1‑90 (MQ=255)   cTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGc                                                                          <  2:724707/90‑1 (MQ=255)    tGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGcc                                                                         <  2:569383/90‑1 (MQ=255)      tttGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCAc                                                                       >  1:723197/1‑90 (MQ=255)               accaTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGa                                                               <  1:201104/89‑1 (MQ=255)               accaTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGa                                                               >  2:201104/1‑89 (MQ=255)               accaTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGAc                                                              >  2:935265/1‑90 (MQ=255)               accaTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGAc                                                              >  1:955560/1‑90 (MQ=255)                   tCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCg                                                          <  1:62211/90‑1 (MQ=255)                     aaaGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGcc                                                        <  1:362432/90‑1 (MQ=255)                          gaTGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTa                                                   >  2:428472/1‑90 (MQ=255)                          gaTGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTa                                                   >  2:143943/1‑90 (MQ=255)                          gaTGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTa                                                   >  2:589782/1‑90 (MQ=255)                                 gaagCGGGTGATACCACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCcgcg                                            <  1:471254/90‑1 (MQ=255)                                 gaagCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCcgcg                                            <  1:681671/90‑1 (MQ=255)                                  aagCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCgg                                           >  1:411775/1‑90 (MQ=255)                                    gCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGt                                         >  2:326868/1‑90 (MQ=255)                                    gCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGt                                         >  2:663777/1‑90 (MQ=255)                                    gCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGt                                         >  1:361064/1‑90 (MQ=255)                                     cGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTc                                        >  1:870254/1‑90 (MQ=255)                                         tGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGa                                    <  2:688979/90‑1 (MQ=255)                                               cacCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGtt                              >  1:406929/1‑90 (MQ=255)                                                 ccATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCcgcg                                            >  2:746513/1‑74 (MQ=255)                                                 ccATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCcgcg                                            <  1:746513/74‑1 (MQ=255)                                                 ccATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCg                            >  1:874867/1‑90 (MQ=255)                                                 ccATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCg                            >  2:604598/1‑90 (MQ=255)                                                 ccATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCg                            >  2:786571/1‑90 (MQ=255)                                                 ccATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCg                            >  2:253534/1‑90 (MQ=255)                                                 ccATCGTCTTCACCAACTTCGTCGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCg                            >  2:279428/1‑90 (MQ=255)                                                   aTCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCcgcg                                            <  2:721902/72‑1 (MQ=255)                                                   aTCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCcgcg                                            >  1:721902/1‑72 (MQ=255)                                                   aTCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTc                                        <  1:621197/76‑1 (MQ=255)                                                   aTCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTc                                        >  2:621197/1‑76 (MQ=255)                                                   aTCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGAc                          >  1:394092/1‑90 (MQ=255)                                                   aTCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGAc                          >  1:375333/1‑90 (MQ=255)                                                    tCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGAcc                         >  1:889746/1‑90 (MQ=255)                                                      gTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGAccgc                       <  2:467531/90‑1 (MQ=255)                                                      gTCTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGAccgc                       <  1:663902/90‑1 (MQ=255)                                                        cTTCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGAccgccg                     >  2:161448/1‑90 (MQ=255)                                                         ttCACCAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGt                    <  1:663777/90‑1 (MQ=255)                                                              cAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGcc               <  2:889746/90‑1 (MQ=255)                                                              cAACTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGcc               <  2:212820/90‑1 (MQ=255)                                                                aCTTCGTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCgg             >  1:441230/1‑90 (MQ=255)                                                                   tCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGc          <  2:874867/90‑1 (MQ=255)                                                                     gTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACt                                 >  2:155156/1‑65 (MQ=255)                                                                     gTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACt                                 <  1:155156/65‑1 (MQ=255)                                                                     gTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGc          >  1:895149/1‑88 (MQ=255)                                                                     gTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGc          <  2:895149/88‑1 (MQ=255)                                                                     gTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTg        >  2:695538/1‑90 (MQ=255)                                                                     gTTGACTTCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTg        <  2:870254/90‑1 (MQ=255)                                                                           ttCGACTCTTCCCGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGtct  >  2:63896/1‑90 (MQ=255)                                                                                       |                                                                              CGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCT  >  NZ_CP009273/4610258‑4610421

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