Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,786,467	T→C	I193V (ATT→GTT)	yibB ←	protein YibB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,786,467	0	T	C	94.0%	144.6 / 1.6	50	I193V (ATT→GTT)	yibB	protein YibB
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (39/8);  total (41/9)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGG  >  NZ_CP009273/3786383‑3786551                                                                                     |                                                                                     tttCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGa                                                                                 >  2:432690/1‑90 (MQ=255) tttCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGa                                                                                 >  1:55526/1‑90 (MQ=255)  ttCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGaa                                                                                >  1:41434/1‑90 (MQ=255)            aaTTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAAcc                                                                      >  2:788104/1‑90 (MQ=255)               ttATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCtgt                                                                   <  2:579664/90‑1 (MQ=255)               ttATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCtgt                                                                   <  1:270604/90‑1 (MQ=255)                    gAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGttgttg                                                              <  1:405369/90‑1 (MQ=255)                           ttCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGAtgtg                                                       >  2:607784/1‑90 (MQ=255)                                    tGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGc                                              <  2:788908/90‑1 (MQ=255)                                     gCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGcc                                             <  1:240993/90‑1 (MQ=255)                                       tgtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCt                                           >  1:333960/1‑90 (MQ=255)                                       tgtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCt                                           >  2:544192/1‑90 (MQ=255)                                        gtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCtt                                          >  2:806206/1‑90 (MQ=255)                                        gtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCtt                                          >  1:11147/1‑90 (MQ=255)                                        gtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCtt                                          >  1:574450/1‑90 (MQ=255)                                               aacaaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCt                                           >  2:339828/1‑82 (MQ=255)                                               aacaaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCt                                           <  1:339828/82‑1 (MQ=255)                                                  aaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGt                                >  2:160210/1‑90 (MQ=255)                                                    ttGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGa                              >  1:680222/1‑90 (MQ=255)                                                            ccAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGc                      <  1:515665/90‑1 (MQ=255)                                                             cAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCa                     >  1:764435/1‑90 (MQ=255)                                                                     tataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt             <  2:631910/90‑1 (MQ=255)                                                                     tataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt             >  2:11506/1‑90 (MQ=255)                                                                     tataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt             >  1:424150/1‑90 (MQ=255)                                                                     tataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt             >  2:336071/1‑90 (MQ=255)                                                                      ataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACttt            >  2:379748/1‑90 (MQ=255)                                                                       taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc           >  1:258729/1‑90 (MQ=255)                                                                       taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc           >  2:712877/1‑90 (MQ=255)                                                                         cATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcat         >  2:528455/1‑90 (MQ=255)                                                                         cATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcat         >  1:286913/1‑90 (MQ=255)                                                                         cATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcat         >  2:236456/1‑90 (MQ=255)                                                                           tGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcatca       >  1:432230/1‑90 (MQ=255)                                                                            gATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa      >  2:258784/1‑90 (MQ=255)                                                                            gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa      >  1:141722/1‑90 (MQ=255)                                                                            gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa      >  2:749092/1‑90 (MQ=255)                                                                            gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa      >  2:626866/1‑90 (MQ=255)                                                                            gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa      >  1:591695/1‑90 (MQ=255)                                                                            gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa      >  1:611320/1‑90 (MQ=255)                                                                            gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCAACTTACTTTCATCaa      >  2:67771/1‑90 (MQ=255)                                                                              ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGcc                                             <  2:67034/49‑1 (MQ=255)                                                                              ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGcc                                             >  1:67034/1‑49 (MQ=255)                                                                              ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt    >  1:566500/1‑90 (MQ=255)                                                                              ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt    >  1:312326/1‑90 (MQ=255)                                                                              ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt    >  2:174764/1‑90 (MQ=255)                                                                              ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt    >  1:658186/1‑90 (MQ=255)                                                                              ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt    >  1:672755/1‑90 (MQ=255)                                                                              ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt    >  1:141991/1‑90 (MQ=255)                                                                                gCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATgg  >  2:188424/1‑90 (MQ=255)                                                                                gCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATgg  >  2:579207/1‑90 (MQ=255)                                                                                gCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACCTTCATCAAATgg  >  1:756876/1‑90 (MQ=255)                                                                                     |                                                                                     TTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGG  >  NZ_CP009273/3786383‑3786551

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