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

GGCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGAC  >  NC_000913/155805‑156065                                                                                                                                   |                                                                                                                                   ggCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAATACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATc                                                                                                                            >  1:182555/1‑139 (MQ=255)      tGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAc                                                                                                                       >  2:682189/1‑139 (MQ=255)                        gCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                     <  2:962077/139‑1 (MQ=255)                                  ttGGTGGAACTTCCAATACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                           <  2:182555/139‑1 (MQ=255)                                         aaCTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTa                                                                                    >  1:163179/1‑139 (MQ=255)                                         aaCTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTa                                                                                    >  1:508297/1‑139 (MQ=255)                                                                    gCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTgg                                                                                                             <  1:286790/87‑1 (MQ=255)                                                                    gCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTgg                                                                                                             >  2:286790/1‑87 (MQ=255)                                                                      tctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTgg                                                                                                             >  1:400928/1‑85 (MQ=255)                                                                      tctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTgg                                                                                                             <  2:400928/85‑1 (MQ=255)                                                                      tctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGcc                                                       <  2:163179/139‑1 (MQ=255)                                                                      tctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGcc                                                       <  2:691866/139‑1 (MQ=255)                                                                         ctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                     >  2:434158/1‑90 (MQ=255)                                                                         ctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGt                                                    >  1:729761/1‑139 (MQ=255)                                                                         ctctGTCTTTAGGCAGTGAGGGGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                     <  1:434158/90‑1 (MQ=255)                                                                          tctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTa                                                   >  2:427458/1‑139 (MQ=255)                                                                                  tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTc                                                             <  1:833332/121‑1 (MQ=255)                                                                                  tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTc                                                             >  2:833332/1‑121 (MQ=255)                                                                                  tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCa                                           >  1:806223/1‑139 (MQ=255)                                                                                        gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGa                                                                                         <  1:880392/87‑1 (MQ=255)                                                                                        gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGa                                                                                         >  2:880392/1‑87 (MQ=255)                                                                                        gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGa                                                                                         <  1:856999/87‑1 (MQ=255)                                                                                        gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGa                                                                                         >  2:856999/1‑87 (MQ=255)                                                                                        gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCa                                                            >  1:1277/1‑116 (MQ=255)                                                                                        gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCa                                                            <  2:1277/116‑1 (MQ=255)                                                                                         tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAcc                                            <  2:296578/131‑1 (MQ=255)                                                                                         tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAcc                                            >  1:296578/1‑131 (MQ=255)                                                                                          gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGt                                                                                                       >  1:90825/1‑71 (MQ=255)                                                                                          gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGt                                                                                                       <  2:90825/71‑1 (MQ=255)                                                                                                       tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                     <  1:827663/60‑1 (MQ=255)                                                                                                       tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                     >  2:827663/1‑60 (MQ=255)                                                                                                       tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                     >  2:549470/1‑60 (MQ=255)                                                                                                       tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                     <  1:549470/60‑1 (MQ=255)                                                                                                       tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGcc                                                                            >  1:806431/1‑85 (MQ=255)                                                                                                       tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGcc                                                                            <  2:806431/85‑1 (MQ=255)                                                                                                       tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGcc                                                                            >  2:625264/1‑85 (MQ=255)                                                                                                       tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGcc                                                                            <  1:625264/85‑1 (MQ=255)                                                                                                        gtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAcg                                                                                                                      <  1:907636/42‑1 (MQ=255)                                                                                                        gtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAcg                                                                                                                      >  2:907636/1‑42 (MQ=255)                                                                                                        gtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGa                                                                                         >  1:121143/1‑71 (MQ=255)                                                                                                        gtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGa                                                                                         <  2:121143/71‑1 (MQ=255)                                                                                                        gtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCTGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGa                                   >  2:89808/1‑125 (MQ=255)                                                                                                        gtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCTGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGa                                   <  1:89808/125‑1 (MQ=255)                                                                                                                 tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTa            <  1:75232/139‑1 (MQ=255)                                                                                                                 tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTa            >  1:226240/1‑139 (MQ=255)                                                                                                                 tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTa            >  1:129585/1‑139 (MQ=255)                                                                                                                     aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTgg                                                                                                             <  2:936230/38‑1 (MQ=255)                                                                                                                     aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTgg                                                                                                             >  1:936230/1‑38 (MQ=255)                                                                                                                          ttGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGa   <  1:682189/139‑1 (MQ=255)                                                                                                                           tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGac  >  2:739557/1‑139 (MQ=255)                                                                                                                            gCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                     <  2:732667/39‑1 (MQ=38)                                                                                                                            gCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                     >  1:732667/1‑39 (MQ=38)                                                                                                                                   |                                                                                                                                   GGCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGAC  >  NC_000913/155805‑156065

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: ‑