Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,882,139	A→G	E17G (GAG→GGG)	tnaA →	tryptophanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,882,139	0	A	G	98.4%	234.8 / NA	64	E17G (GAG→GGG)	tnaA	tryptophanase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (40/23);  total (41/23)
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.98e-01

TGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGC  >  NZ_CP009273/3882004‑3882270                                                                                                                                        |                                                                                                                                    tGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGc                                                                                                                                   >  2:531106/1‑138 (MQ=255)  gCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGc                                                                                                                                   <  2:616204/137‑1 (MQ=255)              ggCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTAc                                                                                                                      >  1:117987/1‑138 (MQ=255)               gctaaTTAATCGGTATAGCAGATGTAATATTAACCGGGATCACTGTAATTAAAATAAATGACGGATTATGGATTGGAAAACTTTAAACATCTCCCTGAACCGTTCGGCATTCGTGTTCTTGGGCCAGTATAACGTAcac                                                                                                                    >  2:547447/4‑137 (MQ=255)               gCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTAc                                                                                                                      <  1:531106/137‑1 (MQ=255)               gCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCa                                                                                                                    >  1:748932/1‑139 (MQ=255)                      aaTCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTcgcg                                                                                                              <  2:193612/138‑1 (MQ=255)                                     tGTAATATTCACAGGGATCACTGTAATTAATATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGc                                                                                              <  1:4885/139‑1 (MQ=255)                                      gTAATATTCACAGGGATCACTGTAATTAATATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGc                                                                                              >  2:4885/1‑138 (MQ=255)                                      gTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGc                                                                                              >  1:108364/1‑138 (MQ=255)                                      gTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGc                                                                                              >  2:248859/1‑138 (MQ=255)                                              cacaGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAAttatta                                                                                      >  2:330199/1‑138 (MQ=255)                                                 aGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATc                                                                                  >  1:395418/1‑139 (MQ=255)                                                 aGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATc                                                                                  >  1:783018/1‑139 (MQ=255)                                                      tCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTa                                                                             >  2:502176/1‑139 (MQ=255)                                                      tCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTa                                                                             >  1:641311/1‑139 (MQ=255)                                                        aCTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTAt                                                                            >  2:794199/1‑138 (MQ=255)                                                           gTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAAc                                                                        >  2:299348/1‑139 (MQ=255)                                                            tAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGaa                                                                         <  2:837667/137‑1 (MQ=255)                                                             aaTTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAAccc                                                                      >  2:522677/1‑139 (MQ=255)                                                              aTTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCg                                                                     >  1:449241/1‑139 (MQ=255)                                                               ttAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGt                                                                    >  1:149426/1‑139 (MQ=255)                                                               ttAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGt                                                                    <  2:748932/139‑1 (MQ=255)                                                               ttAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGt                                                                    >  2:247071/1‑139 (MQ=255)                                                               ttAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGt                                                                    >  2:794087/1‑139 (MQ=255)                                                               ttAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGt                                                                    >  1:247002/1‑139 (MQ=255)                                                                    ataaatGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCctg                                                               <  1:522677/139‑1 (MQ=255)                                                                                        tGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGtttt                                           >  2:191725/1‑139 (MQ=255)                                                                                        tGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGtttt                                           >  2:418216/1‑139 (MQ=255)                                                                                          gAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAattatt                                                                                       <  2:416796/93‑1 (MQ=255)                                                                                          gAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAattatt                                                                                       >  1:416796/1‑93 (MQ=255)                                                                                             aaCTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATc                                       >  1:668436/1‑138 (MQ=255)                                                                                             aaCTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATc                                       >  2:99383/1‑138 (MQ=255)                                                                                             aaCTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATc                                       >  1:309892/1‑138 (MQ=255)                                                                                                   aaaCATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTa                                 <  1:299348/138‑1 (MQ=255)                                                                                                            ccTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGAcc                           <  1:191725/135‑1 (MQ=255)                                                                                                            ccTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGAcc                           <  1:247071/135‑1 (MQ=255)                                                                                                            ccTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGAcc                           <  1:330199/135‑1 (MQ=255)                                                                                                             cTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAg                      <  2:668436/139‑1 (MQ=255)                                                                                                             cTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAg                      <  1:248859/139‑1 (MQ=255)                                                                                                             cTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAg                      <  1:438653/139‑1 (MQ=255)                                                                                                                    gTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCaa                                                                                            >  2:842826/1‑62 (MQ=255)                                                                                                                    gTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCaa                                                                                            <  1:842826/62‑1 (MQ=255)                                                                                                                    gTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTa                                                                             <  2:322499/77‑1 (MQ=255)                                                                                                                    gTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTa                                                                             >  1:322499/1‑77 (MQ=255)                                                                                                                    gTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCAcc               <  1:558213/139‑1 (MQ=255)                                                                                                                    gTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCAcc               >  2:558213/1‑139 (MQ=255)                                                                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATc                                                                                  >  1:520435/1‑71 (MQ=255)                                                                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATc                                                                                  <  2:520435/71‑1 (MQ=255)                                                                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTa                                 <  1:679860/120‑1 (MQ=255)                                                                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTa                                 >  2:679860/1‑120 (MQ=255)                                                                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCg              >  1:805609/1‑139 (MQ=255)                                                                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCg              >  2:174542/1‑139 (MQ=255)                                                                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCg              >  1:739819/1‑139 (MQ=255)                                                                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCg              >  1:629743/1‑139 (MQ=255)                                                                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCg              >  1:463881/1‑139 (MQ=255)                                                                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCg              >  1:319750/1‑139 (MQ=255)                                                                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCg              >  1:116251/1‑139 (MQ=255)                                                                                                                     ttCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCg              >  1:827397/1‑139 (MQ=255)                                                                                                                      tCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCgg             >  1:90833/1‑139 (MQ=255)                                                                                                                      tCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCgg             >  1:449151/1‑139 (MQ=255)                                                                                                                               gtgtTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGAc    <  1:794087/139‑1 (MQ=255)                                                                                                                                tgtTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACg   <  2:793690/139‑1 (MQ=255)                                                                                                                                  ttATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGc  <  2:90833/138‑1 (MQ=255)                                                                                                                                        |                                                                                                                                    TGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGC  >  NZ_CP009273/3882004‑3882270

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: ‑