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%	200.2 / NA	62	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 (33/29);  total (33/29)

TCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGC  >  NC_000913/2997568‑2997824                                                                                                                          |                                                                                                                                        tCTGCGTTTCTCCGGCCCTAAGACATCAATTATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttct                                                                                                                        >  2:1143166/1‑139 (MQ=255)   tGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttcttc                                                                                                                      >  2:1319678/1‑139 (MQ=255)                   tAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGa                                                                                                      <  2:312161/139‑1 (MQ=255)                          tcaatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGAttt                                                                                               >  1:322390/1‑139 (MQ=255)                             atcatcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTa                                                                                            <  1:1319678/139‑1 (MQ=255)                                        tcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                          >  2:569052/1‑130 (MQ=255)                                        tcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                          <  1:569052/130‑1 (MQ=255)                                               aCTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                              <  1:107269/119‑1 (MQ=255)                                               aCTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                              >  2:107269/1‑119 (MQ=255)                                                    tCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTg                                                                                                                    >  1:78695/1‑92 (MQ=255)                                                    tCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTg                                                                                                                    <  2:78695/92‑1 (MQ=255)                                                               gTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCtt                                                          <  2:322390/139‑1 (MQ=255)                                                                       aCTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATc                                                               >  2:559990/1‑126 (MQ=255)                                                                       aCTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATc                                                               <  1:559990/126‑1 (MQ=255)                                                                         tATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                          >  2:147762/1‑97 (MQ=255)                                                                         tATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                          <  1:147762/97‑1 (MQ=255)                                                                          aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                     <  1:233423/101‑1 (MQ=255)                                                                          aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                     >  2:233423/1‑101 (MQ=255)                                                                           tCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGtt                                              >  2:1223623/1‑139 (MQ=255)                                                                                tATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGAttt                                                                                               <  1:308365/85‑1 (MQ=255)                                                                                tATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGAttt                                                                                               >  2:308365/1‑85 (MQ=255)                                                                                     aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttcttc                                                                                                                      >  1:1335579/1‑57 (MQ=255)                                                                                     aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttcttc                                                                                                                      <  2:1335579/57‑1 (MQ=255)                                                                                     aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTc                                    >  2:452161/1‑139 (MQ=255)                                                                                      aGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTcc                                   <  1:452161/139‑1 (MQ=255)                                                                                       gTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataa                                                                                      >  2:407565/1‑87 (MQ=255)                                                                                       gTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataa                                                                                      <  1:407565/87‑1 (MQ=255)                                                                                         gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                          >  1:840488/1‑81 (MQ=255)                                                                                         gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                          <  2:840488/81‑1 (MQ=255)                                                                                            atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTa                                                                                            >  2:545537/1‑76 (MQ=255)                                                                                            atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTa                                                                                            <  1:545537/76‑1 (MQ=255)                                                                                            atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTCTCTTTTTCCAGATaa                             >  2:328236/1‑139 (MQ=255)                                                                                                gTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttcttc                                                                                                                      >  2:387237/1‑46 (MQ=255)                                                                                                gTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttcttc                                                                                                                      <  1:387237/46‑1 (MQ=255)                                                                                                gTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTg                         >  1:1284264/1‑139 (MQ=255)                                                                                                   gTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTgggg                      >  1:778476/1‑139 (MQ=255)                                                                                                    tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCtt                                                                                                                            >  1:288086/1‑36 (MQ=255)                                                                                                    tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCtt                                                                                                                            <  2:288086/36‑1 (MQ=255)                                                                                                    tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttcttc                                                                                                                      <  2:541467/42‑1 (MQ=255)                                                                                                    tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttcttc                                                                                                                      >  1:541467/1‑42 (MQ=255)                                                                                                    tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                              <  1:1103423/66‑1 (MQ=255)                                                                                                    tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                              >  2:1103423/1‑66 (MQ=255)                                                                                                    tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGa                     >  2:1532633/1‑139 (MQ=255)                                                                                                     cTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTg                                                                                                                    >  2:626144/1‑43 (MQ=255)                                                                                                     cTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTg                                                                                                                    <  1:626144/43‑1 (MQ=255)                                                                                                        gagaTTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                          >  1:470480/1‑66 (MQ=255)                                                                                                        gagaTTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                          <  2:470480/66‑1 (MQ=255)                                                                                                        gagaTTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATtct                 <  2:778476/139‑1 (MQ=255)                                                                                                        gagaTTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATtct                 <  2:1284264/139‑1 (MQ=255)                                                                                                             tCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                          <  2:1242362/61‑1 (MQ=255)                                                                                                             tCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                          >  1:1242362/1‑61 (MQ=255)                                                                                                                  gggCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAAt       >  2:702241/1‑139 (MQ=255)                                                                                                                    gCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                          >  1:183705/1‑54 (MQ=255)                                                                                                                    gCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                          <  2:183705/54‑1 (MQ=255)                                                                                                                    gCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAAt       >  2:751098/1‑137 (MQ=255)                                                                                                                    gCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAAt       <  1:751098/137‑1 (MQ=255)                                                                                                                    gCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCa     <  1:1532633/139‑1 (MQ=255)                                                                                                                      cAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCt                                                           >  2:1311979/1‑83 (MQ=255)                                                                                                                      cAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCt                                                           <  1:1311979/83‑1 (MQ=255)                                                                                                                      cAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAg   >  1:1018649/1‑139 (MQ=255)                                                                                                                      cAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAg   <  2:1018649/139‑1 (MQ=255)                                                                                                                       aGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGc  >  2:374125/1‑139 (MQ=255)                                                                                                                          |                                                                                                                                        TCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGC  >  NC_000913/2997568‑2997824

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