Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	746,726	G→A	*264* (TAG→TAA)  A348A (GCC→GCT)	nei → abrB ←	endonuclease VIII and 5‑formyluracil/5‑hydroxymethyluracil DNA glycosylase regulator of aidB expression, inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	746,726	0	G	A	100.0%	125.4 / NA	39	*264* (TAG→TAA)  A348A (GCC→GCT)	nei abrB	endonuclease VIII and 5‑formyluracil/5‑hydroxymethyluracil DNA glycosylase regulator of aidB expression, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/20);  total (19/20)

CATCATGGGGCGCTGTTTCGCTTTAAGGTTTTTCATCGAGATGGCGAACCGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCT  >  NC_000913/746595‑746846                                                                                                                                    |                                                                                                                         catcatGGGGCGCTGTTTCGCTTTAAGGTTTTTCATCGAGATGGCGAACCGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGAcc                                                                                                                   >  2:332535/1‑139 (MQ=255)     atGGGGCGCTGTTTCGCTTTAAGGTTTTTCATCGAGATGGCGAACCGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCtt                                                                                                               >  2:908180/1‑139 (MQ=255)      tGGGGCGCTGTTTCGCTTTAAGGTTTTTCATCGAGATGGCGAACCGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTc                                                                                                              >  1:435911/1‑139 (MQ=255)            gcTGTTTCGCTTTAAGGTTTTTCATCGAGATGGCGAACCGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCa                                                                                                        <  1:686906/139‑1 (MQ=255)                                                  cGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCt                                                                                                                >  2:291229/1‑93 (MQ=255)                                                  cGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCt                                                                                                                <  1:291229/93‑1 (MQ=255)                                                   gTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAgg                                                                                                    <  2:59985/104‑1 (MQ=255)                                                   gTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAgg                                                                                                    >  1:59985/1‑104 (MQ=255)                                                         aCGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATa                                                           <  2:435911/139‑1 (MQ=255)                                                              gtgGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCtt                                                                                                               <  2:973741/82‑1 (MQ=255)                                                              gtgGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCtt                                                                                                               >  1:973741/1‑82 (MQ=255)                                                              gtgGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGaa                                                      <  1:908180/139‑1 (MQ=255)                                                              gtgGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGaa                                                      <  1:43487/139‑1 (MQ=255)                                                               tgGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGaat                                                     <  1:120886/139‑2 (MQ=255)                                                                ggCAGCATCATTTAGTGATCCTCGCTGTCATGTCGCCCGATTCAGTGGTGCCCTGGCTGCCAGCGATAAGCAGACCGCTCCGGCGCATTGGTGGAAATACTCCCAGCAATGACAGGCCCTGGCAGCAAAATACTGAGCa                                                    <  1:882983/139‑1 (MQ=255)                                                                     catcatTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCa                                                                                                        <  1:207007/82‑1 (MQ=255)                                                                     catcatTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCa                                                                                                        >  2:207007/1‑82 (MQ=255)                                                                      atcatTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGaa                                                      <  1:1000595/131‑1 (MQ=255)                                                                      atcatTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGaa                                                      >  2:1000595/1‑131 (MQ=255)                                                                         atTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACcgcg                                                                                     <  1:137738/97‑1 (MQ=255)                                                                         atTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACcgcg                                                                                     >  2:137738/1‑97 (MQ=255)                                                                                  aCCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCg                                  >  2:500685/1‑139 (MQ=255)                                                                                    cACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTg                                       >  1:824036/1‑132 (MQ=255)                                                                                    cACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTg                                       <  2:824036/132‑1 (MQ=255)                                                                                         tGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAAt                           <  1:500685/139‑1 (MQ=255)                                                                                             aTCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAAt                           <  2:984617/135‑1 (MQ=255)                                                                                             aTCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAAt                           >  1:984617/1‑135 (MQ=255)                                                                                                     cGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCa                                                                                                        <  2:238270/50‑1 (MQ=255)                                                                                                     cGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCa                                                                                                        <  1:815921/50‑1 (MQ=255)                                                                                                     cGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCa                                                                                                        >  2:815921/1‑50 (MQ=255)                                                                                                     cGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCa                                                                                                        >  1:238270/1‑50 (MQ=255)                                                                                                     cGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCa                                     <  2:599865/117‑1 (MQ=255)                                                                                                     cGTTTTACTGGTGCCCTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCa                                     >  1:599865/1‑117 (MQ=255)                                                                                                                  cccTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCgcgc                                                                                    <  2:346475/57‑1 (MQ=255)                                                                                                                  cccTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCgcgc                                                                                    >  1:346475/1‑57 (MQ=255)                                                                                                                  cccTGGCTGCCAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCt  >  2:128984/1‑139 (MQ=255)                                                                                                                            cAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGccc                                                                        >  2:612055/1‑59 (MQ=255)                                                                                                                            cAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGccc                                                                        <  1:612055/59‑1 (MQ=255)                                                                                                                            cAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCACTGTCAGCAAAATACTGaa                                                      >  1:676719/1‑77 (MQ=255)                                                                                                                            cAGCATTAAGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCACTGTCAGCAAAATACTGaa                                                      <  2:676719/77‑1 (MQ=255)                                                                                                                                    |                                                                                                                         CATCATGGGGCGCTGTTTCGCTTTAAGGTTTTTCATCGAGATGGCGAACCGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCT  >  NC_000913/746595‑746846

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