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%	177.0 / NA	49	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 (26/23);  total (26/23)

ATGACTTCCGCGCTTTTATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCC  >  NC_000913/115534‑115783                                                                                                                                |                                                                                                                           atgaCTTCCGCGCTTTTATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTgcg                                                                                                                 >  2:139071/1‑139 (MQ=255)  tgtcTTCCGCGCTTTTATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTgcgc                                                                                                                <  1:452982/136‑1 (MQ=255)        ccGCGCTTTTATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGcc                                                                                                          >  2:933277/1‑139 (MQ=255)            gcTTTTATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCg                                                                                                      <  1:933277/139‑1 (MQ=255)            gcTTTTATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCg                                                                                                      <  2:367868/139‑1 (MQ=255)                            cACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTatgccatgcc                                                                                      <  1:694235/139‑1 (MQ=255)                                 cgcAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAgcgc                                                                                 >  2:43603/1‑139 (MQ=255)                                 cgcAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAgcgc                                                                                 >  1:757385/1‑139 (MQ=255)                                                gATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGt                                                                                                 >  2:445138/1‑108 (MQ=255)                                                gATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGt                                                                                                 <  1:445138/108‑1 (MQ=255)                                                    cGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCa                                                              <  1:70148/139‑1 (MQ=255)                                                     gCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCAt                                                             <  1:43603/139‑1 (MQ=255)                                                            gCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCg                                                                                                      <  2:337770/91‑1 (MQ=255)                                                            gCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCg                                                                                                      >  1:337770/1‑91 (MQ=255)                                                                      gTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTatgccatgcc                                                                                      <  2:33101/97‑1 (MQ=255)                                                                      gTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTatgccatgcc                                                                                      >  1:33101/1‑97 (MQ=255)                                                                      gTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTAcg                                            >  1:755907/1‑139 (MQ=255)                                                                            aTCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGAt                                       <  2:918282/138‑1 (MQ=255)                                                                            aTCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGAt                                       >  1:918282/1‑138 (MQ=255)                                                                              cTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTaa                                    >  2:3850/1‑139 (MQ=255)                                                                              cTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTaa                                    <  1:838136/139‑1 (MQ=255)                                                                                 ttGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCg                                                                                                      >  1:485621/1‑70 (MQ=255)                                                                                 ttGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCg                                                                                                      <  2:485621/70‑1 (MQ=255)                                                                                        aGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAg                          <  2:757385/139‑1 (MQ=255)                                                                                            ccATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCg                                                                                                      >  1:120985/1‑59 (MQ=255)                                                                                            ccATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCg                                                                                                      <  2:120985/59‑1 (MQ=255)                                                                                            ccATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAggtgg                      >  2:473262/1‑139 (MQ=255)                                                                                                    aGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGa                                        <  1:728846/113‑1 (MQ=255)                                                                                                    aGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCt              <  1:3850/139‑1 (MQ=255)                                                                                                    aGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCCCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGa                                        >  2:728846/1‑113 (MQ=255)                                                                                                      tAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTaa                                    >  1:890822/1‑115 (MQ=255)                                                                                                      tAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTaa                                    <  2:890822/115‑1 (MQ=255)                                                                                                           tACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTatgccatgcc                                                                                      <  1:191000/60‑1 (MQ=255)                                                                                                           tACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTatgccatgcc                                                                                      >  2:191000/1‑60 (MQ=255)                                                                                                           tACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTa       >  1:820077/1‑139 (MQ=255)                                                                                                             cGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACg     >  1:199251/1‑139 (MQ=255)                                                                                                              ggCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACgg    <  2:199251/139‑1 (MQ=255)                                                                                                              ggCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACgg    >  2:616560/1‑139 (MQ=255)                                                                                                                ctctctGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGcc  >  2:335524/1‑139 (MQ=255)                                                                                                                ctctctGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGcc  >  1:460505/1‑139 (MQ=255)                                                                                                                ctctctGCCCATAGCACCCCATCTTGCGCATTGCCCTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGcc  >  2:672972/1‑139 (MQ=255)                                                                                                                   tctGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAg                                                                                    >  1:881431/1‑54 (MQ=255)                                                                                                                   tctGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAg                                                                                    <  2:881431/54‑1 (MQ=255)                                                                                                                     tGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCAt                                                      >  2:476226/1‑82 (MQ=255)                                                                                                                     tGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCAt                                                      <  1:476226/82‑1 (MQ=255)                                                                                                                            aGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCgg                                                         >  2:671390/1‑72 (MQ=255)                                                                                                                            aGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCgg                                                         <  1:671390/72‑1 (MQ=255)                                                                                                                            aGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAg                          <  1:191437/103‑1 (MQ=255)                                                                                                                            aGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAg                          >  2:191437/1‑103 (MQ=255)                                                                                                                                |                                                                                                                           ATGACTTCCGCGCTTTTATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCC  >  NC_000913/115534‑115783

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