Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,997,689	C→T	intergenic (‑89/+3)	insC1 ← / ← ygeQ	IS2 repressor TnpA/pseudogene, glycosyl hydrolase family 15, part of T3SS PAI ETT2 remnant

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,997,689	0	C	T	100.0%	171.1 / NA	53	intergenic (‑89/+3)	insC1/ygeQ	IS2 repressor TnpA/pseudogene, glycosyl hydrolase family 15, part of T3SS PAI ETT2 remnant
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/26);  total (27/26)

GTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATC  >  NC_000913/2997560‑2997820                                                                                                                                  |                                                                                                                                    gTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGc                                                                                                                            >  2:1429826/1‑139 (MQ=18)     gCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCtt                                                                                                                        <  1:1429826/139‑1 (MQ=37)           tGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttcttc                                                                                                                  <  1:821411/139‑1 (MQ=255)                          cTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTg                                                                                                   >  2:220562/1‑139 (MQ=255)                                  tcaatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGAttt                                                                                           >  2:946923/1‑139 (MQ=255)                                   caatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                          <  1:946923/139‑1 (MQ=255)                                                   aaTGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACAtt                                                                          <  1:220562/139‑1 (MQ=255)                                                    aTGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTa                                                                         >  1:285843/1‑139 (MQ=255)                                                          aGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAg                                                                                               >  1:799611/1‑111 (MQ=255)                                                          aGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAg                                                                                               <  2:799611/111‑1 (MQ=255)                                                                      aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCt                                                       >  2:789248/1‑139 (MQ=255)                                                                      aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCt                                                       >  2:1653378/1‑139 (MQ=255)                                                                            aaGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGctc                                                                                                                          <  1:822184/66‑1 (MQ=255)                                                                            aaGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGctc                                                                                                                          >  2:822184/1‑66 (MQ=255)                                                                             aGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGt                                                <  1:1653378/139‑1 (MQ=255)                                                                                  aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                      >  1:932908/1‑96 (MQ=255)                                                                                  aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                      <  2:932908/96‑1 (MQ=255)                                                                                  aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaata                                                                               <  1:328574/103‑1 (MQ=255)                                                                                  aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaata                                                                               >  2:328574/1‑103 (MQ=255)                                                                                  aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGat                                                                                    >  1:429053/1‑98 (MQ=255)                                                                                  aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGat                                                                                    <  2:429053/98‑1 (MQ=255)                                                                                  aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGt                                           >  2:115183/1‑139 (MQ=255)                                                                                   tCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                      >  2:1517879/1‑95 (MQ=255)                                                                                   tCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                      <  1:1517879/95‑1 (MQ=255)                                                                                         aTTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGa                                                                                                  >  2:1575319/1‑77 (MQ=255)                                                                                         aTTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGa                                                                                                  <  1:1575319/77‑1 (MQ=255)                                                                                           ttAAGTTATATTGGTTGTCTGGAGATTCATGGGGCCAGTTTAATCTTGCTCTTTTCTTCGGGTTTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCtttt                                  <  2:285843/139‑1 (MQ=39)                                                                                                 gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaata                                                                               >  2:1041006/1‑88 (MQ=255)                                                                                                 gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTCTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaatc                                                                               <  1:1041006/88‑2 (MQ=255)                                                                                                   tatTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCttttt                                 <  2:455373/132‑1 (MQ=255)                                                                                                   tatTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCttttt                                 >  1:455373/1‑132 (MQ=255)                                                                                                    atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGaaa                                                                                                >  1:1491334/1‑68 (MQ=255)                                                                                                    atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGaaa                                                                                                <  2:1491334/68‑1 (MQ=255)                                                                                                            tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                          >  1:242601/1‑66 (MQ=255)                                                                                                            tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                          <  2:242601/66‑1 (MQ=255)                                                                                                            tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGc                                                              <  2:269560/94‑1 (MQ=255)                                                                                                            tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGc                                                              >  1:269560/1‑94 (MQ=255)                                                                                                             cTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCAt                                                            >  2:1244949/1‑95 (MQ=255)                                                                                                             cTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCAt                                                            <  1:1244949/95‑1 (MQ=255)                                                                                                                gagaTTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCgg                                                         >  2:198783/1‑95 (MQ=255)                                                                                                                gagaTTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCgg                                                         <  1:198783/95‑1 (MQ=255)                                                                                                                  gaTTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAg                                                                                               >  2:1302238/1‑55 (MQ=255)                                                                                                                  gaTTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAg                                                                                               <  1:1302238/55‑1 (MQ=255)                                                                                                                      cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGt                                                                                         >  2:170085/1‑57 (MQ=255)                                                                                                                      cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGt                                                                                         <  1:170085/57‑1 (MQ=255)                                                                                                                      cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCt       >  2:1578826/1‑139 (MQ=255)                                                                                                                       aGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCa                              >  1:683739/1‑115 (MQ=255)                                                                                                                       aGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCa                              <  2:683739/115‑1 (MQ=255)                                                                                                                       aGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATtct             >  1:1213980/1‑132 (MQ=255)                                                                                                                       aGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGAAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATtct             <  2:1213980/132‑1 (MQ=255)                                                                                                                           ggCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATc  <  1:115183/139‑1 (MQ=255)                                                                                                                              cAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATtct             >  1:436067/1‑125 (MQ=255)                                                                                                                              cAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATtct             <  2:436067/125‑1 (MQ=255)                                                                                                                                  |                                                                                                                                    GTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATC  >  NC_000913/2997560‑2997820

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