Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	115,661	T→C	M22V (ATG→GTG)	hofC ←	assembly protein in type IV pilin biogenesis, transmembrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	115,661	0	T	C	100.0%	118.2 / NA	33	M22V (ATG→GTG)	hofC	assembly protein in type IV pilin biogenesis, transmembrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/13);  total (20/13)

TTTATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGG  >  NC_000913/115548‑115781                                                                                                                  |                                                                                                                         tttATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGc                                                                                                   <  2:725366/137‑1 (MQ=255) tttATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGc                                                                                                   >  1:725366/1‑137 (MQ=255)          tcGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTatgcca                                                                                        >  2:432170/1‑139 (MQ=255)                   cgcAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAgcgc                                                                               >  2:411253/1‑139 (MQ=255)                   cgcAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAgcgc                                                                               <  2:482576/139‑1 (MQ=255)                             aTGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGTACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCt                                                                     >  1:358098/1‑139 (MQ=255)                               ggCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTa                                                                   >  1:191912/1‑139 (MQ=255)                                gCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTAc                                                                  <  1:251132/139‑1 (MQ=255)                                 cGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACt                                                                 <  2:639176/139‑1 (MQ=255)                                  gATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTc                                                                >  1:702650/1‑139 (MQ=255)                                         ttCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGTACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGc                                                         <  2:358098/139‑1 (MQ=255)                                         ttCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGc                                                         >  2:230142/1‑139 (MQ=255)                                                   gTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCa                                               >  1:784844/1‑139 (MQ=255)                                                        gTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCACCATGCGGCATCCCCAGTAcg                                          >  1:760092/1‑139 (MQ=255)                                                                cTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTaa                                  <  1:522311/139‑1 (MQ=255)                                                                  gTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAAct                                >  2:160230/1‑139 (MQ=255)                                                                   ttGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAActc                               <  2:702650/139‑1 (MQ=255)                                                                         tAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTa                                            <  1:152560/120‑1 (MQ=255)                                                                         tAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTa                                            >  2:152560/1‑120 (MQ=255)                                                                              ccATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAggtgg                    >  2:560059/1‑139 (MQ=255)                                                                              ccATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCACCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAggtgg                    <  2:760092/139‑1 (MQ=255)                                                                                aTAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGtt                  <  1:160230/139‑1 (MQ=255)                                                                                   aGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAg               >  1:440096/1‑139 (MQ=255)                                                                                        tAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTg          >  1:626279/1‑139 (MQ=255)                                                                                             tACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTa     >  1:682830/1‑139 (MQ=255)                                                                                             tACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTa     >  2:605062/1‑139 (MQ=255)                                                                                                ggCTCTCTGCCCATAGCACCCCATCTTGCGCATTGc                                                                                                         >  2:381467/1‑36 (MQ=255)                                                                                                ggCTCTCTGCCCATAGCACCCCATCTTGCGCATTGc                                                                                                         <  1:381467/36‑1 (MQ=255)                                                                                                ggCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACgg  >  1:784689/1‑139 (MQ=255)                                                                                                  ctctctGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTa                                                                   <  2:575657/72‑1 (MQ=255)                                                                                                  ctctctGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTa                                                                   >  1:575657/1‑72 (MQ=255)                                                                                                       tGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCa                           <  1:576268/107‑1 (MQ=255)                                                                                                       tGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCa                           >  2:576268/1‑107 (MQ=255)                                                                                                                  |                                                                                                                         TTTATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGG  >  NC_000913/115548‑115781

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