Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,125,097	C→T	100%	E572K (GAG→AAG)	priA ←	primosome factor N'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,125,097	0	C	T	100.0%	90.8 / NA	34	E572K (GAG→AAG)	priA	primosome factor N'
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/19);  total (15/19)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGAT  >  NC_000913/4124957‑4125225                                                                                                                                              |                                                                                                                                 ggCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTg                                                                                                                                <  1:26337/143‑1 (MQ=255) ggCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTg                                                                                                                                >  1:359763/1‑143 (MQ=255)  gCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGc                                                                                                                               <  2:273604/143‑1 (MQ=255)   cAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGGAAAGCTGAGCGAAACGCTTTGcc                                                                                                                              >  1:131536/1‑143 (MQ=255)    aaaGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCATTGCCg                                                                                                                             >  1:386918/1‑143 (MQ=255)      aGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAg                                                                                                                           <  1:590773/143‑1 (MQ=255)        gCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTGTACTCGCACGCCCGGCACGACCGGCGACCGG‑GGTGGAAAGCTGAGCCAAACGCTTTGcccggcg                                                                                                                         >  1:212962/1‑138 (MQ=255)            cgtAGCCTTTATAGAGCAACGTTAGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGaaaa                                                                                                                     >  2:579608/1‑143 (MQ=255)                gCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCgg                                                                                                                 <  2:245796/143‑1 (MQ=255)                              aCGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAgcgcg                                                                                                   >  1:15163/1‑143 (MQ=255)                                     cAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTGGCCGAGCGAAAATCGGCAGAAAACAGCGCGCcgtcca                                                                                            <  1:8711/143‑1 (MQ=255)                                          gagGATGTTCCGTATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCATAAAACAGCGCGCcgtccacgtcc                                                                                       <  2:654456/143‑1 (MQ=255)                                           agGATGTTCCGGATGGTGCGGTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCcgtccacgtcca                                                                                      <  2:172023/143‑1 (MQ=255)                                                ggtCCGGATGCTGCGCTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAAt                                                                                 <  1:661972/141‑1 (MQ=255)                                                gTTCCTGATTGTGCGTTCGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGACACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAAt                                                                                 <  2:377391/143‑1 (MQ=255)                                                gTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAGCGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAAt                                                                                 >  1:507445/1‑143 (MQ=255)                                                gTTCCGGATGGTGCGTTTGCAGCACCACCTCGACTTCCTTCCCCGCACGCCCGGCACGAATGGC‑ACATGCGCTCTAAAACTGTGAGCATCTCTTTCCCTTTTTAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAAt                                                                                 <  1:306187/143‑1 (MQ=255)                                                            gCGTTTGCAGCACCACTTCGCCCTGTTGACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTc                                                                     <  1:510956/143‑1 (MQ=255)                                                             cGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGGca                                                                    >  1:465326/1‑143 (MQ=255)                                                              gtctGCAGCACCACTCCGCCCTGTTTGCCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGGCCACGTCCAGTAATGCAACCAGCGTcac                                                                   <  1:183737/140‑1 (MQ=255)                                                                    agcaCCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCGCATCCgg                                                             >  1:528741/1‑143 (MQ=255)                                                                            ttCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGCAATGGt                                                     >  1:541006/1‑143 (MQ=255)                                                                            ttCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCACCGTCCACGTCCAGTAATGCAACCAGCGTCCCGTCCGGGCAATGGt                                                     >  1:195945/1‑143 (MQ=255)                                                                              cGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGCGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGa                                                   <  1:268451/143‑1 (MQ=255)                                                                                cccTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCAAGTACAGAAACCAACGTCACATCCGGGAAATGGTGAcc                                                 >  1:512106/1‑143 (MQ=255)                                                                                         cccGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAg                                        >  1:457055/1‑143 (MQ=255)                                                                                            gCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGccc                                     >  2:167844/1‑141 (MQ=255)                                                                                                   cGGCGCGCCCGGCGCCCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTa                              <  2:549818/143‑1 (MQ=255)                                                                                                   cGGCACGACCGGCGACCGG‑GGGGCAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGGCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTa                              <  2:550750/143‑1 (MQ=255)                                                                                                      cACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCa                           <  1:407089/143‑1 (MQ=255)                                                                                                            cgcctaCCTG‑GGTGTAAAGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCaaa                     <  1:582571/138‑1 (MQ=255)                                                                                                                 aCCTG‑GGTGTAAAGCTGAGCTAAGCGCTTTGCCGAGCGAGAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCcg                <  2:177923/143‑1 (MQ=255)                                                                                                                            aaaGCTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACCTCCGGGAAGTGGTGACCTTTCGCCACCATTTGTGTACCACTCAAAATCCGCGggcggccg      >  1:426047/1‑143 (MQ=255)                                                                                                                                cTGAGCGAAACGCTTTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGAt  <  2:15163/143‑1 (MQ=255)                                                                                                                                              |                                                                                                                                 GGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTG‑GGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGAT  >  NC_000913/4124957‑4125225

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