Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,271,035	(A)7→8	coding (310/1557 nt)	pdeN →	putative c‑di‑GMP phosphodiesterase PdeN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,271,028	1	.	A	100.0%	80.5 / NA	24	D101E (GAT→GAA)	pdeN	putative c‑di‑GMP phosphodiesterase PdeN
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (12/12);  total (12/12)

ACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑T‑AAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCGGCGCTCGACATTAATAAAAACGTCGATATGGCGATCTTACCCGGCACGCCGATGGTGCCG  >  NC_000913/2270904‑2271151                                                                                                                               |                                                                                                                            aCCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑T‑aaaaaaaa                                                                                                                     >  2:284977/1‑132 (MQ=255) aCCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑T‑aaaaaaaa                                                                                                                     >  1:9408/1‑132 (MQ=255)  ccTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCGTGAATGTCAGAAAGTTTGTGCTGGTGAAAGC‑G‑aaaaaaaac                                                                                                                    >  2:224044/1‑131 (MQ=255)              tactGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATc                                                                                                         >  2:42516/1‑133 (MQ=255)                ctGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTg                                                                                                       <  2:356580/133‑1 (MQ=255)                     cGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGAcc                                                                                                  >  2:55767/1‑133 (MQ=255)                      ggCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCg                                                                                                 <  1:42516/133‑1 (MQ=255)                      ggCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCg                                                                                                 <  1:331383/133‑1 (MQ=255)                          ccAGCCGCTGACCTTGGATACCTGCCAGCACGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTga                                                                                             <  1:224044/133‑1 (MQ=255)                             gCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGAt                                                                                          >  1:434898/1‑133 (MQ=255)                                  tGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGAGA‑AAAAAAAAATTGTGTGCATCTGCGGCCGGGGAGATGGACa                                                                                     >  2:130171/1‑133 (MQ=255)                                               cTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAGCGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTGCAGCTGCTACCGGTGAGATGGACATTCCACTCAATGa                                                                        <  2:493360/133‑1 (MQ=255)                                                tGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCa                                                                                <  1:263387/124‑1 (MQ=255)                                                tGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCa                                                                                >  2:263387/1‑124 (MQ=255)                                                tGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGaa                                                                       >  1:394223/1‑133 (MQ=255)                                                tGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGaa                                                                       >  2:489059/1‑133 (MQ=255)                                                            aaCCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCGgcg                                                           <  2:51403/133‑1 (MQ=255)                                                             acaccGAACTGACCGGCCGCGCAGCGTTTAGCATGAATGTCGGAACGTTTGTGTTGGTGATAGA‑TAAAAAAAACTTTGTGTTCATCTTCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCTGgcgc                                                          >  2:380781/4‑133 (MQ=255)                                                              ccccGAACCGACCGCCCGCGCACCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCGGCGCt                                                         <  1:130171/133‑1 (MQ=255)                                                                   caCTGACCGACCGCGCAGCGTCTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGGCCGGTGAGATGGACATTCCACTCAATGAATTGATTCCGGCGCTCGGCa                                                    <  2:398781/132‑1 (MQ=255)                                                                   aaCTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCg                                                              <  2:87580/123‑1 (MQ=255)                                                                   aaCTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCg                                                              >  1:87580/1‑123 (MQ=255)                                                                   aaCTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCGGCGCTCGACa                                                    >  1:424713/1‑133 (MQ=255)                                                                        aCCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCGGCGCTCGACATtaat                                               >  1:94374/1‑133 (MQ=255)                                                                         ccgcccgcGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCGGCGCTCGACATtaata                                              <  2:434898/133‑1 (MQ=255)                                                                                cgcAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCGGCGCTCGACATTAATAAAAACGt                                       <  2:424713/133‑1 (MQ=255)                                                                                  cAGCGTTTAGCATGACTGTCCGAACGTTTGTGCTGGTGAAACA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCGGCGCTCGACATTAATAAAAACGTCg                                     <  2:132450/133‑1 (MQ=255)                                                                                    gCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATACCGGCGCTCGAGATAAATAAATACGTCGat                                   >  1:160196/1‑133 (MQ=255)                                                                                       tttAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATTAATTGATTCCGGAGCTCGACATTAACAAAAACGTCGATATg                                >  2:216932/1‑133 (MQ=255)                                                                                                             ttGTGCTGGTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCGGCGCTCGACATTAATAAAAACGTCGATATGGCGATCTTACCCGACACGCCGa          <  2:94374/133‑1 (MQ=255)                                                                                                                     gTGAAAGA‑TAAAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCGGCGCTCGACATTAATAAAAACGTCGATATGGCGATCTTACCCGGCACGCCGATGGTGCCg  >  2:212709/1‑133 (MQ=255)                                                                                                                               |                                                                                                                            ACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGCCCGCGCAGCGTTTAGCATGAATGTCCGAACGTTTGTGCTGGTGAAAGA‑T‑AAAAAAACATTCTGTTCATCTGCGACCGGTGAGATGGACATTCCACTCAATGAATTGATTCCGGCGCTCGACATTAATAAAAACGTCGATATGGCGATCTTACCCGGCACGCCGATGGTGCCG  >  NC_000913/2270904‑2271151

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: ‑