Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,465,102	A→G	A7A (GCA→GCG)	gspM →	Type II secretion system protein GspM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,465,102	0	A	G	100.0%	151.6 / NA	44	A7A (GCA→GCG)	gspM	Type II secretion system protein GspM
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/20);  total (24/20)

GAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAA  >  NC_000913/3464975‑3465229                                                                                                                                |                                                                                                                                  gAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTTAGGTATAAGATTAAATCATGGCGGGCGGaata                                                                                                                              >  2:302446/1‑131 (MQ=255)   aGCTAATATCGATCGTTTCTGTGAATTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAGGGAACTCAGGGAGATGATTAAATCATGGTGAGCGGAAAAAt                                                                                                                            >  2:168124/1‑133 (MQ=255)       aaTATCGACCGTTTCAGTCAATTAACACAGCCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCcac                                                                                                                        <  1:168124/133‑1 (MQ=255)         tatCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACAt                                                                                                                      <  2:434519/133‑1 (MQ=255)           tcgatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCg                                                                                                                    <  2:3815/133‑1 (MQ=255)            cgatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATcag                                                                                                                   >  2:635812/1‑131 (MQ=255)              atcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGaa                                                                                                                 >  1:39728/1‑133 (MQ=255)              atcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGATCAGTGGTGTATGGACAGTAAGGAACTCAGGAAAATGATTAAATCATGGTGGGCGGAAAAATCAACATCGGaa                                                                                                                 >  2:108874/1‑133 (MQ=255)               tcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGaaa                                                                                                                >  1:100751/1‑133 (MQ=255)                    tCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAg                                                                                                           >  2:611190/1‑133 (MQ=255)                                  caGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATGCACATCGGAAAGACGGATCGTGGCGGCACt                                                                                             >  1:163137/1‑133 (MQ=255)                                    gTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTgg                                                                                           >  2:527464/1‑133 (MQ=255)                                    gTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTgg                                                                                           >  2:422716/1‑133 (MQ=255)                                      ccTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCg                                                                                         >  1:121035/1‑133 (MQ=255)                                      ccTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATAGGAAAAACAGATGGTTGCGGCACTGGCg                                                                                         >  2:270897/1‑133 (MQ=255)                                               cGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGGTCGTTGCGGCACTGGCGGTTCTctcg                                                                                >  1:595884/1‑133 (MQ=255)                                                             gAAAAAGATCCGGTCCGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTttg                                                                 <  1:257714/133‑1 (MQ=255)                                                             gAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTttg                                                                 <  2:589060/133‑1 (MQ=255)                                                             gAAAAAGATCCGGTCAGTAGTGGATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTttg                                                                 <  1:298220/133‑1 (MQ=255)                                                             gAAAAAGAACCGGTCTGTTGTGTATGGACACTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTttg                                                                 <  1:309765/133‑1 (MQ=255)                                                              aaaaaGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCTCGGC‑Ggttttttt                                                                >  2:577809/1‑131 (MQ=255)                                                               aaaaGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCCCTCGGC‑GGtttttttt                                                               >  1:399244/1‑133 (MQ=255)                                                               aaaaGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GGTTTttgtt                                                               >  2:161543/1‑133 (MQ=255)                                                                   gATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑t                                                          <  1:577809/133‑1 (MQ=255)                                                                             gtggtgTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATa                                                <  2:372072/133‑1 (MQ=255)                                                                              tggtgTATGGACAGTCAGTAACACAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATaa                                               <  1:188051/133‑1 (MQ=255)                                                                              tggtgTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATaa                                               <  2:39728/133‑1 (MQ=255)                                                                                gtgtATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCAGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATaaaa                                             <  1:635812/133‑1 (MQ=255)                                                                                   tATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCCCTCGGCGGTTTTTTTTTGGT‑TTGGGGGGTAAAAAcg                                           >  2:7627/1‑132 (MQ=255)                                                                                      ggACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATAAAACCgata                                       <  2:595884/133‑1 (MQ=255)                                                                                        aCAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTTTTTGC‑TTGGGGGGGTAAAAACcataaa                                     >  1:249244/1‑131 (MQ=255)                                                                                         cAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATAAAACCgatagat                                    <  1:592862/133‑1 (MQ=255)                                                                                         cAGTAAGGAACTCAGCGAAATGATCAAACCATTGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATAAAACCgatagat                                    <  2:163137/133‑1 (MQ=255)                                                                                           gTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTTTCCCTCCGG‑GTTTTTTTTTTGG‑TTGGGGGGGTAAAAACgatagataa                                  >  1:475879/1‑132 (MQ=255)                                                                                            tAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GGTTTTTTTTGGG‑TTTGGGGGATAAAAACCATAGATAcc                                 >  1:405210/1‑133 (MQ=255)                                                                                              aGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTTTTTGC‑TTGGGAGGAGAAAAACGAACGATAACTa                               >  1:277617/1‑133 (MQ=255)                                                                                              aGGAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATAAAACCGATAGATACCTa                               <  2:515844/133‑1 (MQ=255)                                                                                                gAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATAAAACCGATAGATACCTACa                             <  1:527464/133‑1 (MQ=255)                                                                                                gAACTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATAAAACCGATAGATACCTACa                             <  2:121035/133‑1 (MQ=255)                                                                                                 aaCTCAGGTAAATGACTACATCTTGGTGGGAGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATAAAACCGATAGATACCTACAt                            <  1:302446/133‑1 (MQ=255)                                                                                                  aCTCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGG‑GGTTTTTTTTTGG‑TTGGGGGGGTGAAAACGAGAGATACATAAAtt                           >  2:5484/1‑133 (MQ=255)                                                                                                      aGGGAAATGATTAAACCATGGTGGCCGGTAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCag                       <  2:277617/133‑1 (MQ=255)                                                                                                                         ggtggGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCGCTCGGC‑GGTTTTTTTTTGC‑TTGGGGGGGTAAAAACCATAAAGACCCACAATGCAGGGCATCAGGCCCACCCGCa    >  2:235425/1‑133 (MQ=255)                                                                                                                         ggtggGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTTTTGGCTTTGGGGGGTAAAAACCGA‑AGATAACCCCATTTTCGAGCATCACTCCCCCgcgcc    >  2:534137/1‑132 (MQ=255)                                                                                                                         ggtggGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCTCGGCGGTTTTTTTTTGGT‑TTGGGG‑GATAAAAACGAAAGATACCTACATTTCAGAGCATCAGTCCCCCGCGCa    >  2:346937/1‑133 (MQ=255)                                                                                                                           tggGCGGATAAATCCGCAGCGGAAAAACAGGTCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCaaa  <  1:270897/133‑1 (MQ=255)                                                                                                                                |                                                                                                                                  GAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGC‑GTTTTTTGTTGGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAA  >  NC_000913/3464975‑3465229

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