Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,128,358	T→C	E177E (GAA→GAG)	wcaF ←	putative acyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,128,358	0	T	C	100.0%	109.6 / NA	32	E177E (GAA→GAG)	wcaF	putative acyl transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/9);  total (23/9)

CTCGTAACCTTTTTCCAGCAGAAACTCTGCCAGGTAAGAACCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGG  >  NC_000913/2128231‑2128485                                                                                                                                |                                                                                                                                cTCGTAACCTTTTTCCAGCAGAAACTCTGCCAGGTAAGAACCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTcgcg                                                                                                                            >  1:312674/1‑133 (MQ=255) cTCGTAACCTTTTTCCAGCAGAAACTCTGCCAGGTAAGAACCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATGCCTCTGTATTTTAGAATTTATTCAGATTCAACGCGCTcgcg                                                                                                                            >  2:134094/1‑133 (MQ=255)           ttttCCAGCAGAAACTCTGCCAGGTAAGAACCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACt                                                                                                                  >  1:82214/1‑133 (MQ=255)                          tctGCCAGGTAAGAACCGTCTTGTCCGGTTACACCGGCGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCgg                                                                                                   <  1:322849/133‑1 (MQ=255)                             gCCAGGTAAGAACTGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCaa                                                                                                >  1:461223/1‑133 (MQ=255)                                   tAAGAACCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACcaca                                                                                          <  1:353978/133‑1 (MQ=255)                                        aCCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGc                                                                                     >  2:110010/1‑133 (MQ=255)                                        aCCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGc                                                                                     >  1:316653/1‑133 (MQ=255)                                            tCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGa                                                                                 >  1:324402/1‑133 (MQ=255)                                             cTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGaa                                                                                >  1:276727/1‑133 (MQ=255)                                               tatCCGGTTATACCGGTGATGAGTCCTACTTTAGACATGTACTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGc                                                                              <  2:239871/131‑1 (MQ=255)                                                 tCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGa                                                                            >  1:31192/1‑133 (MQ=255)                                                      ttACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTaa                                                                       <  1:275056/133‑1 (MQ=255)                                                           ccGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTaaaa                                                                     <  2:283584/130‑1 (MQ=255)                                                           ccGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTaaaa                                                                     >  1:283584/1‑130 (MQ=255)                                                               tgatgaGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGc                                                              >  1:372891/1‑133 (MQ=255)                                                                      agCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCACACCACATTCGCCGGGACCGATTTAACAACATTGCTTCGTGc                                                       >  2:231226/1‑133 (MQ=255)                                                                                         tattCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCg                                    >  1:121552/1‑133 (MQ=255)                                                                                                ctGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATtg                             >  1:297262/1‑133 (MQ=255)                                                                                                ctGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATtg                             >  1:374863/1‑133 (MQ=255)                                                                                                ctGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATtg                             >  2:25973/1‑133 (MQ=255)                                                                                                ctGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATtg                             >  1:315484/1‑133 (MQ=255)                                                                                                ctGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATtg                             >  2:303443/1‑133 (MQ=255)                                                                                                ctGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATtg                             >  2:341529/1‑133 (MQ=255)                                                                                                ctGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATtg                             >  2:353305/1‑133 (MQ=255)                                                                                                ctGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCGTCGTGCAACCACGACGGTGCCATCGCCGATtg                             >  2:159191/1‑133 (MQ=255)                                                                                                ctGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCATCGGGCACCCACGACGTTGCCATCGCCGATtg                             >  2:162121/1‑133 (MQ=255)                                                                                                 tGTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATtgt                            <  2:121552/133‑1 (MQ=255)                                                                                                  gTATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATtgtg                           <  2:31192/133‑1 (MQ=255)                                                                                                   tATTTTTGAATTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGa                          >  2:295451/1‑133 (MQ=255)                                                                                                       tttGAATTTACTCAGGTTCAACGCGCTCCCGTCTCACCACTTCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGcc                      <  2:155736/133‑1 (MQ=255)                                                                                                            aTTTATTCAGTTTCAACGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTgggg                 <  2:325840/133‑1 (MQ=255)                                                                                                                           aCGCGCTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCgg  >  2:320739/1‑133 (MQ=255)                                                                                                                                |                                                                                                                                CTCGTAACCTTTTTCCAGCAGAAACTCTGCCAGGTAAGAACCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGG  >  NC_000913/2128231‑2128485

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