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%	121.0 / NA	35	D18G (GAT→GGT)	yibB	protein YibB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/24);  total (11/24)

GGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGA  >  NZ_CP009273/3786912‑3787073                                                                                |                                                                                   ggCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTAcc                                                                          <  1:752374/90‑1 (MQ=255)      gACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGt                                                                     <  1:742024/90‑1 (MQ=255)                gTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGc                                                           <  2:669438/90‑1 (MQ=255)                 tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  1:415287/90‑1 (MQ=255)                  aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCa                                                                         <  2:837282/74‑1 (MQ=255)                  aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCa                                                                         >  1:837282/1‑74 (MQ=255)                      aaaaTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGAt                                                     >  2:679831/1‑90 (MQ=255)                        aaTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATaa                                                   >  1:857509/1‑90 (MQ=255)                        aaTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATaa                                                   >  2:275998/1‑90 (MQ=255)                              cATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTgtgg                                             <  2:45994/90‑1 (MQ=255)                              cATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTgtgg                                             >  1:506424/1‑90 (MQ=255)                                  aaCTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGga                                         <  2:286019/90‑1 (MQ=255)                                    cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  2:857509/90‑1 (MQ=255)                                    cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  2:361877/90‑1 (MQ=255)                                    cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  1:36624/90‑1 (MQ=255)                                           aGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCAt                                <  2:199785/90‑1 (MQ=255)                                             cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  2:834505/90‑1 (MQ=255)                                             cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  2:522117/90‑1 (MQ=255)                                             cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  2:773633/90‑1 (MQ=255)                                             cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  1:609254/90‑1 (MQ=255)                                              tAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTAc                             <  2:419198/90‑1 (MQ=255)                                              tAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTAc                             <  2:506424/90‑1 (MQ=255)                                              tAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTAc                             <  1:679831/90‑1 (MQ=255)                                                   tttCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACtaata                        <  1:153870/90‑1 (MQ=255)                                                   tttCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACtaata                        <  2:518009/90‑1 (MQ=255)                                                     tCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAg                      <  1:311514/90‑1 (MQ=255)                                                         ggAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCa                  >  2:445366/1‑90 (MQ=255)                                                                tttATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAatat           <  2:483781/90‑1 (MQ=255)                                                                    ttaGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATaaaa       >  1:411849/1‑90 (MQ=255)                                                                       gCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  2:217998/36‑1 (MQ=39)                                                                       gCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          >  1:217998/1‑36 (MQ=39)                                                                         aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGa  >  1:416317/1‑90 (MQ=255)                                                                         aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGa  >  1:256795/1‑90 (MQ=255)                                                                         aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGCGATAATTGTGGTGGATGATTTCATTa                              <  2:219899/62‑1 (MQ=255)                                                                         aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGCGATAATTGTGGTGGATGATTTCATTa                              >  1:219899/1‑62 (MQ=255)                                                                                |                                                                                   GGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGA  >  NZ_CP009273/3786912‑3787073

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: ‑