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%	132.4 / NA	35	N63D (AAT→GAT)	yraI	putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/20);  total (15/20)

TGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATAATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAAGAGTCGCTGTTCTACATTAATGTTCGT  >  NC_000913/3288129‑3288422                                                                                                                                                    |                                                                                                                                                   tgtAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAgatgat                                                                                                                                                  >  2:409548/1‑150 (MQ=255)           gTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGt                                                                                                                                         <  1:70534/149‑1 (MQ=255)            tATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTaa                                                                                                                                       >  1:545775/1‑150 (MQ=255)                  ccAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGt                                                                                                                                <  1:105762/151‑1 (MQ=255)                    aGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTg                                                                                                                               <  2:715262/150‑1 (MQ=255)                        gcTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGcccc                                                                                                                           >  1:586605/1‑150 (MQ=255)                                 aggaATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTtattat                                                                                                                  >  1:479104/1‑150 (MQ=255)                                      tCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAAcccc                                                                                                            >  1:573127/1‑151 (MQ=255)                                        gTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCAc                                                                                                          >  1:233270/1‑151 (MQ=255)                                             aCAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCt                                                                                                      <  2:880356/150‑1 (MQ=255)                                             aCAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCt                                                                                                      <  1:470756/150‑1 (MQ=255)                                                gCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTAttt                                                                                                  >  2:576671/1‑151 (MQ=255)                                                    aCGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTc                                                                                              <  1:220162/151‑1 (MQ=255)                                                                 atgatgCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTgg                                                                                    <  1:611065/148‑1 (MQ=255)                                                                 atgatgCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTg                                                                                     <  1:73476/147‑1 (MQ=255)                                                                 atgatgCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGc                                                                                   <  1:447078/149‑1 (MQ=255)                                                                 atgatgCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGc                                                                                   <  1:347320/149‑1 (MQ=255)                                                                 atgatgCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCg                                                                                  <  1:377545/150‑1 (MQ=255)                                                                   gatgCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGAt                                                                                <  1:631038/150‑1 (MQ=255)                                                                            cGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCAc                                                                      <  2:527561/151‑1 (MQ=255)                                                                                   agGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAat                                                               <  1:503759/151‑1 (MQ=255)                                                                                   agGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAat                                                               >  1:62907/1‑151 (MQ=255)                                                                                        gCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAattatt                                                           >  2:326564/1‑150 (MQ=255)                                                                                                      aaacaaaGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAAc                                            <  1:140731/151‑1 (MQ=255)                                                                                                          aaaGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAAcctgc                                        <  1:579735/151‑1 (MQ=255)                                                                                                              gcgcAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAAcctgcctg                                     >  1:603798/1‑150 (MQ=255)                                                                                                               cgcAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGaa                                   <  2:612110/151‑1 (MQ=255)                                                                                                                gcAGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGaa                                   >  2:395417/1‑150 (MQ=255)                                                                                                                  aGAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAAt                                >  2:219285/1‑151 (MQ=255)                                                                                                                   gAAACGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAAt                                >  1:610033/1‑150 (MQ=255)                                                                                                                      aCGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTTATAACCTGCCTGAAAATaa                              <  1:331002/149‑1 (MQ=255)                                                                                                                      aCGCCTTATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATaa                              >  1:586383/1‑149 (MQ=255)                                                                                                                            tATTTACTGCAATCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAAGAGTc                        <  1:298301/149‑1 (MQ=255)                                                                                                                                       aTCATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAAGAGTCGCTGTTCTACAtt           <  2:479057/151‑1 (MQ=255)                                                                                                                                         cATGGGTAGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAAGAGTCGCTGTTCTACATTa          >  1:72620/1‑150 (MQ=255)                                                                                                                                                aGATGATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAAGAGTCGCTGTTCTACATTAATGTTCGt  >  1:93157/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACGCGAGTGATCTATGATGCCAGCCGCAAAGAGGCTGCGTTACCTGTCGCAAACAAAGGCGCAGAAACGCCTTATTTACTGCAATCATGGGTAGATAATATAGATGGTAAAAGCCGTGCCCCATTTATTATAACCCCACCGCTATTTCGTCTTGAGGCTGGCGATGACTCATCACTGCGAATTATTAAAACAGCTGATAACCTGCCTGAAAATAAAGAGTCGCTGTTCTACATTAATGTTCGT  >  NC_000913/3288129‑3288422

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