Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,926,857	Δ1 bp	coding (147/459 nt)	asnC ←	DNA‑binding transcriptional dual regulator AsnC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,926,854	0	C	.	100.0%	131.5 / NA	32	coding (150/459 nt)	asnC	DNA‑binding transcriptional dual regulator AsnC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (17/15);  total (17/15)

TTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACG  >  NC_000913/3926722‑3926973                                                                                                                                     |                                                                                                                        ttCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATc                                                                                                                         >  2:125902/1‑133 (MQ=255)         ggCTTTCCAGCTTTGCCAGCGCGGAAGGTTAGACTTTGGCGTTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGc                                                                                                                <  2:265696/133‑1 (MQ=255)         ggCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATTAATACGCGCCCGGGGAATGAT‑GCCGCCTGc                                                                                                                >  1:277242/1‑133 (MQ=255)            tttCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTc                                                                                                             <  1:349925/133‑1 (MQ=255)                        ccAGCGCGCAAGGGTCGTCTTTCGCGCTCTCTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACt                                                                                                 <  1:173682/133‑1 (MQ=255)                              cgGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCACCGGTAATGAT‑CTCGCGTGCTTCATATTCTTTACTCGAACg                                                                                           >  2:173682/1‑133 (MQ=255)                                   gggTAGTCTTTGGCGTTCTTTAGTATAATGGCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAAt                                                                                      <  1:259051/133‑1 (MQ=255)                                                  ctctTTAATAAAATGCCGATAAAGCAGCATACGTCATAACCTAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGAc                                                                       <  2:277242/133‑1 (MQ=255)                                                   tctTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACg                                                                      <  1:125902/133‑1 (MQ=255)                                                               tGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTc                                                          >  2:151023/1‑133 (MQ=255)                                                               tGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTc                                                          >  2:156174/1‑133 (MQ=255)                                                               tGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTc                                                          >  1:518189/1‑133 (MQ=255)                                                               tGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTc                                                          >  1:604087/1‑133 (MQ=255)                                                                gCCGATAAAGCAGCCTACGTCATACCCGAGCTGCTTCGGGCTGACACCAATACGCCCCCCGGTAAGGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCg                                                         <  2:499172/133‑1 (MQ=255)                                                                 ccGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGc                                                        >  2:617889/1‑133 (MQ=255)                                                                 ccGATAAAGCAGCCTACGTCATAACCGAGCTGCTTC‑GGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGcc                                                       >  1:136000/1‑133 (MQ=255)                                                                             ccTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTa                                            >  2:643936/1‑133 (MQ=255)                                                                                  gTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCgg                                       >  2:75129/1‑133 (MQ=255)                                                                                     aTAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTgc                                    <  1:617889/133‑1 (MQ=255)                                                                                       aaCCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTgcgc                                  >  1:87518/1‑133 (MQ=255)                                                                                          cGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCa                               <  2:604087/133‑1 (MQ=255)                                                                                          cGAGCTGCGTCGGGCAGAAATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCa                               <  1:249806/133‑1 (MQ=255)                                                                                            aGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCAtt                             >  1:585711/1‑133 (MQ=255)                                                                                              ctgctTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCAAAATTGTTCCGCCAGTTCGGCGTAAGCGGTTCGCGCATTGc                           >  2:244424/1‑133 (MQ=255)                                                                                                 ctTCGGGCTGACTTCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCa                        >  1:322895/1‑133 (MQ=255)                                                                                                 ctTCGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGACAGCTCGGCGTAAGCGCTGCGAGCATTGCCCa                        >  2:260952/1‑133 (MQ=255)                                                                                                    tcggCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTa                     <  2:136000/131‑1 (MQ=255)                                                                                                    cGGGCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTa                     <  1:643936/133‑1 (MQ=255)                                                                                                     gggCTGACATCAATACGCGCCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTaa                    <  2:585711/133‑1 (MQ=255)                                                                                                                     gcgcCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGa         >  1:514458/1‑128 (MQ=255)                                                                                                                     gcgcCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGa         <  2:514458/128‑1 (MQ=255)                                                                                                                     gcgcCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCa    <  1:75129/133‑1 (MQ=255)                                                                                                                       gcCCCGGTAATGAT‑CCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACg  >  1:147256/1‑133 (MQ=255)                                                                                                                                     |                                                                                                                        TTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACG  >  NC_000913/3926722‑3926973

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