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%	154.8 / NA	49	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/22);  total (27/22)

TGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCA  >  NC_000913/2997558‑2997821                                                                                                                                    |                                                                                                                                     tgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACtt                                                                                                                               >  1:426441/1‑139 (MQ=14)       gCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCtt                                                                                                                         <  2:882095/139‑1 (MQ=37)            cTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttctt                                                                                                                    <  2:426441/139‑1 (MQ=255)                     tcCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGttt                                                                                                           >  1:907408/1‑139 (MQ=255)                      ccGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTa                                                                                                          <  2:709676/139‑1 (MQ=255)                       cGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTAt                                                                                                         >  2:51479/1‑139 (MQ=255)                             tAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGa                                                                                                   <  1:326462/139‑1 (MQ=255)                              aaGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGaa                                                                                                  <  2:907408/139‑1 (MQ=255)                                                  tcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAAc                                                                              >  1:749001/1‑139 (MQ=255)                                                      aTGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTa                                                                          >  1:128173/1‑139 (MQ=255)                                                            aGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                           <  1:506442/116‑1 (MQ=255)                                                            aGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                           >  2:506442/1‑116 (MQ=255)                                                             gTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTa                                                                                         <  1:430338/117‑1 (MQ=255)                                                             gTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTa                                                                                         >  2:430338/1‑117 (MQ=255)                                                                        aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCt                                                        >  1:120028/1‑139 (MQ=255)                                                                        aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCt                                                        >  1:314247/1‑139 (MQ=255)                                                                         gTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCtt                                                       <  2:120028/139‑1 (MQ=255)                                                                                 aCTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGACAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCgg                                                          >  1:162824/1‑128 (MQ=255)                                                                                 aCTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCgg                                                          <  2:162824/128‑1 (MQ=255)                                                                                    aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGt                                            >  1:875100/1‑139 (MQ=255)                                                                                    aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGt                                            >  2:774771/1‑139 (MQ=255)                                                                                    aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGt                                            >  2:119878/1‑139 (MQ=255)                                                                                           aTTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                  <  2:756021/94‑1 (MQ=255)                                                                                           aTTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                  >  1:756021/1‑94 (MQ=255)                                                                                             ttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaat                                                                                 >  2:120258/1‑93 (MQ=255)                                                                                             ttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaat                                                                                 <  1:120258/93‑1 (MQ=255)                                                                                             ttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                  <  2:69316/92‑1 (MQ=255)                                                                                             ttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                  >  1:69316/1‑92 (MQ=255)                                                                                               aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTc                                 >  2:3533/1‑139 (MQ=255)                                                                                                   gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCtt                                                       <  2:321878/113‑1 (MQ=255)                                                                                                   gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCtt                                                       >  1:321878/1‑113 (MQ=255)                                                                                                   gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGa                             >  2:186390/1‑139 (MQ=255)                                                                                                   gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGa                             <  1:186390/139‑1 (MQ=255)                                                                                                      atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTa                                                                                                          >  2:750377/1‑59 (MQ=255)                                                                                                      atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTa                                                                                                          <  1:750377/59‑1 (MQ=255)                                                                                                      atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATaa                          <  2:749001/139‑1 (MQ=255)                                                                                                             gTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTgggg                   >  1:245892/1‑139 (MQ=255)                                                                                                              tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGAt                            <  2:381205/129‑1 (MQ=255)                                                                                                              tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGAt                            >  1:381205/1‑129 (MQ=255)                                                                                                              tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGa                  >  1:88981/1‑139 (MQ=255)                                                                                                                    gaTTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                           >  1:648538/1‑60 (MQ=255)                                                                                                                    gaTTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                           <  2:648538/60‑1 (MQ=255)                                                                                                                        cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGt                                                                       <  1:706578/76‑1 (MQ=255)                                                                                                                        cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGt                                                                       >  2:706578/1‑76 (MQ=255)                                                                                                                        cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATc                                                            >  1:327834/1‑87 (MQ=255)                                                                                                                        cAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATc                                                            <  2:327834/87‑1 (MQ=255)                                                                                                                              gCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCa  <  1:3533/139‑1 (MQ=255)                                                                                                                               ccAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGcc                                                                 <  2:755628/75‑1 (MQ=255)                                                                                                                               ccAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGcc                                                                 >  1:755628/1‑75 (MQ=255)                                                                                                                                    |                                                                                                                                     TGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCA  >  NC_000913/2997558‑2997821

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