Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	197,199	T→C	R218R (CGT→CGC)	rseP →	intramembrane zinc metalloprotease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	197,199	0	T	C	100.0%	168.8 / NA	48	R218R (CGT→CGC)	rseP	intramembrane zinc metalloprotease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/13);  total (35/13)

GAAAGCACCACCATTACAGTAGCGCCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTT  >  NC_000913/197071‑197334                                                                                                                                 |                                                                                                                                        gAAAGCACCACCATTACAGTAGCGCCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACGGGGGGTTTGAGCCTGATAAACAAAAGCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAAt                                                                                                                              >  2:501536/1‑140 (MQ=255)    aGCACCACCATTACAGTAGCGCCATTTGGCAGCGACCAACGGCGTGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGa                                                                                                                           <  2:663492/140‑1 (MQ=255)      cagcaccaTTACAGTAGCGCCATTTGGCGGCGCCCAACGGCGGGATGTAAAGCTCGTTTTCCGTCCCGGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCCCGGGGGCCGCAAATTGAAc                                                                                                                         <  2:905058/137‑1 (MQ=255)      caccacCATTACAGTAGCGCCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAAc                                                                                                                         >  2:452482/1‑140 (MQ=255)            caTTACAGTAGCGCCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTAc                                                                                                                   >  1:709563/1‑140 (MQ=255)              ttACAGTAGCGCCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTg                                                                                                                 >  2:953105/1‑140 (MQ=255)                    tAGCGCCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGaactg                                                                                                           >  1:905058/1‑137 (MQ=255)                      gcgcCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAAc                                                                                                                         <  1:990485/124‑1 (MQ=255)                      gcgcCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAAc                                                                                                                         >  2:990485/1‑124 (MQ=255)                      gcgcCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCt                                                                                                                       <  1:146048/126‑1 (MQ=255)                      gcgcCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCt                                                                                                                       >  2:146048/1‑126 (MQ=255)                      gcgcCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAAtgt                                                                                                         >  2:293430/1‑140 (MQ=255)                      gcgcCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAAtgt                                                                                                         <  2:830564/140‑1 (MQ=255)                      gcgcCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAAtgt                                                                                                         >  1:131489/1‑140 (MQ=255)                      gcgcCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAAtgt                                                                                                         >  1:830564/1‑140 (MQ=255)                            tttGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGcc                                                                                                   >  1:47214/1‑140 (MQ=255)                                               gggATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGcaagcaag                                                                                >  2:809778/1‑140 (MQ=255)                                                         gCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGc                                                                      <  1:809778/140‑1 (MQ=255)                                                          cTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCa                                                                     >  2:314592/1‑140 (MQ=255)                                                                         tGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGAt                                                      >  2:591448/1‑140 (MQ=255)                                                                         tGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGAt                                                      >  1:635948/1‑140 (MQ=255)                                                                           ggCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCg                                                    >  1:298321/1‑140 (MQ=255)                                                                             cGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGtt                                                  <  1:953105/140‑1 (MQ=255)                                                                               tttGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTaa                                                >  2:31379/1‑140 (MQ=255)                                                                                ttGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTaaa                                               >  2:211404/1‑140 (MQ=255)                                                                                   aGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTc                                            >  1:264122/1‑140 (MQ=255)                                                                                   aGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTc                                            >  1:310986/1‑140 (MQ=255)                                                                                    gCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCg                                           <  2:106659/140‑1 (MQ=255)                                                                                         aTAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGt                                      <  2:298321/140‑1 (MQ=255)                                                                                           aaaGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCa                                    <  1:314592/140‑1 (MQ=255)                                                                                                   tCCGGTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTaa                            >  1:957302/1‑140 (MQ=255)                                                                                                      ggTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGc                         <  2:131489/140‑1 (MQ=255)                                                                                                       gTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCa                                    >  1:843263/1‑128 (MQ=255)                                                                                                       gTATCTTCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCa                                    <  2:843263/128‑1 (MQ=255)                                                                                                            ttCGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTggg                   >  1:563243/1‑140 (MQ=255)                                                                                                              cGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTg                 >  1:527328/1‑140 (MQ=255)                                                                                                              cGCTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTg                 >  2:534217/1‑140 (MQ=255)                                                                                                                cTGGGGATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGAc               <  2:957302/140‑1 (MQ=255)                                                                                                                    ggATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGt                  >  1:965314/1‑133 (MQ=255)                                                                                                                    ggATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCttt           >  1:680277/1‑140 (MQ=255)                                                                                                                    ggATTCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCttt           >  2:294869/1‑140 (MQ=255)                                                                                                                    ggATTCGTCCTCGCGGCCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGt                  <  2:965314/133‑1 (MQ=255)                                                                                                                    ggATTAGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCttt           >  1:389000/1‑140 (MQ=255)                                                                                                                       ttCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTtgtg        >  1:961478/1‑140 (MQ=255)                                                                                                                        tCGTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCATGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGa       >  1:962603/1‑140 (MQ=255)                                                                                                                          gTCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATg     >  2:742337/1‑140 (MQ=255)                                                                                                                           tCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGc    <  2:680277/140‑1 (MQ=255)                                                                                                                           tCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGc    >  1:40179/1‑140 (MQ=255)                                                                                                                           tCCTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGc    >  2:41289/1‑140 (MQ=255)                                                                                                                             cTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAGGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCtt  >  1:535094/1‑140 (MQ=255)                                                                                                                             cTCGCGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGACAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCtt  >  1:780340/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                        GAAAGCACCACCATTACAGTAGCGCCATTTGGCAGCGACCAACGGCGGGATGTAAAGCTCGATTTACGTCACTGGGCGTTTGAGCCTGATAAAGAAGATCCGGTATCTTCGCTGGGGATTCGTCCTCGTGGGCCGCAAATTGAACCTGTACTGGAAAATGTGCAGCCAAACTCGGCGGCAAGCAAGGCAGGTTTGCAAGCAGGCGACAGGATCGTTAAAGTCGATGGTCAGCCCTTAACGCAGTGGGTGACCTTTGTGATGCTT  >  NC_000913/197071‑197334

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