Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,288,276	A→G	N63D (AAT→GAT)	yraI →	putative periplasmic pilin chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,288,276	0	A	G	100.0%	135.0 / NA	36	N63D (AAT→GAT)	yraI	putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/20);  total (16/20)

GCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATAATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAAGAGTCGCTGTTCTAC  >  NC_000913/3288133‑3288410                                                                                                                                                |                                                                                                                                       gCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGa                                                                                                                                 <  2:241419/151‑1 (MQ=255)  cTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGatatcgct                                                                                                                                >  1:48327/1‑147 (MQ=255)       gTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTaa                                                                                                                           <  2:33379/151‑1 (MQ=255)             gCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCg                                                                                                                     <  2:306261/151‑1 (MQ=255)              ccAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGt                                                                                                                    >  1:250255/1‑151 (MQ=255)               cAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTg                                                                                                                   <  1:81523/151‑1 (MQ=255)                        gCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATtta                                                                                                          <  1:240223/151‑1 (MQ=255)                                aaTCGTTTTACAGCGAACGCGAGTGATCTATGATGCCTGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCTAACATGGGTAGATGATATAGATGGTAAAAGCCGTGCACCATTTATTATATcc                                                                                                  <  1:766417/151‑1 (MQ=255)                                aaTCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAAcc                                                                                                  <  2:371857/151‑1 (MQ=255)                                   cGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAAcccc                                                                                                <  1:703820/150‑1 (MQ=255)                                   cGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAAcccc                                                                                                <  2:263007/150‑1 (MQ=255)                                    gTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCAc                                                                                              <  2:233983/151‑1 (MQ=255)                                      tttACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCAcc                                                                                             <  1:90413/150‑1 (MQ=255)                                      tttACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCg                                                                                            >  2:308591/1‑151 (MQ=255)                                         aCAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTa                                                                                         <  2:237279/151‑1 (MQ=255)                                          cAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTAt                                                                                        <  1:388141/151‑1 (MQ=255)                                              gAACGCGAGGGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCg                                                                                    >  2:384029/1‑151 (MQ=255)                                                     aGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAg                                                                             <  2:594069/151‑1 (MQ=255)                                                       tGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGc                                                                           >  2:625682/1‑151 (MQ=255)                                                        gATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCt                                                                          >  2:770103/1‑151 (MQ=255)                                                            tATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCg                                                                      >  2:375959/1‑151 (MQ=255)                                                                atgCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATg                                                                   >  1:742406/1‑150 (MQ=255)                                                                        cGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACtcatca                                                           <  2:406252/150‑1 (MQ=255)                                                                           aaaGAGGCTGCGTTACCTGTCGCACACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCAc                                                          >  2:492462/1‑148 (MQ=255)                                                                             agagGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCg                                                      >  2:575154/1‑150 (MQ=255)                                                                                   tGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAattatt                                               >  1:751976/1‑151 (MQ=255)                                                                                         aCCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACa                                         <  2:137889/151‑1 (MQ=255)                                                                                             gTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTg                                     <  1:201933/151‑1 (MQ=255)                                                                                                               gAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAAt                    >  1:205150/1‑150 (MQ=255)                                                                                                                 aaCGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATaa                  >  2:426843/1‑150 (MQ=255)                                                                                                                  aCGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAag                <  2:631603/151‑1 (MQ=255)                                                                                                                    gCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAagag              <  2:411082/151‑1 (MQ=255)                                                                                                                      cTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAAGAGTc            >  1:42879/1‑151 (MQ=255)                                                                                                                       ttATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAAGAGTCg           >  1:581244/1‑151 (MQ=255)                                                                                                                         aTTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAAGAGTCGCt         <  1:318294/151‑1 (MQ=255)                                                                                                                                 cAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAAGAGTCGCTGTTCTAc  >  1:158828/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                       GCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATAATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAAGAGTCGCTGTTCTAC  >  NC_000913/3288133‑3288410

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