Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,493,505	A→G	S346G (AGC→GGC)	trg →	methyl‑accepting chemotaxis protein Trg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,493,505	0	A	G	100.0%	124.9 / NA	36	S346G (AGC→GGC)	trg	methyl‑accepting chemotaxis protein Trg
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/17);  total (19/17)

CACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCT  >  NC_000913/1493375‑1493639                                                                                                                                   |                                                                                                                                       cACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCg                                                                                                                               >  1:160386/1‑140 (MQ=255)              gcggcgGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTc                                                                                                                  <  2:1120743/140‑1 (MQ=255)                ggcggcTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCg                                                                                                                <  2:996736/140‑1 (MQ=255)                 gcggcTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCgg                                                                                                               <  1:365541/140‑1 (MQ=255)                 gcggcTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCgg                                                                                                               >  2:821129/1‑140 (MQ=255)                 gcggcTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCgg                                                                                                               <  1:1061444/140‑1 (MQ=255)                   ggcTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGgtg                                                                                                             <  1:895768/140‑1 (MQ=255)                       aTCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTgtagt                                                                                                         <  1:58109/140‑1 (MQ=255)                                aCTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGAt                                                                                                <  2:160386/140‑1 (MQ=255)                                     cgccAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGgcg                                                                                           <  1:340417/140‑1 (MQ=255)                                                aGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACg                                                                                >  2:967881/1‑140 (MQ=255)                                                aGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACg                                                                                >  2:231843/1‑140 (MQ=255)                                                     cTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTc                                                                           >  2:696495/1‑140 (MQ=255)                                                             gACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGaaaa                                                                   >  2:1044549/1‑140 (MQ=255)                                                                  tGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCt                                                              >  1:389164/1‑140 (MQ=255)                                                                       cAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCg                                                                          <  2:942056/123‑1 (MQ=255)                                                                       cAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCg                                                                          >  1:942056/1‑123 (MQ=255)                                                                          aaTGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCAc                                                      >  1:512738/1‑140 (MQ=255)                                                                            tGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCAccg                                                    <  1:1044549/140‑1 (MQ=255)                                                                            tGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCAccg                                                    >  1:445339/1‑140 (MQ=255)                                                                                    aCGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGtcatc                                            >  2:571013/1‑140 (MQ=255)                                                                                    aCGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGtcatc                                            <  2:445339/140‑1 (MQ=255)                                                                                      gcgcATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCaa                                          <  1:231843/140‑1 (MQ=255)                                                                                         catcatGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAg                                       >  2:190010/1‑140 (MQ=255)                                                                                          atcatGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGt                                      >  1:547764/1‑140 (MQ=255)                                                                                                   gCAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTc                             <  2:512738/140‑1 (MQ=255)                                                                                                    cAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTcc                            >  2:436203/1‑140 (MQ=255)                                                                                                    cAAACTGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTcc                            >  1:343055/1‑140 (MQ=255)                                                                                                         tGGCGCAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACg                       >  2:573335/1‑140 (MQ=255)                                                                                                           gcgcAAGAGGCTTCTATTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGaa                     <  2:343055/140‑1 (MQ=255)                                                                                                                          aTTAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGaa      <  1:190010/140‑1 (MQ=255)                                                                                                                           ttAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAAt     >  2:469216/1‑140 (MQ=255)                                                                                                                            tAAAGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAAtg    <  1:469216/140‑1 (MQ=255)                                                                                                                             aaaGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAAtgc   >  2:695243/1‑140 (MQ=255)                                                                                                                              aaGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAAtgct  >  2:893871/1‑140 (MQ=255)                                                                                                                              aaGCCGGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAAtgct  <  1:695243/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       CACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCT  >  NC_000913/1493375‑1493639

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