Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	607,884	G→A	R209* (CGA→TGA)	dck ←	deoxynucleoside kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	607,884	0	G	A	96.7%	182.4 / NA	60	R209* (CGA→TGA)	dck	deoxynucleoside kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (34/24);  minor base C (1/1);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGTATT  >  CP000730/607744‑608024                                                                                                                                             |                                                                                                                                             gtatataTCTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATATCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTGAAGCAAtt                                                                                                                                       >  2:578655/4‑148 (MQ=255)  tgtATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTTAATAATTCAAGCAAttt                                                                                                                                      >  1:62625/1‑148 (MQ=255)  tgtATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAttt                                                                                                                                      >  1:49999/1‑148 (MQ=255)     atatGTGAATTAATGTATTCGATAGAATTGATAAAATAATCAATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGt                                                                                                                                   <  1:50107/148‑1 (MQ=255)     atatGTGAATTAATGTAGTCGATAGAATTGATAAAATAATCGATAATATCTATGGACTGTACATGGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGt                                                                                                                                   <  1:299537/148‑1 (MQ=255)        tgtgAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCCTGTATCTCCTTGTCGATATATTTGAATAATTCAAGCAATTTTGTTTa                                                                                                                                >  1:588961/1‑148 (MQ=255)         gtgAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTAt                                                                                                                               >  2:331324/1‑148 (MQ=255)         gggAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTAt                                                                                                                               >  2:419596/3‑148 (MQ=255)          tgAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATc                                                                                                                              <  2:384324/148‑1 (MQ=255)                    tATTCGAGAGAAGTGATAAAAAAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGAATCTACGTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTa                                                                                                                    <  2:175147/148‑1 (MQ=255)                    caTTCGATAGAATTGATAAAATAATCTATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTa                                                                                                                    <  1:203822/147‑1 (MQ=255)                     aTTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTaa                                                                                                                   >  1:61037/1‑148 (MQ=255)                       tCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTCTGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCCATTTTGTTTATCATAGGATTTAAAg                                                                                                                 >  1:385228/1‑148 (MQ=255)                           tAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTc                                                                                                             >  2:81094/1‑148 (MQ=255)                              aattaattAAATAATCGATAATATCTATGGAATATACATCGAACATAAACATATACGCTCTGTCAATAGTCTTTTAGCGTGTATCTAATCGTCGATTGGTTTCAATAACATAAGCCATTTTGTTGTTCATAGGGTTTAAAGAGTCgcg                                                                                                          >  2:264628/1‑146 (MQ=255)                                 taaTAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTc                                                                                                       <  2:556194/146‑1 (MQ=255)                                  taTAAAAAAATCGATAATATCTATTGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCTATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTcc                                                                                                      <  2:364178/147‑1 (MQ=255)                                      aaaTAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTAt                                                                                                  <  2:525547/148‑1 (MQ=255)                                            tCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGatat                                                                                            >  1:579528/1‑148 (MQ=255)                                            tCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGatat                                                                                            >  1:617878/1‑148 (MQ=255)                                            tCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGGCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCCAGCAATTTTGTTTATCATAGGATTTAGAGATTAGGGGTCCTTATGAatgt                                                                                            >  1:109820/1‑146 (MQ=255)                                             cGATAATATCTATGGAATATACATTGAATATAAACATATACGATGTGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCCTAGGATTTAAAGATTCGGCGTCCTTATGGATTTc                                                                                           >  1:23832/1‑148 (MQ=255)                                                   tatCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTc                                                                                     <  1:576688/148‑1 (MQ=255)                                                     tCTATGGAAGATACATGGAAAATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTAGGGGTCCTTATGGATATCATATTCAt                                                                                   <  1:335000/148‑1 (MQ=255)                                                          tgAATATACAGTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAata                                                                              <  1:300145/147‑1 (MQ=255)                                                              tataCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGa                                                                          <  2:61037/148‑1 (MQ=255)                                                              tataCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGACTTAAAGATTAGGGGTCCTTATGGATATCATATTCATTAATATTGa                                                                          <  1:505337/148‑1 (MQ=255)                                                              tataCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAAGAATTACCGCCTTTCTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGa                                                                          <  1:301537/148‑1 (MQ=255)                                                              taaaCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGACTCGGGGTCCTTATGGATATCATATTCATTAATATTGa                                                                          <  1:12180/145‑1 (MQ=255)                                                                 aCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGATCTTATGGATATCATATTCATTAATATTCATAc                                                                       >  2:121094/1‑148 (MQ=255)                                                                  cATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGTATATCATATTCATTAATATTGATACg                                                                      <  1:598620/148‑1 (MQ=255)                                                                  cATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACg                                                                      <  1:529187/148‑1 (MQ=255)                                                                  cATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACg                                                                      >  1:476608/1‑148 (MQ=255)                                                                  cATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACg                                                                      >  1:215825/1‑148 (MQ=255)                                                                  cATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACg                                                                      >  1:366489/1‑148 (MQ=255)                                                                             taCAGATACGCTATGTCTTTTGTCTTTTATCGTTTTTCTACTTTTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATCAATATTGATACGTACAACTGGAc                                                           >  2:394651/2‑148 (MQ=255)                                                                             aaCATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGAc                                                           >  1:476605/1‑148 (MQ=255)                                                                                agataCGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATCTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATg                                                        <  1:419236/146‑1 (MQ=255)                                                                                  ataCGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCa                                                      >  1:577724/1‑148 (MQ=255)                                                                                  ataCGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATAAATATTGATACGTACAACAGGACATCCa                                                      >  2:604286/1‑148 (MQ=255)                                                                                   taCGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCAt                                                     >  1:3465/1‑148 (MQ=255)                                                                                   taCGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCAt                                                     >  1:556915/1‑148 (MQ=255)                                                                                    aCGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCAtt                                                    >  2:371978/1‑148 (MQ=255)                                                                                      gCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAGTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTaa                                                  >  1:421892/1‑148 (MQ=255)                                                                                        tATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAg                                                >  1:178532/1‑148 (MQ=255)                                                                                           gTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTTAAGCtc                                             >  2:424640/1‑147 (MQ=255)                                                                                            tCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTCAAGCTAt                                            <  2:6912/148‑1 (MQ=255)                                                                                                  gTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGATATTAATcc                                      >  2:166618/1‑148 (MQ=255)                                                                                                  gTCTTTTATCGTGTATCGATTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATcc                                      <  1:80707/148‑1 (MQ=255)                                                                                                    cTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCGTTAATATTGGTACGTACAACTGGACATGCATTATAGCTACTAATCCaa                                    >  2:244174/1‑148 (MQ=255)                                                                                                      tttATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATc                                  >  2:615721/1‑148 (MQ=255)                                                                                                      tttATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAAATGGACATGCATTAAAGCTATTAATCCAATc                                  >  2:604160/1‑148 (MQ=255)                                                                                                           cGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTCATCATAAGATTTAAAGATTCGGGGTACTTATGCATATCATATTCATTTATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCAt                             >  2:461030/1‑148 (MQ=255)                                                                                                                 tCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATTGATATCATATTCATTAATATTGATAGGTACAACCGGGTCAGTGGTAAAGCTATTAATCCGCTCGTCATAGCGtt                       <  1:463429/148‑1 (MQ=255)                                                                                                                      ttGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAAt                  >  1:493387/1‑148 (MQ=255)                                                                                                                          cGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCt              <  1:81094/148‑1 (MQ=255)                                                                                                                             tatGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCtttt           <  2:139193/148‑1 (MQ=255)                                                                                                                              atGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCttttt          <  2:49999/148‑1 (MQ=255)                                                                                                                               tGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTc         <  2:476605/148‑1 (MQ=255)                                                                                                                                      ataatTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGTAtt  >  1:314430/1‑148 (MQ=255)                                                                                                                                             |                                                                                                                                             GTGTATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGTATT  >  CP000730/607744‑608024

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