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%	199.8 / NA	56	K89K (AAA→AAG)	yadV	putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/23);  total (33/23)

TGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACG  >  NC_000913/155810‑156066                                                                                                                                |                                                                                                                                    tGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAc                                                                                                                        >  1:672079/1‑139 (MQ=255) tGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAc                                                                                                                        <  1:569207/139‑1 (MQ=255)  gATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAcg                                                                                                                       >  2:954061/1‑139 (MQ=255)      cgcgACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGAt                                                                                                                   >  1:103674/1‑139 (MQ=255)                         ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCa                                                                                                                            <  1:832472/111‑1 (MQ=255)                         ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCa                                                                                                                            >  2:832472/1‑111 (MQ=255)                         ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTa                                                                                                       <  2:801649/132‑1 (MQ=255)                         ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTa                                                                                                       >  1:801649/1‑132 (MQ=255)                         ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaa                                                                                                >  1:269115/1‑139 (MQ=255)                           ttttGGTGGAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaaa                                                                                              <  2:672079/139‑1 (MQ=255)                              tggtggAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAggg                                                                                           >  2:442129/1‑139 (MQ=255)                              tggtggAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTGGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAggg                                                                                           <  1:954061/139‑1 (MQ=255)                               ggtggAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                            <  2:1011037/139‑1 (MQ=255)                                gtggAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGAc                                                                                         <  2:791208/139‑1 (MQ=255)                                              aCGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAg                                                                           >  1:301816/1‑139 (MQ=255)                                                  tAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAc                                                                                                                        >  1:442801/1‑90 (MQ=255)                                                  tAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAc                                                                                                                        <  2:442801/90‑1 (MQ=255)                                                  tAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTc                                                                       >  1:985491/1‑139 (MQ=255)                                                          aaCACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGt                                                               >  1:205171/1‑139 (MQ=255)                                                             acG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCAt                                                            <  2:830013/139‑1 (MQ=255)                                                              cG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATc                                                           >  1:898640/1‑139 (MQ=255)                                                               g‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCg                                                          >  2:830101/1‑139 (MQ=255)                                                               g‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCg                                                          >  2:986189/1‑139 (MQ=255)                                                                   tctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGcc                                                        <  2:114445/139‑1 (MQ=255)                                                                        ctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTAt                                                   >  2:336296/1‑139 (MQ=255)                                                                             tttAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAc                                              <  2:205171/139‑1 (MQ=255)                                                                               tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCAt                                                                                                                              <  2:204930/57‑1 (MQ=255)                                                                               tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCAt                                                                                                                              >  1:204930/1‑57 (MQ=255)                                                                                aGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCCCCAGTATCTAACCaa                                           >  1:249762/1‑139 (MQ=255)                                                                                 ggCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAAc                                          >  2:110940/1‑139 (MQ=255)                                                                                     gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTg                                      >  1:628081/1‑139 (MQ=255)                                                                                     gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTg                                      >  1:632994/1‑139 (MQ=255)                                                                                      tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGt                                                                                                        >  1:64418/1‑72 (MQ=255)                                                                                      tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGt                                                                                                        <  2:64418/72‑1 (MQ=255)                                                                                      tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTg                                                                               >  1:173778/1‑97 (MQ=255)                                                                                      tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTg                                                                               <  2:173778/97‑1 (MQ=255)                                                                                       gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGa                                    >  2:958660/1‑139 (MQ=255)                                                                                       gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGa                                    >  2:1009931/1‑139 (MQ=255)                                                                                        aGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCGTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGAc                                   <  2:632994/139‑1 (MQ=255)                                                                                         ggTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACa                                  <  1:442129/139‑1 (MQ=255)                                                                                             ctggctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCa                              <  1:211946/139‑1 (MQ=255)                                                                                              tggctgTGTACATTAATTTGATTGTTGGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCaa                             <  2:285251/139‑1 (MQ=255)                                                                                                gctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGc                                                                                                                                >  1:80820/1‑38 (MQ=255)                                                                                                gctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGc                                                                                                                                <  2:80820/38‑1 (MQ=255)                                                                                                    tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGtt                       >  1:1021787/1‑139 (MQ=255)                                                                                                              tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTggcgg                                                                                                           >  2:687080/1‑45 (MQ=255)                                                                                                              tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTggcgg                                                                                                           <  1:687080/45‑1 (MQ=255)                                                                                                              tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTa             <  2:628081/139‑1 (MQ=255)                                                                                                                tGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTAtt           >  2:607823/1‑139 (MQ=255)                                                                                                                tGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTAtt           >  2:157119/1‑139 (MQ=255)                                                                                                                      tttGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAg     >  2:63948/1‑139 (MQ=255)                                                                                                                       ttGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGa    >  1:715179/1‑139 (MQ=255)                                                                                                                        tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                      >  2:661404/1‑40 (MQ=255)                                                                                                                        tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                      <  1:661404/40‑1 (MQ=255)                                                                                                                        tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGac   >  2:1030735/1‑139 (MQ=255)                                                                                                                         gCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGacg  <  1:157119/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                    TGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACG‑‑CTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACG  >  NC_000913/155810‑156066

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