Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,113,662	C→T	A231T (GCC→ACC)	pppA ←	bifunctional prepilin leader peptidase/ methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,113,662	0	C	T	100.0%	176.0 / NA	56	A231T (GCC→ACC)	pppA	bifunctional prepilin leader peptidase/ methylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/25);  total (31/25)

GGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGA  >  NC_000913/3113535‑3113796                                                                                                                                |                                                                                                                                       ggTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAg                                                                                                                                >  1:849936/1‑136 (MQ=255) ggTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAg                                                                                                                                <  2:849936/136‑1 (MQ=255)             aaCAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCg                                                                                                                 >  2:722935/1‑139 (MQ=255)             aaCAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCg                                                                                                                 >  1:443420/1‑139 (MQ=255)                  tGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAGCATTgggaaacgcaaat                                                                                                            >  1:1334922/1‑129 (MQ=255)                  tGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAAt                                                                                                            >  2:1710134/1‑139 (MQ=255)                   gCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATg                                                                                                           >  1:874774/1‑139 (MQ=255)                    ccTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGc                                                                                                          <  1:1710134/139‑1 (MQ=255)                     cTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGcc                                                                                                         <  2:443420/139‑1 (MQ=255)                     cTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGcc                                                                                                         <  2:1631793/139‑1 (MQ=255)                        tAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAg                                                                                                                                >  1:1358594/1‑113 (MQ=255)                        tAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAg                                                                                                                                <  2:1358594/113‑1 (MQ=255)                             aaaGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCcccaccca                                                                                                 >  2:272679/1‑139 (MQ=255)                               aGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGc                                                                                               >  2:938968/1‑139 (MQ=255)                                gtgtTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGcc                                                                                              >  2:687888/1‑139 (MQ=255)                                  gtTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCAc                                                                                            >  2:21486/1‑139 (MQ=255)                                    tGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCt                                                                                          >  2:1489575/1‑139 (MQ=255)                                     gCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTa                                                                                         >  2:313224/1‑139 (MQ=255)                                         taCCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACa                                                                                                              >  1:465435/1‑114 (MQ=255)                                         taCCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACa                                                                                                              <  2:465435/114‑1 (MQ=255)                                             gCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACgcggcg                                                                                 >  1:1137551/1‑139 (MQ=255)                                             gCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACgcggcg                                                                                 <  2:1137551/139‑1 (MQ=255)                                             gCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACgcggcg                                                                                 <  1:313224/139‑1 (MQ=255)                                               ccAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGaa                                                                               >  1:780600/1‑139 (MQ=255)                                               ccAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGaa                                                                               <  2:780600/139‑1 (MQ=255)                                                   aCTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGAt                                                                                                                                     <  2:1163861/81‑1 (MQ=255)                                                   aCTTAAACACGGTCCAAAAGGCAGTGTGATTGATACTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGAt                                                                                                                                     >  1:1163861/1‑81 (MQ=255)                                                        aaCACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATAc                                                                      >  2:1177141/1‑139 (MQ=255)                                                         acacGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACa                                                                     <  1:1124676/139‑1 (MQ=255)                                                         acacGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACa                                                                     <  1:938968/139‑1 (MQ=255)                                                             ggTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGc                                                                 >  1:530944/1‑139 (MQ=255)                                                              gTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGcc                                                                <  2:530944/139‑1 (MQ=255)                                                               tCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGccc                                                               <  1:917488/139‑1 (MQ=255)                                                                      ggCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTa                                                        <  1:21486/139‑1 (MQ=255)                                                                         aGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATg                                                     <  1:1177141/139‑1 (MQ=255)                                                                          gtgtGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACAt                                                                                                                          >  1:208506/1‑69 (MQ=255)                                                                          gtgtGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACAt                                                                                                                          <  2:208506/69‑1 (MQ=255)                                                                           tgtgGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGct                                                   >  1:77547/1‑139 (MQ=255)                                                                             tgGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACgcggcg                                                                                 >  2:960482/1‑107 (MQ=255)                                                                             tgGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACgcggcg                                                                                 <  1:960482/107‑1 (MQ=255)                                                                              ggTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGcttct                                                <  1:146376/139‑1 (MQ=255)                                                                                        ctTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAg                                                                                                                                <  1:424544/49‑1 (MQ=255)                                                                                        ctTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAg                                                                                                                                >  2:424544/1‑49 (MQ=255)                                                                                          ttttGTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACt                                    >  2:15129/1‑139 (MQ=255)                                                                                              gTAATAACGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTc                                >  1:136098/1‑139 (MQ=255)                                                                                                   aaCGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCt                           <  2:1334922/139‑1 (MQ=255)                                                                                                     cGGCATATATCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTAt                         >  1:569286/1‑139 (MQ=255)                                                                                                             atCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCa                 <  1:15129/139‑1 (MQ=255)                                                                                                              tCAGGCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCag                <  2:874774/139‑1 (MQ=255)                                                                                                                 ggCCGCAGCATGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCagaga             >  2:71792/1‑139 (MQ=255)                                                                                                                      cagcaTGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCccca                                                                                                     >  2:738140/1‑46 (MQ=255)                                                                                                                      cagcaTGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCccca                                                                                                     <  1:738140/46‑1 (MQ=255)                                                                                                                      cagcaTGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTaaa        >  2:1439536/1‑139 (MQ=255)                                                                                                                      cagcaTGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTaaa        >  1:1611471/1‑139 (MQ=255)                                                                                                                       agcaTGAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTaaaa       <  1:272679/139‑1 (MQ=255)                                                                                                                            gAGGTGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGa  >  2:1567341/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                       GGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGA  >  NC_000913/3113535‑3113796

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