Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	53,752	A→G	G317G (GGT→GGC)	surA ←	peptidyl‑prolyl cis‑trans isomerase SurA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	53,752	0	A	G	100.0%	74.4 / NA	21	G317G (GGT→GGC)	surA	peptidyl‑prolyl cis‑trans isomerase SurA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/11);  total (10/11)

CTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTCAGCAGAATATGGCGAGCATGAACTTCGGTCACCGAGATATTTTTGCT  >  NC_000913/53616‑53874                                                                                                                                         |                                                                                                                           cTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCa                                                                                                                         <  1:144192/140‑1 (MQ=255)    tGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCActc                                                                                                                      >  1:73298/1‑140 (MQ=255)    tGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCActc                                                                                                                      >  1:9867/1‑140 (MQ=255)           gCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGatat                                                                                                               <  2:392224/140‑1 (MQ=255)                           cGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTc                                                                                               <  2:325314/140‑1 (MQ=255)                                  ccGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTc                                                                                        <  1:111478/140‑1 (MQ=255)                                   cGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTca                                                                                       >  1:270251/1‑140 (MQ=255)                                        cGAAAATATCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACgg                                                                                  >  2:286986/1‑140 (MQ=255)                                          aaaaTATCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGc                                                                                <  2:73298/140‑1 (MQ=255)                                                tCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGc                                                                                >  2:523865/1‑134 (MQ=255)                                                tCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGc                                                                                <  1:523865/134‑1 (MQ=255)                                                  tGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCgt                                                                        <  1:169313/140‑1 (MQ=255)                                                              gCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCCGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCg                                                            >  1:334632/1‑140 (MQ=255)                                                                 gagaTCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCg                                                         <  2:140398/140‑1 (MQ=255)                                                                 gagaTCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCg                                                         <  1:34252/140‑1 (MQ=255)                                                                                          ccTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTCAGCAGAATATGGCGAGc                                >  1:175490/1‑140 (MQ=255)                                                                                            tGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTCAGCAGAATATGGCGAGCAt                              <  2:175490/140‑1 (MQ=255)                                                                                             ggATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTCAGCAGAATATGGCGAGCATg                             >  2:8508/1‑140 (MQ=255)                                                                                               aTCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTCAGCAGAATATGGCGAGCATGaa                           >  2:384929/1‑140 (MQ=255)                                                                                                      gagaACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTCAGCAGAATATGGCGAGCATGAACTTCGGt                    <  2:9867/140‑1 (MQ=255)                                                                                                                        gcagcaGAAGTCGTTTTGCCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTCAGCAGAATATGGCGAGCATGAACTTCGGTCACCGAGATATTTTTGCt  >  2:176942/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                           CTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCGCCGCCCTGGTTAGCAGAGCCTGGATCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTCAGCAGAATATGGCGAGCATGAACTTCGGTCACCGAGATATTTTTGCT  >  NC_000913/53616‑53874

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