Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	155,935	T→C	K89K (AAA→AAG)	yadV ←	putative periplasmic pilin chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	155,935	0	T	C	100.0%	212.7 / NA	59	K89K (AAA→AAG)	yadV	putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/28);  total (31/28)

AGCAGGCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACG  >  NC_000913/155801‑156066                                                                                                                                       |                                                                                                                                    agcagGCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTgg                                                                                                                                 >  1:383109/1‑139 (MQ=255)                gcgACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCAt                                                                                                                              <  1:599762/127‑1 (MQ=255)                gcgACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCAt                                                                                                                              >  2:599762/1‑127 (MQ=255)                  gACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACt                                                                                                                <  1:609188/139‑1 (MQ=255)                                  ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACt                                                                                                                <  1:121440/123‑1 (MQ=255)                                  ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACt                                                                                                                >  2:121440/1‑123 (MQ=255)                                  ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaa                                                                                                >  1:480357/1‑139 (MQ=255)                                      ttGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                            <  2:674967/139‑1 (MQ=255)                                              aCTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTaa                                                                                    >  2:453240/1‑139 (MQ=255)                                                    aGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGc                                                                              <  2:480357/139‑1 (MQ=255)                                                       aCGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTgg                                                                                                              <  2:316168/104‑1 (MQ=255)                                                       aCGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTgg                                                                                                              >  1:316168/1‑104 (MQ=255)                                                          ttAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                      >  1:776437/1‑109 (MQ=255)                                                          ttAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                      <  2:776437/109‑1 (MQ=255)                                                           tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                            >  2:710353/1‑118 (MQ=255)                                                           tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                            <  1:710353/118‑1 (MQ=255)                                                           tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTc                                                                       >  2:455140/1‑139 (MQ=255)                                                                 aGAACGCGTTCTCTCTGCCTTGAGGCGGTGAGGTGCTTGCTGTGAACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTGGCATTAAAAGGGACTGTAATTCTGCCAGGCTCAGCGtt                                                                 <  2:159506/139‑1 (MQ=255)                                                                   aaCACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                            >  2:218359/1‑110 (MQ=255)                                                                   aaCACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                            <  1:218359/110‑1 (MQ=255)                                                                    acacGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTc                                                              <  1:453240/139‑1 (MQ=255)                                                                         ctctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                      <  2:809590/94‑1 (MQ=255)                                                                         ctctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                      >  1:809590/1‑94 (MQ=255)                                                                         ctctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGc                                                         <  1:455140/139‑1 (MQ=255)                                                                          tctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGcc                                                        <  2:383109/139‑1 (MQ=255)                                                                                          cAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGTCTGTAATACTGCCAGGCTCAGCGTTGTCAt                                                            <  1:609075/119‑1 (MQ=255)                                                                                          cAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGTCTGTAATACTGCCAGGCTCAGCGTTGTCAt                                                            >  2:609075/1‑119 (MQ=255)                                                                                            gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAcg                                                                                                                       <  1:28769/58‑1 (MQ=255)                                                                                            gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAcg                                                                                                                       >  2:28769/1‑58 (MQ=255)                                                                                            gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGc                                                                              <  2:831975/99‑1 (MQ=255)                                                                                            gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGc                                                                              >  1:831975/1‑99 (MQ=255)                                                                                            gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGcc                                                        >  1:468950/1‑121 (MQ=255)                                                                                            gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGcc                                                        <  2:468950/121‑1 (MQ=255)                                                                                            gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAcc                                             <  1:771580/132‑1 (MQ=255)                                                                                            gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAcc                                             >  2:771580/1‑132 (MQ=255)                                                                                             tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTaa                                                                                    >  2:423389/1‑92 (MQ=255)                                                                                             tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTaa                                                                                    <  1:423389/92‑1 (MQ=255)                                                                                              gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAcgc                                                                                                                      <  2:324226/57‑1 (MQ=255)                                                                                              gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAcgc                                                                                                                      >  1:324226/1‑57 (MQ=255)                                                                                              gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                            >  1:550070/1‑83 (MQ=255)                                                                                              gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                            <  2:550070/83‑1 (MQ=255)                                                                                                         tgtgTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                      <  2:396434/62‑1 (MQ=255)                                                                                                         tgtgTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                      >  1:396434/1‑62 (MQ=255)                                                                                                           tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGa                                    >  1:744547/1‑126 (MQ=255)                                                                                                           tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGa                                    <  2:744547/126‑1 (MQ=255)                                                                                                                 ttAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCa                                                                                                    <  1:734927/56‑1 (MQ=255)                                                                                                                 ttAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCa                                                                                                    >  2:734927/1‑56 (MQ=255)                                                                                                                  tAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCt                >  2:473712/1‑139 (MQ=255)                                                                                                                       tGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGc                                                                              <  1:556522/72‑1 (MQ=255)                                                                                                                       tGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGc                                                                              >  2:556522/1‑72 (MQ=255)                                                                                                                       tGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAAGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTAtt           >  2:718949/1‑139 (MQ=255)                                                                                                                       tGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAAGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTAtt           >  2:40296/1‑139 (MQ=255)                                                                                                                          ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTc        >  2:236800/1‑139 (MQ=255)                                                                                                                             tttGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAg     >  2:809728/1‑139 (MQ=255)                                                                                                                              ttGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCa                              <  2:496075/113‑1 (MQ=255)                                                                                                                              ttGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCa                              >  1:496075/1‑113 (MQ=255)                                                                                                                                gCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGacg  >  1:387519/1‑139 (MQ=255)                                                                                                                                   cACGCTTGGCATCAATACGTGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGtt                                                                 <  1:380692/73‑1 (MQ=255)                                                                                                                                   cACGCTTGGCATCAATACGTGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGtt                                                                 >  2:380692/1‑73 (MQ=255)                                                                                                                                       |                                                                                                                                    AGCAGGCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACG  >  NC_000913/155801‑156066

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