Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	657,377	T→C	V29A (GTG→GCG)	cspE →	constitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	657,377	0	T	C	100.0%	123.7 / NA	34	V29A (GTG→GCG)	cspE	constitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/20);  total (14/20)

GTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGA  >  NC_000913/657238‑657516                                                                                                                                            |                                                                                                                                            gtACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACt                                                                                                                                  >  1:46485/1‑151 (MQ=255)         aaCCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCa                                                                                                                          <  2:454048/151‑1 (MQ=255)            cGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATc                                                                                                                       <  2:315748/151‑1 (MQ=255)                       ttGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATg                                                                                                             <  2:109359/150‑1 (MQ=255)                           gtCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                                                                                                         <  2:253185/150‑1 (MQ=255)                              tATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa                                                                                                      >  2:550949/1‑150 (MQ=255)                               aTTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa                                                                                                      <  1:219464/149‑1 (MQ=255)                                  tttCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAct                                                                                                  <  1:193587/150‑1 (MQ=255)                                      aTGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGc                                                                                             <  2:174016/151‑1 (MQ=255)                                       tGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCt                                                                                            >  1:308422/1‑151 (MQ=255)                                       tGTAAAGGTAATTTTGATGGCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCt                                                                                            >  2:238935/1‑151 (MQ=255)                                         tAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg                                                                                           >  1:434127/1‑150 (MQ=255)                                          aaaGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa                                                                                          <  2:23353/150‑1 (MQ=255)                                           aaGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa                                                                                          >  2:174725/1‑149 (MQ=255)                                               tAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTc                                                                                     >  2:601382/1‑150 (MQ=255)                                                aaTTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTc                                                                                     >  1:296402/1‑149 (MQ=255)                                                aaTTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAg                                                                                   <  2:629581/151‑1 (MQ=255)                                                   tttGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg                                                                                 <  1:173276/150‑1 (MQ=255)                                                      gATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGt                                                                              <  1:561987/150‑1 (MQ=255)                                                      gATGTCTAAGATTAAAGGCAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGt                                                                              >  2:577069/1‑150 (MQ=255)                                                             aaGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCg                                                                      >  1:375684/1‑151 (MQ=255)                                                                  tAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATc                                                                 >  2:586355/1‑151 (MQ=255)                                                                   aaaGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATc                                                                 >  2:472389/1‑150 (MQ=255)                                                                    aaGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCAc                                                               <  2:543587/151‑1 (MQ=255)                                                                    aaGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCAc                                                               >  2:563588/1‑151 (MQ=255)                                                                    aaGGCAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCAc                                                               >  1:392008/1‑151 (MQ=255)                                                                       gTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTaa                                                            <  2:402605/151‑1 (MQ=255)                                                                                                  aaGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCaaa                                 <  2:68486/151‑1 (MQ=255)                                                                                                          ggTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCg                         <  1:656450/151‑1 (MQ=255)                                                                                                                      ccGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATa             <  2:398307/151‑1 (MQ=255)                                                                                                                         gAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTATTCGCTCTGTAAGATACGt          <  2:682716/151‑1 (MQ=255)                                                                                                                          aaGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTc         <  1:398863/151‑1 (MQ=255)                                                                                                                               ggCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCaa    <  1:168999/151‑1 (MQ=255)                                                                                                                                  agcaAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGa  <  1:457758/150‑1 (MQ=255)                                                                                                                                            |                                                                                                                                            GTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGA  >  NC_000913/657238‑657516

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