Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,786,991	T→C	D18G (GAT→GGT)	yibB ←	protein YibB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,786,991	0	T	C	100.0%	179.6 / NA	50	D18G (GAT→GGT)	yibB	protein YibB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/32);  total (18/32)

TCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGA  >  NZ_CP009273/3786904‑3787073                                                                                        |                                                                                   tCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAAccc                                                                                  <  1:613558/90‑1 (MQ=255) tCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAAccc                                                                                  <  1:211333/90‑1 (MQ=255)   aaGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACcccc                                                                                >  1:392314/1‑90 (MQ=255)      tGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctct                                                                             >  1:756588/1‑90 (MQ=255)      tGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctct                                                                             >  2:381403/1‑90 (MQ=255)       gCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTa                                                                            >  2:29659/1‑90 (MQ=255)          gCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCa                                                                         <  1:684845/90‑1 (MQ=255)           ccAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCaa                                                                        <  1:167672/90‑1 (MQ=255)                cGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCa                                                                   <  2:392314/90‑1 (MQ=255)                cGCACAAAGTACCTAAAATAAACATCACCCCAACGATCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCa                                                                   <  2:379807/90‑1 (MQ=255)                     aaaGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAata                                                              <  2:633559/90‑1 (MQ=255)                       aGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATg                                                            <  1:231774/90‑1 (MQ=255)                        gTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGc                                                           <  1:652963/90‑1 (MQ=255)                         tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  2:477798/90‑1 (MQ=255)                         tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  2:367728/90‑1 (MQ=255)                         tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  1:582812/90‑1 (MQ=255)                         tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  1:600132/90‑1 (MQ=255)                           aCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGt                                                        >  1:628353/1‑90 (MQ=255)                            cTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       <  2:678886/90‑1 (MQ=255)                              aaaaTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGAt                                                     <  1:29659/90‑1 (MQ=255)                                aaTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATaa                                                   >  1:250940/1‑90 (MQ=255)                                 aTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAAt                                                  <  1:686994/90‑1 (MQ=255)                                 aTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAAt                                                  <  2:893807/90‑1 (MQ=255)                                 aTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAAt                                                  <  2:359110/90‑1 (MQ=255)                                 aTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAAt                                                  <  1:179730/90‑1 (MQ=255)                                        tCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTgtggtg                                           >  2:255141/1‑90 (MQ=255)                                            cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  1:682987/90‑1 (MQ=255)                                            cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  1:798534/90‑1 (MQ=255)                                             tGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatga                                      >  2:345732/1‑90 (MQ=255)                                              gAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatgat                                     >  1:743570/1‑90 (MQ=255)                                                 cGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTc                                  <  1:381403/90‑1 (MQ=255)                                                   aGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCAt                                >  1:664266/1‑90 (MQ=255)                                                   aGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCAt                                <  1:658286/90‑1 (MQ=255)                                                     cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  2:98857/90‑1 (MQ=255)                                                     cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  1:345732/90‑1 (MQ=255)                                                     cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  2:302998/90‑1 (MQ=255)                                                      tAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTAc                             <  2:664266/90‑1 (MQ=255)                                                           tttCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACtaata                        <  1:225781/90‑1 (MQ=255)                                                             tCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAAt                                                                       >  1:451739/1‑41 (MQ=255)                                                             tCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAAt                                                                       <  2:451739/41‑1 (MQ=255)                                                             tCACGGAACCCTTTATTAGCAGTCCAACCCCCTCGACCAATGt                                                                     >  1:138165/1‑43 (MQ=255)                                                             tCACCGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGt                                                                     <  2:138165/43‑1 (MQ=255)                                                              cACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGa                     <  1:565746/90‑1 (MQ=255)                                                                               gCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATaaaaaaa    >  1:598952/1‑90 (MQ=255)                                                                                cAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       <  2:836013/38‑1 (MQ=39)                                                                                cAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       >  1:836013/1‑38 (MQ=255)                                                                                 aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       <  1:429141/37‑1 (MQ=39)                                                                                 aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       >  2:429141/1‑37 (MQ=39)                                                                                 aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGa  >  1:442373/1‑90 (MQ=255)                                                                                 aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGa  >  2:118190/1‑90 (MQ=255)                                                                                  gTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  2:540518/52‑1 (MQ=255)                                                                                  gTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       >  1:540518/1‑52 (MQ=255)                                                                                        |                                                                                   TCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGA  >  NZ_CP009273/3786904‑3787073

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