Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,809,775	A→T	D551V (GAT→GTT)	katE →	catalase HPII

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,809,775	0	A	T	98.2%	210.9 / ‑3.4	57	D551V (GAT→GTT)	katE	catalase HPII
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (28/28);  total (29/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 = 7.18e-01

ATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGG  >  NZ_CP009273/1809643‑1809910                                                                                                                                     |                                                                                                                                        aTCTGCGTCTGTTCTGGCTAAGTCAGACGCGATTTGAACAGCGCCATATTGACGATGGATTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACTAGCTGGCGCATATTGTTCTCAc                                                                                                                                   >  2:1282628/1‑139 (MQ=255) aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTtct                                                                                                                                      >  1:135067/1‑136 (MQ=255) aTCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTtct                                                                                                                                      <  2:135067/136‑1 (MQ=255)  tCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCAct                                                                                                                                  >  1:1218423/1‑139 (MQ=255)  tCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCAct                                                                                                                                  >  1:345878/1‑139 (MQ=255)      cgTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTgg                                                                                                                              >  2:689480/1‑139 (MQ=255)       gTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGc                                                                                                                             >  1:830880/1‑139 (MQ=255)            ttCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTtct                                                                                                                                      >  2:856993/1‑125 (MQ=255)            ttCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTtct                                                                                                                                      <  1:856993/125‑1 (MQ=255)                 gCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTg                                                                                                                   >  2:739013/1‑139 (MQ=255)                      ggcAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGa                                                                                                               >  1:291256/3‑138 (MQ=255)                      gTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCAct                                                                                                                                  <  2:929212/119‑1 (MQ=255)                      gTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCAct                                                                                                                                  >  1:929212/1‑119 (MQ=255)                      gTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGa                                                                                                               >  1:125217/1‑138 (MQ=255)                      gTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGa                                                                                                               >  1:835173/1‑138 (MQ=255)                            cGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAAtct                                                                                                        <  1:679198/139‑1 (MQ=255)                            cGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAAtct                                                                                                        <  1:291315/139‑1 (MQ=255)                               cATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCgg                                                                                                     >  1:396048/1‑139 (MQ=255)                                 tttGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTa                                                                                                   >  2:711545/1‑139 (MQ=255)                                      gcagcGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGaa                                                                                              <  2:479955/139‑1 (MQ=255)                                        agcGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAact                                                                                            <  2:1125590/139‑1 (MQ=255)                                          cgcCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAactga                                                                                          >  2:590019/1‑139 (MQ=255)                                          cgcCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAactga                                                                                          >  1:260813/1‑139 (MQ=255)                                                   gTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTgacgac                                                                                  <  1:475034/138‑1 (MQ=255)                                                   gTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTgacgac                                                                                  <  1:1270480/138‑1 (MQ=255)                                                    tCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTgacgac                                                                                  >  2:1270480/1‑137 (MQ=255)                                                    tCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTgacgac                                                                                  >  1:1086151/1‑137 (MQ=255)                                                      gATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAg                                                                               <  2:1086151/138‑1 (MQ=255)                                                              cAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAatat                                                                       <  2:260813/138‑1 (MQ=255)                                                              cAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAatat                                                                       >  1:177122/1‑138 (MQ=255)                                                                 ttttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCAcc                                                                   <  1:711545/139‑1 (MQ=255)                                                                   ttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCAccc                                                                  <  2:835173/138‑1 (MQ=255)                                                                   ttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCAccc                                                                  <  2:345878/138‑1 (MQ=255)                                                                   ttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCAccc                                                                  <  2:1218423/138‑1 (MQ=255)                                                                   ttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTAACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCAccc                                                                  <  2:729122/138‑1 (MQ=255)                                                                        ttAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTTACGACCAGCTGAATATCACCCCACCTc                                                            <  2:291256/139‑1 (MQ=255)                                                                                aGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGt                                                     >  1:378949/1‑138 (MQ=255)                                                                                    gttcgtCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACg                                                <  1:1235012/139‑1 (MQ=255)                                                                                     ttcgtCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACgg                                               >  1:987316/1‑139 (MQ=255)                                                                                       cgtCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGt                                              <  2:396048/138‑1 (MQ=255)                                                                                            gTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGaaa                                        <  1:682175/139‑1 (MQ=255)                                                                                            gTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGaaa                                        <  1:764718/139‑1 (MQ=255)                                                                                             tataTTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGaaaa                                       <  1:887675/139‑1 (MQ=255)                                                                                                  tCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGACAAAGGac                                  >  2:749522/1‑138 (MQ=255)                                                                                                   cGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGat                                  <  1:590019/138‑1 (MQ=255)                                                                                                    gTGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGTACGTCAACGGTCTGAAAAAGGatcc                                <  2:1004667/139‑1 (MQ=255)                                                                                                     tGAGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGatcca                               >  2:46315/1‑139 (MQ=255)                                                                                                       aGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCAct                                                                                                                                  >  2:218593/1‑38 (MQ=255)                                                                                                       aGCGCGTTGTTGACCAGCTGGCGCATATTGTTCTCAct                                                                                                                                  <  1:218593/38‑1 (MQ=255)                                                                                                         cgcgTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCt                           >  2:710055/1‑139 (MQ=255)                                                                                                         cgcgTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCt                           >  1:23273/1‑139 (MQ=255)                                                                                                           cgTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGatc                                 <  1:24758/131‑1 (MQ=255)                                                                                                           cgTTGTTGACCAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGatc                                 >  2:24758/1‑131 (MQ=255)                                                                                                                     cAGCTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACg               <  2:1169064/139‑1 (MQ=255)                                                                                                                        cTGGCGCATATTGTTCTCACTCTGGCCCAGGCGGCGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATTCATCCTTAAGTTTGTACGCCa            >  2:1134553/1‑139 (MQ=255)                                                                                                                                  ttGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACgg  >  2:330082/1‑139 (MQ=255)                                                                                                                                  ttGTTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACgg  <  1:991276/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        ATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGG  >  NZ_CP009273/1809643‑1809910

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑