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%	182.8 / NA	51	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 (27/24);  total (27/24)

TATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCA  >  NC_000913/115550‑115788                                                                                                                |                                                                                                                                tATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGt                                                                                                      <  1:394334/139‑1 (MQ=255)            cACTGCGCAGAATTGATTGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTatgccatgcc                                                                                           <  1:595803/139‑1 (MQ=255)                 cgcAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAgcgc                                                                                      >  1:1829601/1‑139 (MQ=255)                    aGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCag                                                                                   >  2:74418/1‑139 (MQ=255)                      aaTTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCagag                                                                                 >  1:1197857/1‑139 (MQ=255)                      aaTTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCagag                                                                                 >  2:591854/1‑139 (MQ=255)                           aTGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTatgccatgc                                                                                            >  2:1292499/1‑123 (MQ=255)                           aTGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTatgccatgc                                                                                            <  1:1292499/123‑1 (MQ=255)                                gATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTc                                                                       >  2:1421014/1‑139 (MQ=255)                                 aTTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCg                                                                      >  1:1509470/1‑139 (MQ=255)                                     gCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCAt                                                                  <  2:1197857/139‑1 (MQ=255)                                     gCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCAt                                                                  <  1:74418/139‑1 (MQ=255)                                          aGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGc                                                             <  2:1829601/139‑1 (MQ=255)                                                      gTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTAtgccatgcca                                                                                          >  1:1863983/1‑98 (MQ=255)                                                      gTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTAtgccatgcca                                                                                          <  2:1863983/98‑1 (MQ=255)                                                      gTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCagag                                                                                 <  1:236475/107‑1 (MQ=255)                                                      gTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCagag                                                                                 >  2:236475/1‑107 (MQ=255)                                                            aTCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGc                                                                >  2:598719/1‑118 (MQ=255)                                                            aTCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGc                                                                <  1:598719/118‑1 (MQ=255)                                                            aTCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGAtt                                           <  2:89150/139‑1 (MQ=255)                                                                gTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAAct                                       >  2:716292/1‑139 (MQ=255)                                                                       tAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCaaa                                >  1:951639/1‑139 (MQ=255)                                                                        aGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAgcgc                                                                                      <  2:771504/84‑1 (MQ=255)                                                                        aGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAgcgc                                                                                      >  1:771504/1‑84 (MQ=255)                                                                        aGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAg                               <  1:359500/139‑1 (MQ=255)                                                                         gTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCa                                                                                    >  2:1650328/1‑85 (MQ=255)                                                                         gTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCa                                                                                    <  1:1650328/85‑1 (MQ=255)                                                                           gCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAggtg                            >  2:760031/1‑139 (MQ=255)                                                                            ccATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAggtgg                           >  2:600628/1‑139 (MQ=255)                                                                                    aGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCg                                                                                                           >  2:962986/1‑51 (MQ=255)                                                                                    aGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCg                                                                                                           <  1:962986/51‑1 (MQ=255)                                                                                    aGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGc                                                             >  1:1175388/1‑97 (MQ=255)                                                                                    aGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGc                                                             <  2:1175388/97‑1 (MQ=255)                                                                                    aGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCt                   <  1:701079/139‑1 (MQ=255)                                                                                         aGTACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGctc              <  1:591854/139‑1 (MQ=255)                                                                                           tACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCAt                                                                  >  1:1313197/1‑85 (MQ=255)                                                                                           tACGGCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCAt                                                                  <  2:1313197/85‑1 (MQ=255)                                                                                              ggCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCAt                                                                  >  1:664280/1‑82 (MQ=255)                                                                                              ggCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCAt                                                                  <  2:664280/82‑1 (MQ=255)                                                                                              ggCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTAcg                                                 >  1:1124045/1‑99 (MQ=255)                                                                                              ggCTCTCTGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTAcg                                                 <  2:1124045/99‑1 (MQ=255)                                                                                                ctctctGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCtt                  >  1:1607327/1‑128 (MQ=255)                                                                                                ctctctGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCtt                  <  2:1607327/128‑1 (MQ=255)                                                                                                ctctctGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGcc       >  1:1752411/1‑139 (MQ=255)                                                                                                   tctGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGt                                                                                                      <  2:1854125/41‑1 (MQ=255)                                                                                                   tctGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGt                                                                                                      >  1:1854125/1‑41 (MQ=255)                                                                                                     tGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTatgcc                                                                                                >  2:421758/1‑45 (MQ=255)                                                                                                     tGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTatgcc                                                                                                <  1:421758/45‑1 (MQ=255)                                                                                                     tGCCCATAGCACCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCCGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCa  >  2:1821043/1‑139 (MQ=255)                                                                                                            aGCACCCCATCTTGCGCATTGCCGTCGCCGGTTAtgccatgcca                                                                                          >  1:1716347/1‑44 (MQ=255)                                                                                                            aGCACCCCATCTTGCGCATTGCCGTCGCCGGTTAtgccatgcca                                                                                          <  2:1716347/44‑1 (MQ=255)                                                                                                                |                                                                                                                                TATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCA  >  NC_000913/115550‑115788

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