Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,406,008	A→G	D22D (GAT→GAC)	yfbU ←	YfbU family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,406,008	0	A	G	100.0%	196.1 / NA	51	D22D (GAT→GAC)	yfbU	YfbU family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/28);  total (23/28)

GCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACT  >  NZ_CP009273/2405927‑2406092                                                                                  |                                                                                     gcgATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATa                                                                             >  2:811468/1‑91 (MQ=255)      ccAATTCTCGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAAGGTTCAGCATTTGCCGGGTCGAGCATAGtcat                                                                        <  1:831690/91‑1 (MQ=255)      ccAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGtcat                                                                        <  1:112586/91‑1 (MQ=255)      ccAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGtcat                                                                        <  1:770284/91‑1 (MQ=255)      ccAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGtcat                                                                        <  1:406564/91‑1 (MQ=255)        aaTTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGtcatca                                                                      <  1:611988/91‑1 (MQ=255)          ttCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCatcatc                                                                    >  1:619872/1‑91 (MQ=255)              cGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGtcat                                                                        <  2:32740/83‑1 (MQ=255)              cGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGtcat                                                                        >  1:32740/1‑83 (MQ=255)              cGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGt                                                                >  2:114215/1‑91 (MQ=255)               gCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTa                                                               <  1:63228/91‑1 (MQ=255)                    tGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGAt                                                          <  2:839623/91‑1 (MQ=255)                     gTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGAtt                                                         <  1:361909/91‑1 (MQ=255)                          ccGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCAt                                                                              <  2:702289/65‑1 (MQ=255)                          ccGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCAt                                                                              >  1:702289/1‑65 (MQ=255)                          ccGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACa                                                    >  1:712748/1‑91 (MQ=255)                          ccGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACa                                                    >  2:623724/1‑91 (MQ=255)                            gTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGtcat                                                                        <  2:357124/69‑1 (MQ=255)                            gTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGtcat                                                                        >  1:357124/1‑69 (MQ=255)                            gTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACaaa                                                  >  1:309487/1‑91 (MQ=255)                            gTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACaaa                                                  >  2:28489/1‑91 (MQ=255)                            gTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACaaa                                                  >  2:65100/1‑91 (MQ=255)                            gTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACaaa                                                  >  1:343134/1‑91 (MQ=255)                            gTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACaaa                                                  >  2:564572/1‑91 (MQ=255)                            gTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGGCGAGCATAGCCATCATCTCGTACTGATTAGACaaa                                                  >  2:431187/1‑91 (MQ=255)                             tAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACaaaa                                                 >  2:3480/1‑91 (MQ=255)                                ccACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATaa                                              <  1:673362/91‑1 (MQ=255)                                  aCGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGa                                            <  1:741107/91‑1 (MQ=255)                                  aCGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGa                                            >  2:188603/1‑91 (MQ=255)                                       cAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTg                                       <  2:72407/91‑1 (MQ=255)                                               gTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGAc                                                     <  2:416108/69‑1 (MQ=255)                                               gTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGAc                                                     >  1:416108/1‑69 (MQ=255)                                               gTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTgg                               >  1:46636/1‑91 (MQ=255)                                                 ttGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTc                             >  2:63106/1‑91 (MQ=255)                                                 ttGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTc                             >  1:266584/1‑91 (MQ=255)                                                 ttGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTc                             >  2:332038/1‑91 (MQ=255)                                                    cAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAtt                          <  1:764613/91‑1 (MQ=255)                                                         gACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAt                     <  1:564572/91‑1 (MQ=255)                                                         gACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAt                     >  1:790851/1‑91 (MQ=255)                                                              tAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttgttg                <  1:162101/91‑1 (MQ=255)                                                              tAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttgttg                <  2:343134/91‑1 (MQ=255)                                                              tAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttgttg                <  1:3480/91‑1 (MQ=255)                                                              tAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttgttg                <  1:332038/91‑1 (MQ=255)                                                                 cGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACt             <  2:695723/91‑1 (MQ=255)                                                                 cGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACt             <  1:361028/91‑1 (MQ=255)                                                                     cAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTg         >  2:568846/1‑91 (MQ=255)                                                                         aTTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTTATTTCCATTGTTGACTCCTGTATc     <  1:606652/91‑1 (MQ=255)                                                                         aTTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATc     <  2:644232/91‑1 (MQ=255)                                                                         aTTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATc     <  2:46636/91‑1 (MQ=255)                                                                         aTTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATc     <  2:266584/91‑1 (MQ=255)                                                                         aTTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGCTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATc     <  2:799255/91‑1 (MQ=255)                                                                           ttGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCAc   <  2:758028/91‑1 (MQ=255)                                                                            tGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCAct  >  2:770448/1‑91 (MQ=255)                                                                                  |                                                                                     GCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACT  >  NZ_CP009273/2405927‑2406092

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