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	93.5%	175.3 / 4.4	62	I193V (ATT→GTT)	yibB	protein YibB
Reads supporting (aligned to +/- strand):  ref base T (4/0);  new base C (41/17);  total (45/17)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

TTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGA  >  NZ_CP009273/3786384‑3786553                                                                                    |                                                                                       ttCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGaa                                                                                  <  2:1371395/90‑1 (MQ=255) ttCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGaa                                                                                  <  1:1451275/90‑1 (MQ=255)           aaTTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAAcc                                                                        >  2:1511773/1‑90 (MQ=255)           aaTTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAAcc                                                                        >  2:389159/1‑90 (MQ=255)           aaTTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAAcc                                                                        <  1:1434083/90‑1 (MQ=255)              ttATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCtgt                                                                     <  1:1362512/90‑1 (MQ=255)              ttATAGAATACTTCCCATTTACGCTGTTAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAACCCtgt                                                                     <  2:937596/90‑1 (MQ=255)                  aGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTgttgtt                                                                 <  1:462160/90‑1 (MQ=255)                     aTTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTgag                                                              >  2:1208305/1‑90 (MQ=255)                     aTTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTgag                                                              >  1:1060988/1‑90 (MQ=255)                                      tgtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCt                                             <  1:127724/90‑1 (MQ=255)                                      tgtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCt                                             >  1:726454/1‑90 (MQ=255)                                       gtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGcc                                               <  2:867408/87‑1 (MQ=255)                                       gtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGcc                                               >  1:867408/1‑87 (MQ=255)                                       gtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCtt                                            >  1:77822/1‑90 (MQ=255)                                       gtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCtt                                            >  1:472820/1‑90 (MQ=255)                                       gtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCtt                                            >  1:1142328/1‑90 (MQ=255)                                       gtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCtt                                            >  2:909313/1‑90 (MQ=255)                                         gAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGcc                                               >  1:1152425/1‑85 (MQ=255)                                         gAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGcc                                               <  2:1152425/85‑1 (MQ=255)                                         gAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCtttt                                          >  2:1507980/1‑90 (MQ=255)                                              aacagTTGCGCCCCCCATCATATATGCATGCTTGCCCACCATCAGCTCCTCCCCCTGTTGTTGAGATGTGTCAGTCAGCCCtttttttct                                     <  1:1411682/90‑3 (MQ=255)                                              aacaaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTtat                                     <  1:205534/90‑1 (MQ=255)                                               acaaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTtata                                    >  2:138661/1‑90 (MQ=255)                                                   ttGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGa                                <  1:282620/90‑1 (MQ=255)                                                   ttGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGa                                <  1:554306/90‑1 (MQ=255)                                                         ccccAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTtat                                     >  2:1339094/1‑79 (MQ=255)                                                         ccccAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTtat                                     <  1:1339094/79‑1 (MQ=255)                                                                  tatataCATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTAc                 >  2:1525822/1‑90 (MQ=255)                                                                    tataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt               >  2:404322/1‑90 (MQ=255)                                                                      taCATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc             >  1:447521/1‑90 (MQ=255)                                                                      taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc             >  2:503913/1‑90 (MQ=255)                                                                      taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc             >  2:582651/1‑90 (MQ=255)                                                                      taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc             >  1:176302/1‑90 (MQ=255)                                                                      taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc             >  2:920889/1‑90 (MQ=255)                                                                      taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc             >  2:996072/1‑90 (MQ=255)                                                                        cATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcat           >  2:232880/1‑90 (MQ=255)                                                                        cATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcat           >  2:31739/1‑90 (MQ=255)                                                                        cATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcat           >  2:1315738/1‑90 (MQ=255)                                                                          tGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcatca         >  2:835544/1‑90 (MQ=255)                                                                          tGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcatca         >  2:518423/1‑90 (MQ=255)                                                                           gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa        >  2:207812/1‑90 (MQ=255)                                                                           gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa        >  1:1280496/1‑90 (MQ=255)                                                                           gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa        >  1:1552671/1‑90 (MQ=255)                                                                           gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa        >  2:438223/1‑90 (MQ=255)                                                                           gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa        >  2:686001/1‑90 (MQ=255)                                                                             ttGCCAACCATGAAATCAAAAACCTGTTGTTgaga                                                             <  2:1200593/35‑1 (MQ=38)                                                                             ttGCCAACCATGAAATCAAAAACCTGTTGTTgaga                                                             >  1:1200593/1‑35 (MQ=38)                                                                             ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATg                         >  2:1154050/1‑71 (MQ=255)                                                                             ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATg                         <  1:1154050/71‑1 (MQ=255)                                                                             ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt      >  1:110/1‑90 (MQ=255)                                                                             ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt      >  2:77659/1‑90 (MQ=255)                                                                             ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt      >  1:268391/1‑90 (MQ=255)                                                                             ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt      >  1:798157/1‑90 (MQ=255)                                                                             ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt      >  2:395294/1‑90 (MQ=255)                                                                             ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt      >  2:390184/1‑90 (MQ=255)                                                                               gCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATgg    >  2:755660/1‑90 (MQ=255)                                                                               gCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATgg    >  2:410992/1‑90 (MQ=255)                                                                               gCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATgg    >  2:400257/1‑90 (MQ=255)                                                                               gCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATgg    >  2:94534/1‑90 (MQ=255)                                                                               gCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATgg    >  2:1454716/1‑90 (MQ=255)                                                                                 cAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGa  <  1:611642/90‑1 (MQ=255)                                                                                    |                                                                                       TTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGA  >  NZ_CP009273/3786384‑3786553

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