Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,963,498	G→A	T73T (ACC→ACT)	ppdB ←	putative prepilin peptidase‑dependent protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,963,498	0	G	A	100.0%	110.2 / NA	32	T73T (ACC→ACT)	ppdB	putative prepilin peptidase‑dependent protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/12);  total (20/12)

TGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCGGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCG  >  NC_000913/2963370‑2963632                                                                                                                                 |                                                                                                                                       tGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGc                                                                                                                              >  1:362632/1‑139 (MQ=255)        tCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACa                                                                                                                       >  1:628374/1‑139 (MQ=255)            aCGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACa                                                                                                                                   >  1:297593/1‑123 (MQ=255)            aCGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACa                                                                                                                                   <  2:297593/123‑1 (MQ=255)                   ccAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGccc                                                                                                            <  1:628479/139‑1 (MQ=255)                                    ttttACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTcgcg                                                                                           <  2:628374/139‑1 (MQ=255)                                      ttACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGAc                                                                                         <  2:255092/139‑1 (MQ=255)                                              tCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATa                                                                                 >  2:89047/1‑139 (MQ=255)                                                              cGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTc                                                                 <  2:345335/139‑1 (MQ=255)                                                              cGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTc                                                                 >  2:335291/1‑139 (MQ=255)                                                              cGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTc                                                                 >  1:345335/1‑139 (MQ=255)                                                                ttACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGa                                                                                          >  2:492588/1‑112 (MQ=255)                                                                ttACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGa                                                                                          <  1:492588/112‑1 (MQ=255)                                                                ttACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCAt                                                               >  1:152393/1‑139 (MQ=255)                                                                      ttCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCa                                                         >  2:169982/1‑139 (MQ=255)                                                                      ttCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCa                                                         >  2:346747/1‑139 (MQ=255)                                                                        cGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCATCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAgc                                                       >  1:519846/1‑139 (MQ=255)                                                                        cGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAAcc                                                                                                                 >  1:581327/1‑81 (MQ=255)                                                                        cGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAAcc                                                                                                                 <  2:581327/81‑1 (MQ=255)                                                                          cATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAgcgc                                                     <  2:519846/139‑1 (MQ=255)                                                                           aTCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAgcgc                                                     >  2:528179/1‑138 (MQ=255)                                                                           aTCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAgcgc                                                     <  1:528179/138‑1 (MQ=255)                                                                                aCTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCtt                                               >  2:216949/1‑139 (MQ=255)                                                                                  tGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTAc                                             >  2:29924/1‑139 (MQ=255)                                                                                       aaTGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTg                                        <  1:335291/139‑1 (MQ=255)                                                                                                            gACAATTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGt                   >  2:183786/1‑139 (MQ=255)                                                                                                                aTTTCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACg               >  2:582696/1‑139 (MQ=255)                                                                                                                 tttCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACgc              <  1:346747/139‑1 (MQ=255)                                                                                                                   tCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCg                                                                             >  1:603737/1‑74 (MQ=255)                                                                                                                   tCCAGCCCTTCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCg                                                                             <  2:603737/74‑1 (MQ=255)                                                                                                                            tCGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGc   >  1:596149/1‑139 (MQ=255)                                                                                                                             cGCCAGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCg  >  2:1742/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                       TGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCGGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCG  >  NC_000913/2963370‑2963632

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