Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,522,277	A→G	20.3%	M84V (ATG→GTG)	yrfG →	GMP/IMP nucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,522,277	0	A	G	20.3%	80.7 / 2.7	51	M84V (ATG→GTG)	yrfG	GMP/IMP nucleotidase
Reads supporting (aligned to +/- strand):  ref base A (28/9);  new base G (2/10);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 4.61e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGAGCATACCGGTC  >  CP009273/3522129‑3522421                                                                                                                                                     |                                                                                                                                                 ggTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGa                                                                                                                                                  <  2:736967/149‑1 (MQ=255)  gTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCCATATCACGACGTACAGCATACGCTAACCTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGAt                                                                                                                                                 >  1:932019/1‑149 (MQ=255)   tcgtcgAAACATGGGGCGCGAAAAGCGGGGTTACGCCACAGGAAGCGATTGAATAGATGTGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATg                                                                                                                                                <  2:674081/144‑1 (MQ=255)     ccTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGGCTGGATATCTGTGCGATGac                                                                                                                                              >  2:50913/1‑149 (MQ=255)     ccTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGGGAGCAACTGGGTCTGGATATCTGTGCGGTGac                                                                                                                                              >  1:187899/1‑149 (MQ=255)             aTGGGGCGGGAAAAACGGCGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCACCTGGGTCTGGATATCTGTGCGATGACCACCgaga                                                                                                                                      <  1:136312/149‑1 (MQ=255)                   cttgAAAAACGAGTTTACGTTACAGGAAGCAAGTGAATATCTTCGCCAGCAATATCACGACGTACAGCATACGCTAATCTAGTACTGTCTTGATTACTGGAGTGAGCAACTTTTTTTTTTTATTGGTGCGGCGACCACCGAGATGGGAc                                                                                                                                <  2:698510/146‑1 (MQ=255)                     cgAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATAAGCTAAACTGGTACGGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACGACCGAGATGGGACcg                                                                                                                              >  1:598699/1‑149 (MQ=255)                      gAAAACCGGGTTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTTGTACTGTCTTGATTACTGGAGTGAGCAACTTGGTCTTGTTAATGATACGGCGACCACCGAGATGGGACcgc                                                                                                                             <  2:918674/149‑1 (MQ=255)                       aaaaaCGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACGGCATACGCTAAACTGGTCCTGTCTCGATAACTGGAGTGGGCAACTGGGTCTGGAACTCTGTCCGACGACCACCGAGATGGGACagcg                                                                                                                            >  2:219241/1‑149 (MQ=255)                          aaCGGAGTTACGCCACTGGAACCGGTGTAATATAGGCGCCAGCAATATCACGACGTACAGCTTACGCTAAACTTGTACTGTCTTGATTACTTGAGTGAGCAACTTGTTTTTTTTAATGATGCGGCGACCACCGAGATGGGACCGCGTGc                                                                                                                         <  2:1138829/149‑1 (MQ=255)                              gggTTACGCCACAGGAAGCGATTGGACATATGGGCCAGCAATAGCACGACGTACATCATACGCGAAACTGGTACTTGCTTGATTACTGGAGTGGGCAACTGTGTCTTGTTAATGATGCGGCGACCACCGAGATGGGACCGCGTGCCGTa                                                                                                                     <  2:137921/149‑1 (MQ=255)                              gggTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCGGTa                                                                                                                     >  1:462817/1‑149 (MQ=255)                                 ttACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTg                                                                                                                  >  1:878356/1‑149 (MQ=255)                                 ttACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGGCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGGGGCGTACtt                                                                                                                  >  2:149369/1‑148 (MQ=255)                                         caGGAAGCGATGGAATATATGCGCCAGCGATATCACGACGTACAGCATACACCAAACTGGTACGGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATGTCTGTGCCATGACCACGCAGGTCGGGCCGCGTGCCGTGCTGCGTGagca                                                                                                          >  2:698640/1‑146 (MQ=255)                                                 gATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTc                                                                                                  >  1:719135/1‑149 (MQ=255)                                                 gATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATAACATTc                                                                                                  >  1:650641/1‑149 (MQ=255)                                                   tGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTCCTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGCTGACCACCGAGATGGGACCACGTGCCGCGCTGTGTGAAGATCCCATTctg                                                                                                >  1:90917/1‑147 (MQ=255)                                                            tGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGTTTTGTTTAATGGTGCGGCGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAg                                                                                       <  2:905280/149‑1 (MQ=255)                                                             gcgcCGGCAATATCACGACGTACAGTGTGGGTTTAACTGGAGCTGTCTTGATTACTGGGGTGAGCAACTGGGTCTGGATATCTGTGCGGTGACCACCGGGATGGGACCGCGTGCCGTAGTGGCTGTAGAGACACTTGGGTTTATTGAgg                                                                                      >  2:211655/1‑149 (MQ=255)                                                             gcgcCAGCAATCTCACGACGTACAGCATACGCTAAACTGGTACTGTCTCGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGCCCACCGCCGGGGGACCGCGAGCCCTCCGACGCGAAGACACCACTCCGttcttttagt                                                                                      >  2:187369/1‑141 (MQ=255)                                                             gcgcCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGAGGGGACCGCGGGCCGTACTGCGTGAAGACACCATTCCGTTTCTTGacg                                                                                      >  2:1178612/1‑147 (MQ=255)                                                                   gCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGCTGGGACCGCGTGGCGTACTGCGTGGAGATAGCATTCCGTTTCTTGAGGCCCTgg                                                                                >  1:404519/1‑148 (MQ=255)                                                                   gCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGa                                                                                >  1:1020878/1‑149 (MQ=255)                                                                       tatCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTTCGATGACCACCGAGACGGGACCGCGTGCCGTACTGCGAGAAGATACCAGTCCGCATCTTTAGGCAcataacgc                                                                            >  2:356104/1‑142 (MQ=255)                                                                            cgacgTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCg                                                                       <  1:945268/149‑1 (MQ=255)                                                                               cgTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTa                                                                    >  1:705207/1‑149 (MQ=255)                                                                                gTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATAACGATGCGGTGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTaa                                                                   <  2:694242/149‑1 (MQ=255)                                                                                  aCAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTGagc                                                                 >  1:355918/1‑149 (MQ=255)                                                                                           gCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCCCCGAAAGCCAGCGGTAAGCAGCGAAttt                                                        >  1:454447/1‑149 (MQ=255)                                                                                             tAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCACTTCTTTGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTg                                                      <  1:961413/149‑1 (MQ=255)                                                                                                          tCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGGGATGGGACCGCGTGCCGTTCTCCGTGCAGAAACCATTCCGCTTCTCGCGGGAACGAGAGCCTGCGGGTCGCAGCGAACTCTGCGCACAATTgctc                                         >  2:588533/1‑147 (MQ=255)                                                                                                               atcaCTGCAGTGAGCAACTGGGTCTGGATACCTGTACGAGGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCg                                    <  2:38160/146‑1 (MQ=255)                                                                                                                      gAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAAcc                             >  1:800096/1‑149 (MQ=255)                                                                                                                        gtgtGCAACTGGGTCTGGATATCTGTGCGGTGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGCAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTg                           <  2:663778/149‑1 (MQ=255)                                                                                                                         tGAGCAACTGGGTCTGGATATCTGTGCGATGACCACGGAGATGGGACCGCGTGCCGTACTGCGTGAAGACACCATTCCGTCTCTTGAGGCACTGAAAGACAGCGGTAAGCAGCGAATCTAGCTCACCAAAGCGCAGCAGGACAACCCgg                          >  1:255196/1‑149 (MQ=255)                                                                                                                         tGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGGTGGGACCGCGTGCCGTACTGTGTGAAGGTACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCGGCGAATCTTGCTCACCAATGCGCATCCGCCCAACCTgg                          >  1:545667/1‑149 (MQ=255)                                                                                                                          gAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGCGAAGGTACCATTCCGTTTCTTGAGCCACTAAAAGCCCGGGGTACGGCGCGAATTTTGGTCTCGAAGGCGAATCCggaacacgtggc                         >  2:426952/1‑138 (MQ=255)                                                                                                                           ggcacgtGGGTCTGGTTATCTGTGCGGTGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCg                        <  2:131303/143‑1 (MQ=255)                                                                                                                           aGCAACTGGGTCTGGTTAACTGTGCGGTGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCg                        <  2:1147413/149‑1 (MQ=255)                                                                                                                           aGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCGGCGGTAAGGAGGGAATTTTGCTCACCAGTGCGCATCGGCCGAACCTGGCg                        >  1:190951/1‑149 (MQ=255)                                                                                                                                    gTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAGCtt               >  1:206997/1‑149 (MQ=255)                                                                                                                                      cttgttatCTGTGCGGCGACCACCGAGATGGGACCGCGTGCCGTACTGCTTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGGAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGa             <  2:812644/144‑1 (MQ=255)                                                                                                                                      cttgttatCGGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGa             <  2:702786/144‑1 (MQ=255)                                                                                                                                      ctgtatagtggtgCGACGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGa             <  2:998729/139‑1 (MQ=255)                                                                                                                                      cTGGATATCTGTGCGATGACCTCCGAGTTGGGGCCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGa             <  1:333186/149‑1 (MQ=255)                                                                                                                                      cTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGa             <  1:441176/149‑1 (MQ=255)                                                                                                                                       tgtatatCGGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGAg            <  2:1006157/146‑1 (MQ=255)                                                                                                                                       tGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGAg            >  1:86374/1‑149 (MQ=255)                                                                                                                                                 tgCGATGACCACCGAGATGGGACCGCGTGCCGTACGGGGTGAAGATACCAATCGGTATGCAGCGAAACTGAAAGCCAAAGGTGAGCAGGGAATTTTGGTCACTAATACGCATCCGCACAACGCGGGGGTCAAACCTGAGCTTATCggct  >  2:1132332/1‑147 (MQ=255)                                                                                                                                                 tgCGATGACCACCGAGATGGGACCGCGCGCCGTACTGCGTGCAGAAACCCATCCGGCACCTGAGGCACTGAAAGCCAGCGGCAAGCAGCGAATTTTGCGCACCAATGCGCATCCGCAAAAGCGGGCGGCAAAACGTGACCATAACGGtt  >  2:793705/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACTGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGAGCATACCGGTC  >  CP009273/3522129‑3522421

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