Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,962,441	C→T	*108* (TAG→TAA)	ppdC ←	putative prepilin peptidase‑dependent protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,962,441	0	C	T	100.0%	137.3 / NA	43	*108* (TAG→TAA)	ppdC	putative prepilin peptidase‑dependent protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/19);  total (24/19)

TACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCG  >  NC_000913/2962305‑2962574                                                                                                                                         |                                                                                                                                      tACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGatta                                                                                                                                     >  2:487894/1‑139 (MQ=255)      gTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGAc                                                                                                                                >  2:455708/1‑139 (MQ=255)      gTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGAc                                                                                                                                >  2:523030/1‑139 (MQ=255)            tGCACCAGAATCATCTCGGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCg                                                                                                                          >  2:603878/1‑139 (MQ=255)              cACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCggg                                                                                                                        <  2:822822/139‑1 (MQ=255)              cACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCggg                                                                                                                        >  1:655922/1‑139 (MQ=255)                               gTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCAt                                                                                                       >  2:938723/1‑139 (MQ=255)                                    aaaGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGc                                                                                                  >  2:112122/1‑139 (MQ=255)                                                aGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGcccc                                                                                      <  1:603878/139‑1 (MQ=255)                                                     ttcgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                                 <  1:1016050/139‑1 (MQ=255)                                                     ttcgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                                 >  2:794406/1‑139 (MQ=255)                                                       cgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGa                                                                               >  2:342711/1‑139 (MQ=255)                                                                cAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGt                                                                      >  2:1994/1‑139 (MQ=255)                                                                    aaaCTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACg                                                                  <  1:206872/139‑1 (MQ=255)                                                                      aCTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCt                                                                <  1:455708/139‑1 (MQ=255)                                                                          cATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATg                                                            <  1:966804/139‑1 (MQ=255)                                                                          cATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATg                                                            >  2:1022506/1‑139 (MQ=255)                                                                           aTCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGc                                                           >  1:940715/1‑139 (MQ=255)                                                                               aCGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGattact                                                                                                                                   >  1:301921/1‑63 (MQ=255)                                                                               aCGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGattact                                                                                                                                   <  2:301921/63‑1 (MQ=255)                                                                                 gCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCggg                                                                                  <  2:718539/110‑1 (MQ=255)                                                                                 gCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCggg                                                                                  >  1:718539/1‑110 (MQ=255)                                                                                 gCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTAGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGacac                                                     >  1:630326/1‑139 (MQ=255)                                                                                   ttCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACAt                                                   <  1:938723/139‑1 (MQ=255)                                                                                        gCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGc                                              <  2:940715/139‑1 (MQ=255)                                                                                        gCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGc                                              <  1:246265/139‑1 (MQ=255)                                                                                             tCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                         <  2:630326/139‑1 (MQ=255)                                                                                             tCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                         <  2:1020805/139‑1 (MQ=255)                                                                                                 gACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGt                                     <  1:14818/139‑1 (MQ=255)                                                                                                      ttGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCa                                <  1:62573/139‑1 (MQ=255)                                                                                                       tGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCAt                               >  1:704695/1‑139 (MQ=255)                                                                                                            tGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGt                          <  1:571945/139‑1 (MQ=255)                                                                                                             ggTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGtt                         <  1:798785/139‑1 (MQ=255)                                                                                                                    ccTTAACATAGCGGCTCCTTATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGc                  <  1:472619/139‑1 (MQ=255)                                                                                                                        aaCATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGt              >  1:80654/1‑139 (MQ=255)                                                                                                                           aTAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTgg           <  1:794406/139‑1 (MQ=255)                                                                                                                             aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCa         >  2:471318/1‑139 (MQ=255)                                                                                                                              gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACg                                                                  <  2:224571/81‑1 (MQ=255)                                                                                                                              gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACg                                                                  >  1:224571/1‑81 (MQ=255)                                                                                                                              gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAg        >  2:326407/1‑139 (MQ=255)                                                                                                                              gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAg        >  1:83271/1‑139 (MQ=255)                                                                                                                              gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAg        >  1:883729/1‑139 (MQ=255)                                                                                                                                ggCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGt      >  2:445276/1‑139 (MQ=255)                                                                                                                                  cTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTgg    >  2:445285/1‑139 (MQ=255)                                                                                                                                   tCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGc   >  1:900989/1‑139 (MQ=255)                                                                                                                                    ccTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCg  >  2:317910/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                      TACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCG  >  NC_000913/2962305‑2962574

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