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	92.9%	76.9 / ‑0.5	28	I193V (ATT→GTT)	yibB	protein YibB
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (18/8);  total (19/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

AATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGA  >  NZ_CP009273/3786394‑3786553                                                                          |                                                                                       aaTTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAAcc                                                                        >  1:13383/1‑90 (MQ=255) aaTTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAAcc                                                                        >  1:44197/1‑90 (MQ=255)         gAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGttgttg                                                                <  2:13383/90‑1 (MQ=255)           aTTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTgag                                                              >  1:339930/1‑90 (MQ=255)                     tttATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAAcc                                                                        <  2:674119/70‑1 (MQ=255)                     tttATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAAcc                                                                        >  1:674119/1‑70 (MQ=255)                      ttATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTa                                                   >  1:663692/1‑90 (MQ=255)                          gCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGcc                                               <  2:94882/90‑1 (MQ=255)                             gtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCtt                                            >  2:31665/1‑90 (MQ=255)                                  ccAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCttttttt                                       >  2:259735/1‑90 (MQ=255)                                    aacaaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTtat                                     <  2:83162/90‑1 (MQ=255)                                    aacaaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTtat                                     <  2:107342/90‑1 (MQ=255)                                      caaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAg                                   <  1:134756/90‑1 (MQ=255)                                      caaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAg                                   <  2:184068/90‑1 (MQ=255)                                               ccccAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGAt                          <  1:194298/90‑1 (MQ=255)                                                          tataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt               >  2:12681/1‑90 (MQ=255)                                                          tataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt               >  2:182482/1‑90 (MQ=255)                                                            taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc             >  1:355451/1‑90 (MQ=255)                                                            taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc             >  2:312119/1‑90 (MQ=255)                                                            taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc             >  2:677533/1‑90 (MQ=255)                                                               aTGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcatc          >  2:154904/1‑90 (MQ=255)                                                                tGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcatca         >  1:522070/1‑90 (MQ=255)                                                                 gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGGTGCATCTTACTTTCATCaa        >  2:600797/1‑90 (MQ=255)                                                                 gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa        >  1:239236/1‑90 (MQ=255)                                                                   ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt      >  1:268917/1‑90 (MQ=255)                                                                   ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt      >  1:13929/1‑90 (MQ=255)                                                                     gCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATgg    >  2:130285/1‑90 (MQ=255)                                                                       cAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGa  <  1:391469/90‑1 (MQ=255)                                                                          |                                                                                       AATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGA  >  NZ_CP009273/3786394‑3786553

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