Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,073,163	A→T	D26V (GAT→GTT)	BW25113_RS20185 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,073,163	0	A	T	97.4%	146.0 / ‑3.5	39	D26V (GAT→GTT)	BW25113_RS20185	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (21/17);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

GATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCA  >  NZ_CP009273/4073080‑4073248                                                                                    |                                                                                      gaTATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGa                                                                                 >  2:92517/1‑90 (MQ=255)      ttATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                            >  1:703214/1‑90 (MQ=255)            aCTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTa                                                                      <  1:632716/90‑1 (MQ=255)               ctctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATc                                                                   >  1:782238/1‑90 (MQ=255)               ctctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATc                                                                   >  2:501710/1‑90 (MQ=255)                   ggCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTcc                                                               <  1:501710/90‑1 (MQ=255)                       ggTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGt                                                           >  2:34817/1‑90 (MQ=255)                         tataGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGc                                                         >  2:306874/1‑90 (MQ=255)                           taGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCgg                                                       <  2:31297/90‑1 (MQ=255)                             gATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAc                                                     <  1:240913/90‑1 (MQ=255)                             gATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAc                                                     <  2:703214/90‑1 (MQ=255)                              atatGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAcc                                                    <  1:254846/90‑1 (MQ=255)                                 tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                                 >  1:471897/1‑90 (MQ=255)                                 tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                                 >  1:587378/1‑90 (MQ=255)                                 tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                                 >  2:163840/1‑90 (MQ=255)                                 tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                                 >  2:182184/1‑90 (MQ=255)                                  gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  2:655781/1‑90 (MQ=255)                                  gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  2:64029/1‑90 (MQ=255)                                  gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  2:526142/1‑90 (MQ=255)                                  gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  1:740952/1‑90 (MQ=255)                                  gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  1:156897/1‑90 (MQ=255)                                  gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  1:251878/1‑90 (MQ=255)                                   ggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTa                                               <  2:339362/90‑1 (MQ=255)                                    gCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTa                                               <  1:630371/89‑1 (MQ=255)                                    gCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTa                                               >  2:630371/1‑89 (MQ=255)                                       gAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAggg                                           >  2:282320/1‑90 (MQ=255)                                       gAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAggg                                           >  1:333355/1‑90 (MQ=255)                                        aaTTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGa                                          <  1:282320/90‑1 (MQ=255)                                         aTTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGaa                                         <  2:70215/90‑1 (MQ=255)                                           tttACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGc                                       <  2:42306/90‑1 (MQ=255)                                              aCTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGt                                    <  2:156897/90‑1 (MQ=255)                                              aCTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGt                                    <  1:220538/90‑1 (MQ=255)                                               cTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTa                                   <  2:414205/90‑1 (MQ=255)                                                        tCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATAc                          <  1:320064/90‑1 (MQ=255)                                                            aTGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCTTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGa                      >  1:194662/1‑90 (MQ=255)                                                             tGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCgg                                                       <  2:563858/56‑1 (MQ=255)                                                             tGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCgg                                                       >  1:563858/1‑56 (MQ=255)                                                                            aCTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAAc      <  2:332702/90‑1 (MQ=255)                                                                            aCTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAAc      >  1:551080/1‑90 (MQ=255)                                                                                gatgatCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCa  <  1:403596/90‑1 (MQ=255)                                                                                    |                                                                                      GATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCA  >  NZ_CP009273/4073080‑4073248

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