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

ACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAG  >  NC_000913/2962312‑2962568                                                                                                                                  |                                                                                                                                aCTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGa                                                                                                                        >  2:272542/1‑139 (MQ=255)       cACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCggg                                                                                                                  >  1:416714/1‑139 (MQ=255)            gAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTg                                                                                                                            <  1:12711/124‑1 (MQ=255)            gAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTg                                                                                                                            >  2:12711/1‑124 (MQ=255)             aaTCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGc                                                                                                            >  1:580960/1‑139 (MQ=255)               tcatcTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCa                                                                                                           >  2:607215/1‑138 (MQ=255)               tcatcTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCa                                                                                                           <  1:607215/138‑1 (MQ=255)                                  aTCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTtct                                                                                       >  1:206505/1‑139 (MQ=255)                                         aGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGcccc                                                                                <  1:88834/139‑1 (MQ=255)                                              ttcgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                           >  1:552868/1‑139 (MQ=255)                                                       gACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGc                                                                  >  1:296343/1‑139 (MQ=255)                                                                   cATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCg                                                                              >  2:65521/1‑115 (MQ=255)                                                                   cATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCg                                                                              <  1:65521/115‑1 (MQ=255)                                                                   cATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATg                                                      >  2:727842/1‑139 (MQ=255)                                                                   cATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATg                                                      <  2:53571/139‑1 (MQ=255)                                                                              ccAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATcc                                           <  2:552868/139‑1 (MQ=255)                                                                                aGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCg                                         >  1:456568/1‑139 (MQ=255)                                                                                aGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCg                                         <  2:580960/139‑1 (MQ=255)                                                                                     gTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGc                                    <  1:272542/139‑1 (MQ=255)                                                                                      ttcacaCGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTAACGATGCTGGCAATGCAGGCGGGTCATCCCGCCTGCTCTGCCCCCGGGTGAAACTAGCGTAACGCGGATGCTGACACCTCCCGCCTGCg                                   <  2:124958/134‑1 (MQ=255)                                                                                       ccAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGa                                  <  2:676555/139‑1 (MQ=255)                                                                                          gACGATTGGAATGGTAGACGCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATTCTGACACATCCCGCCTGCGATGt                               >  1:514371/1‑139 (MQ=255)                                                                                               ttGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCa                          >  2:189671/1‑139 (MQ=255)                                                                                                             ccTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                           <  2:618616/76‑1 (MQ=255)                                                                                                             ccTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                           >  1:618616/1‑76 (MQ=255)                                                                                                             ccTTAACATAGCGGCTCCTGATTACCGACGATTCGGGCAATGCCGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGc            <  2:334936/139‑1 (MQ=255)                                                                                                                 aaCATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGt        >  1:50702/1‑139 (MQ=255)                                                                                                                 aaCATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGt        >  2:231270/1‑139 (MQ=255)                                                                                                                     tAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGc    <  2:456568/139‑1 (MQ=255)                                                                                                                      aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCa   >  2:362086/1‑139 (MQ=255)                                                                                                                       gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTctct                                                                                     <  1:217286/56‑1 (MQ=255)                                                                                                                       gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTctct                                                                                     >  2:217286/1‑56 (MQ=255)                                                                                                                       gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATg                                                      <  2:383640/87‑1 (MQ=255)                                                                                                                       gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATg                                                      >  1:383640/1‑87 (MQ=255)                                                                                                                       gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAg  >  2:693589/1‑139 (MQ=255)                                                                                                                       gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATACGGTTGACCTGCCAGTTGGCAg  >  1:146960/1‑139 (MQ=255)                                                                                                                              cTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGc            <  1:345233/122‑1 (MQ=255)                                                                                                                              cTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGc            >  2:345233/1‑122 (MQ=255)                                                                                                                                  |                                                                                                                                ACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAG  >  NC_000913/2962312‑2962568

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