Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,465,938	A→T	V132V (GTA→GTT)	gspO →	bifunctional prepilin leader peptidase/ methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,465,938	0	A	T	100.0%	149.9 / NA	47	V132V (GTA→GTT)	gspO	bifunctional prepilin leader peptidase/ methylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (24/23);  total (24/23)

ACGCTCAGTGTTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTAGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCAATGA  >  NC_000913/3465807‑3466071                                                                                                                                    |                                                                                                                                      aCGCTCAGTGTTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCt                                                                                                                                <  2:664414/139‑1 (MQ=255) aCGCTCAGTGTTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCCTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGTGGTTTACGGCGCGGTTGCGGGCt                                                                                                                                >  2:554769/1‑139 (MQ=255)      cAGTGTTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATggg                                                                                                                           <  2:266910/139‑1 (MQ=255)       aGTGTTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATgggg                                                                                                                          >  2:337889/1‑139 (MQ=255)       aGTGTTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATgggg                                                                                                                          >  2:244302/1‑139 (MQ=255)         tgtTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTg                                                                                                                        <  2:16677/139‑1 (MQ=255)              ttGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCtgtg                                                                                                                   <  2:140720/139‑1 (MQ=255)                   tttCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATg                                                                                                                             >  1:144190/1‑124 (MQ=255)                   tttCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATg                                                                                                                             <  2:144190/124‑1 (MQ=255)                           tCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTgg                                                                                                      <  1:244302/139‑1 (MQ=255)                             aGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTgggg                                                                                                    <  1:309451/139‑1 (MQ=255)                               gttcTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGGTTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCg                                                                                                  <  1:632727/136‑1 (MQ=255)                                     cccGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGt                                                                                            <  1:746084/139‑1 (MQ=255)                                                         gctgctCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGa                                                                        <  2:711269/139‑1 (MQ=255)                                                           tgctCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACaaa                                                                              <  1:158590/131‑1 (MQ=255)                                                           tgctCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACaaa                                                                              >  2:158590/1‑131 (MQ=255)                                                             ctctGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGTCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTgg                                                                    >  2:551169/1‑139 (MQ=255)                                                             ctctGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAaga                                                                            >  2:394668/1‑131 (MQ=255)                                                             ctctGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAaga                                                                            <  1:394668/131‑1 (MQ=255)                                                                   cttggcttggTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACg                                                              <  2:277672/139‑1 (MQ=255)                                                                            gTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCa                                                     >  1:363858/1‑139 (MQ=255)                                                                                  aaTGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTAc                                               >  2:676272/1‑139 (MQ=255)                                                                                    tGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTg                                             >  1:606779/1‑139 (MQ=255)                                                                                      cgcAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATggg                                                                                                                           <  1:31922/59‑1 (MQ=255)                                                                                      cgcAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATggg                                                                                                                           >  2:31922/1‑59 (MQ=255)                                                                                      cgcAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGc                                           <  1:676272/139‑1 (MQ=255)                                                                                       gcAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATg                                                                                                                             <  1:341289/56‑1 (MQ=255)                                                                                       gcAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATg                                                                                                                             >  2:341289/1‑56 (MQ=255)                                                                                                aTTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCggg                                                                                                                                  <  2:556101/42‑1 (MQ=255)                                                                                                aTTAATTGATTTACATGATGCGGTTTACGGCGCGGTTGCggg                                                                                                                                  >  1:556101/1‑42 (MQ=255)                                                                                                     ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGttt                                                                                                           <  1:110743/60‑1 (MQ=255)                                                                                                     ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGttt                                                                                                           >  2:110743/1‑60 (MQ=255)                                                                                                     ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGc                                           <  2:93025/124‑1 (MQ=255)                                                                                                     ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGc                                           >  1:93025/1‑124 (MQ=255)                                                                                                     ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATg                            >  1:483181/1‑139 (MQ=255)                                                                                                     ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATg                            >  2:687485/1‑139 (MQ=255)                                                                                                     ttGATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATg                            >  2:352802/1‑139 (MQ=255)                                                                                                       gATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGgt                          <  2:483181/139‑1 (MQ=255)                                                                                                       gATTTACATGATGCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGgt                          <  2:606779/139‑1 (MQ=255)                                                                                                          ttACATGATGCGGTTGACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGtgtg                       <  2:410710/139‑1 (MQ=255)                                                                                                                  tgCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCaa               >  2:647702/1‑139 (MQ=255)                                                                                                                  tgCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCaa               >  2:676812/1‑139 (MQ=255)                                                                                                                  tgCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCaa               >  2:302238/1‑139 (MQ=255)                                                                                                                   gCGGTTTACGGCGCGGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGATTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCaaa              >  2:248466/1‑139 (MQ=255)                                                                                                                               gcgGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCAATGa  >  2:633615/1‑139 (MQ=255)                                                                                                                                cgGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTa                                                <  1:93424/92‑1 (MQ=255)                                                                                                                                cgGTTGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTa                                                >  2:93424/1‑92 (MQ=255)                                                                                                                                    |                                                                                                                                      ACGCTCAGTGTTATTGATTTTCGCACTCAGCTCCTGCCCGACAAACTCACCTTACCGCTGCTCTGGCTTGGCTTGGTATTTAATGCGCAGTATGGATTAATTGATTTACATGATGCGGTTTACGGCGCGGTAGCGGGCTATGGGGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCAATGA  >  NC_000913/3465807‑3466071

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