Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	49,769	G→A	intergenic (+138/‑54)	kefC → / → folA	glutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	49,769	0	G	A	97.9%	192.0 / ‑4.3	48	intergenic (+138/‑54)	kefC/folA	glutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (25/22);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

ACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCC  >  NZ_CP009273/49641‑49900                                                                                                                                 |                                                                                                                                    aCGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACg                                                                                                                           <  1:701320/139‑1 (MQ=255) aCGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACg                                                                                                                           >  1:436900/1‑139 (MQ=255)                        aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAAtt                                                                                                    >  1:630993/1‑139 (MQ=255)                          tgtttgttTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAAtttt                                                                                                  >  1:77457/1‑139 (MQ=255)                                             gTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCa                                                                                <  1:520704/138‑1 (MQ=255)                                                     gcgAGTCCAGGGAGAGATCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGt                                                                       >  2:82077/1‑139 (MQ=255)                                                     gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGt                                                                       <  2:630993/139‑1 (MQ=255)                                                     gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGt                                                                       >  2:594954/1‑139 (MQ=255)                                                     gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGt                                                                       <  1:233449/139‑1 (MQ=255)                                                        aGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGa                                                                                       <  2:497591/120‑1 (MQ=255)                                                        aGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGa                                                                                       >  1:497591/1‑120 (MQ=255)                                                        aGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAg                                                                        <  2:321375/135‑1 (MQ=255)                                                        aGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAg                                                                        >  1:321375/1‑135 (MQ=255)                                                              gggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCATTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcg                                                              <  2:97516/139‑1 (MQ=255)                                                                  gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                          <  1:16650/139‑1 (MQ=255)                                                                   agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTAc                                                                                                                            <  1:567399/72‑1 (MQ=255)                                                                   agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTAc                                                                                                                            >  2:567399/1‑72 (MQ=255)                                                                   agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCATTCTGATTGCGGCGt                                                          >  2:121196/1‑138 (MQ=255)                                                                    gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                          >  1:352324/1‑137 (MQ=255)                                                                    gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                          >  1:427092/1‑137 (MQ=255)                                                                    gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                          >  2:770894/1‑137 (MQ=255)                                                                    gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                          >  2:527084/1‑137 (MQ=255)                                                                    gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                          >  1:246878/1‑137 (MQ=255)                                                                    gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                          >  2:360877/1‑137 (MQ=255)                                                                         tGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTa                                                                                             >  1:193864/1‑97 (MQ=255)                                                                         tGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTa                                                                                             <  2:193864/97‑1 (MQ=255)                                                                         tGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGAt                                                                           <  1:122722/115‑1 (MQ=255)                                                                         tGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGAt                                                                           >  2:122722/1‑115 (MQ=255)                                                                         tGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                          <  2:272859/132‑1 (MQ=255)                                                                         tGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                          >  1:272859/1‑132 (MQ=255)                                                                                gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                          <  2:16284/125‑1 (MQ=255)                                                                                gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                          >  1:16284/1‑125 (MQ=255)                                                                                  cagcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGt                                          >  1:324597/1‑139 (MQ=255)                                                                                     cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                       <  2:174414/139‑1 (MQ=255)                                                                                          tataAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCa                                  >  2:820460/1‑139 (MQ=255)                                                                                            taAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCt                                                                                                                         <  2:168852/50‑1 (MQ=255)                                                                                            taAAATTTTCCTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCt                                                                                                                         >  1:168852/1‑50 (MQ=255)                                                                                                      cTCAACATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGcc                       <  2:697042/138‑1 (MQ=255)                                                                                                         aaCATCATCCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCAt                     <  1:770894/137‑1 (MQ=255)                                                                                                                tcCTCGCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa            <  1:669307/139‑1 (MQ=255)                                                                                                                     gCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCa                      >  1:233645/1‑124 (MQ=255)                                                                                                                     gCACCAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCa                      <  2:233645/124‑1 (MQ=255)                                                                                                                        ccAGTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcctg    <  2:436900/139‑1 (MQ=255)                                                                                                                           gTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCtgcctgcc  >  1:26871/1‑138 (MQ=255)                                                                                                                           gTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGa                              <  1:244680/110‑1 (MQ=255)                                                                                                                           gTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGa                              >  2:244680/1‑110 (MQ=255)                                                                                                                           gTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCa                      <  2:530121/118‑1 (MQ=255)                                                                                                                           gTCGACAACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCa                      >  1:530121/1‑118 (MQ=255)                                                                                                                                 |                                                                                                                                    ACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCC  >  NZ_CP009273/49641‑49900

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