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%	180.4 / NA	56	*108* (TAG→TAA)	ppdC	putative prepilin peptidase‑dependent protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/25);  total (31/25)

GTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGG  >  NC_000913/2962310‑2962572                                                                                                                                    |                                                                                                                                    gTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGAc                                                                                                                              >  2:1021075/1‑139 (MQ=255)    cTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGTCTCCTGATTACTGACGAt                                                                                                                           <  1:1021075/139‑1 (MQ=255)                          gTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCAt                                                                                                     <  2:460479/139‑1 (MQ=255)                               aaaGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGc                                                                                                >  1:955228/1‑139 (MQ=255)                                           aGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGcccc                                                                                    <  1:1030713/139‑1 (MQ=255)                                              cgttcgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCgg                                                                                 >  2:713080/1‑139 (MQ=255)                                                ttcgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                               >  2:1073723/1‑139 (MQ=255)                                                          aCAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCg                                                                     <  2:533969/139‑1 (MQ=255)                                                                      aTCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGAt                                                                                                                           <  2:1816426/73‑1 (MQ=255)                                                                      aTCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGAt                                                                                                                           >  1:1816426/1‑73 (MQ=255)                                                                      aTCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGc                                                         >  2:1352927/1‑139 (MQ=255)                                                                              ttCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACAt                                                 <  1:1352927/139‑1 (MQ=255)                                                                               tCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCAt                                                                                                     <  2:754153/86‑1 (MQ=255)                                                                               tCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCAt                                                                                                     >  1:754153/1‑86 (MQ=255)                                                                               tCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATc                                                >  1:308437/1‑139 (MQ=255)                                                                                    cACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCggg                                                                                                                      >  2:993677/1‑64 (MQ=255)                                                                                    cACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCggg                                                                                                                      <  1:993677/64‑1 (MQ=255)                                                                                        tCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                       <  1:1073723/139‑1 (MQ=255)                                                                                        tCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                       >  1:1676954/1‑139 (MQ=255)                                                                                         ccAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGa                                      >  2:941410/1‑139 (MQ=255)                                                                                                 ttGGAATGGTAGACCCTTAACATAGCGGCTCCTGGTTACTGACGATTCggg                                                                                                                      >  1:1274317/1‑51 (MQ=255)                                                                                                 ttGGAATGGTAGACCCTTAACATAGCGGCTCCTGGTTACTGACGATTCggg                                                                                                                      <  2:1274317/51‑1 (MQ=255)                                                                                                 ttGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCACCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCa                              <  2:955228/139‑1 (MQ=255)                                                                                                       tGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGt                        <  2:1676954/139‑1 (MQ=255)                                                                                                        ggTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGtt                       <  1:941410/139‑1 (MQ=255)                                                                                                         gTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTc                                                                                                       <  1:832287/58‑1 (MQ=255)                                                                                                         gTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTc                                                                                                       >  2:832287/1‑58 (MQ=255)                                                                                                            gACCCTTAACATAGCGGCTCCTGATTACTGACGATTCggg                                                                                                                      <  1:166544/40‑1 (MQ=255)                                                                                                            gACCCTTAACATAGCGGCTCCTGATTACTGACGATTCggg                                                                                                                      >  2:166544/1‑40 (MQ=255)                                                                                                               ccTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCa                                                                                                               >  1:991771/1‑44 (MQ=255)                                                                                                               ccTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCa                                                                                                               <  2:991771/44‑1 (MQ=255)                                                                                                                  tAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTtct                                                                                           <  2:667300/61‑1 (MQ=255)                                                                                                                  tAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTtct                                                                                           >  1:667300/1‑61 (MQ=255)                                                                                                                   aaCATAGCGGCTCCTGATTACTGACGATTCGGGCAAt                                                                                                                  >  1:1436395/1‑37 (MQ=255)                                                                                                                   aaCATAGCGGCTCCTGATTACTGACGATTCGGGCAAt                                                                                                                  <  2:1436395/37‑1 (MQ=255)                                                                                                                       tAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCg                                                                     <  1:1776987/78‑1 (MQ=255)                                                                                                                       tAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCg                                                                     >  2:1776987/1‑78 (MQ=255)                                                                                                                        aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGcccc                                                                                    >  1:162785/1‑62 (MQ=255)                                                                                                                        aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGcccc                                                                                    <  2:162785/62‑1 (MQ=255)                                                                                                                        aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACAt                                                 <  1:1528357/97‑1 (MQ=255)                                                                                                                        aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACAt                                                 >  2:1528357/1‑97 (MQ=255)                                                                                                                        aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                       >  2:1633334/1‑107 (MQ=255)                                                                                                                        aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                       <  1:1633334/107‑1 (MQ=255)                                                                                                                        aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCa       >  1:682148/1‑139 (MQ=255)                                                                                                                         gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGTGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAg      >  1:527051/1‑139 (MQ=255)                                                                                                                         gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTtct                                                                                           <  2:120331/54‑1 (MQ=255)                                                                                                                         gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTtct                                                                                           >  1:120331/1‑54 (MQ=255)                                                                                                                         gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGa                                                      >  2:800909/1‑91 (MQ=255)                                                                                                                         gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGa                                                      <  1:800909/91‑1 (MQ=255)                                                                                                                         gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAg      >  2:1337839/1‑139 (MQ=255)                                                                                                                         gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAg      >  1:527052/1‑139 (MQ=255)                                                                                                                          cGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAgg     >  2:497138/1‑139 (MQ=255)                                                                                                                             cTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTgg  >  2:1115885/1‑139 (MQ=255)                                                                                                                             cTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTgg  >  1:1207438/1‑139 (MQ=255)                                                                                                                                cTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGa                                                      >  2:411889/1‑84 (MQ=255)                                                                                                                                cTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGa                                                      <  1:411889/84‑1 (MQ=255)                                                                                                                                    |                                                                                                                                    GTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGG  >  NC_000913/2962310‑2962572

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