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%	193.2 / NA	61	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 (32/29);  total (32/29)

GGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCA  >  NC_000913/2997559‑2997821                                                                                                                                   |                                                                                                                                     ggTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTg                                                                                                                              <  2:1216164/139‑1 (MQ=17)        cGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCtttt                                                                                                                       >  2:248249/1‑139 (MQ=255)           cTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGCTATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttctt                                                                                                                    >  1:177525/1‑139 (MQ=255)            tGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttcttc                                                                                                                   <  1:1071137/139‑1 (MQ=255)            tGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttcttc                                                                                                                   >  1:1496496/1‑139 (MQ=255)                                        catcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                       <  2:177525/139‑1 (MQ=255)                                                 tcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAAc                                                                              >  1:1668017/1‑139 (MQ=255)                                                       gACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAg                                                                        <  1:248249/139‑1 (MQ=255)                                                        aCTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGaaa                                                                                                 <  1:924222/113‑1 (MQ=255)                                                        aCTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGaaa                                                                                                 >  2:924222/1‑113 (MQ=255)                                                        aCTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGt                                                                       >  2:289616/1‑139 (MQ=255)                                                             tCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATAGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGc                                                                  >  1:1133938/1‑139 (MQ=255)                                                                     cTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCgg                                                          <  2:1496496/139‑1 (MQ=255)                                                                       aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCATTTTAAACTTGCTCTTttcttc                                                                                                                   >  1:260765/1‑80 (MQ=255)                                                                       aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTg                                                                                                                              <  2:946992/69‑1 (MQ=25)                                                                       aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTg                                                                                                                              >  1:946992/1‑69 (MQ=25)                                                                       aGTATTAAGACTATCACTTAATTAAGTGATATTGGTTTTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTttcttc                                                                                                                   <  2:260765/80‑1 (MQ=17)                                                                        gTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCtt                                                       <  1:289616/139‑1 (MQ=255)                                                                              aGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                       >  1:1214847/1‑101 (MQ=255)                                                                              aGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                       <  2:1214847/101‑1 (MQ=255)                                                                                   aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACa                                                                             >  2:56419/1‑106 (MQ=255)                                                                                   aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACa                                                                             <  1:56419/106‑1 (MQ=255)                                                                                     cACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGttt                                          >  2:875567/1‑139 (MQ=255)                                                                                          aTTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCtt                                     >  1:1351910/1‑139 (MQ=255)                                                                                           tttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTg                                                                                                                              >  2:107260/1‑49 (MQ=255)                                                                                           tttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTg                                                                                                                              <  1:107260/49‑1 (MQ=255)                                                                                               aGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTtct                                                     >  1:208338/1‑118 (MQ=255)                                                                                               aGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTtct                                                     <  2:208338/118‑1 (MQ=255)                                                                                                gTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTg                                                                                                                              <  2:1299336/44‑1 (MQ=255)                                                                                                gTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTg                                                                                                                              >  1:1299336/1‑44 (MQ=255)                                                                                                 tGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTCGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAg                              <  2:1133938/139‑1 (MQ=255)                                                                                                  gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                       <  1:448938/81‑1 (MQ=255)                                                                                                  gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                       >  2:448938/1‑81 (MQ=255)                                                                                                  gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                  <  2:1167606/86‑1 (MQ=255)                                                                                                  gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaa                                                                                  >  1:1167606/1‑86 (MQ=255)                                                                                                     atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCtt                                                                                                                         <  1:565976/44‑1 (MQ=255)                                                                                                     atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCtt                                                                                                                         >  2:565976/1‑44 (MQ=255)                                                                                                     atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGa                                                                                                   >  1:59480/1‑66 (MQ=255)                                                                                                     atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGa                                                                                                   <  2:59480/66‑1 (MQ=255)                                                                                                     atTGGTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATaa                          >  2:405953/1‑139 (MQ=255)                                                                                                         gTTGTCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTg                      >  1:449235/1‑139 (MQ=255)                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaata                                                                                >  1:90705/1‑77 (MQ=255)                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaata                                                                                <  2:90705/77‑1 (MQ=255)                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACa                                                                             >  1:598133/1‑80 (MQ=255)                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACa                                                                             <  2:598133/80‑1 (MQ=255)                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGc                                                         >  2:1715113/1‑100 (MQ=255)                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGc                                                         <  1:1715113/100‑1 (MQ=255)                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTc                                 >  2:742688/1‑124 (MQ=255)                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTc                                 <  1:742688/124‑1 (MQ=255)                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGa                             >  1:1302876/1‑128 (MQ=255)                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGa                             <  2:1302876/128‑1 (MQ=255)                                                                                                             tCTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGa                  >  2:1540152/1‑139 (MQ=255)                                                                                                              cTGGAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGAt                 >  1:1484371/1‑139 (MQ=255)                                                                                                                ggAGATTCAGGGGGCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATtc               <  2:449235/139‑1 (MQ=255)                                                                                                                           gggCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAAt    >  2:950768/1‑139 (MQ=255)                                                                                                                           gggCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAAt    <  2:1484371/139‑1 (MQ=255)                                                                                                                             gCCAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCa  <  1:950768/139‑1 (MQ=255)                                                                                                                              ccAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaat                                                                                 <  2:325356/59‑1 (MQ=255)                                                                                                                              ccAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGataaat                                                                                 >  1:325356/1‑59 (MQ=255)                                                                                                                              ccAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACa                                                                             >  2:581576/1‑63 (MQ=255)                                                                                                                              ccAGTTTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACa                                                                             <  1:581576/63‑1 (MQ=255)                                                                                                                                   |                                                                                                                                     GGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCA  >  NC_000913/2997559‑2997821

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: ‑