Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,143,407	T→C	F48L (TTC→CTC)	dgcE →	putative diguanylate cyclase DgcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,143,407	0	T	C	100.0%	130.8 / NA	38	F48L (TTC→CTC)	dgcE	putative diguanylate cyclase DgcE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/12);  total (26/12)

CCCTGCCCCACCCGCTGCTTCACCTGGTCAGTTTAGGTTTAGTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGTTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGACCT  >  NC_000913/2143297‑2143542                                                                                                               |                                                                                                                                        ccctgcccCACCCGCTGCTTCACCTGGTCAGTTTAGGTTTAGTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGttt                                                                                                            >  1:821582/1‑140 (MQ=255) ccctgcccCACCCGCTGCTTCACCTGGTCAGTTTAGGTTTAGTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGttt                                                                                                            >  1:771208/1‑140 (MQ=255)           cccGCTGCTTCACCTGGTCAGTTTAGGTTTAGTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATg                                                                                                  >  2:145819/1‑140 (MQ=255)           cccGCTGCTTCACCTGGACAGTTTAGGTTTAGTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATg                                                                                                  >  2:256886/1‑140 (MQ=255)                tgctTCACCTGGTCAGTTTAGGTTTAGTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCggg                                                                                             >  1:271377/1‑140 (MQ=255)                tgctTCACCTGGTCAGTTTAGGTTTAGTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCggg                                                                                             >  2:581363/1‑140 (MQ=255)                     cACCTGGTCAGTTTAGGTTTAGTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGCGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCAt                                                                                        >  1:156955/1‑140 (MQ=255)                            tCAGTTTAGGTTTAGTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCg                                                                                 >  2:85031/1‑140 (MQ=255)                                tttaggtttagTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGaa                                                                             <  2:271377/140‑1 (MQ=255)                                tttaggtttagTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGaa                                                                             <  1:120811/140‑1 (MQ=255)                                tttaggtttagTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGaa                                                                             >  2:650109/1‑140 (MQ=255)                                     gtttagTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTgcg                                                                        >  2:439625/1‑140 (MQ=255)                                     gtttagTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTgcg                                                                        <  1:650109/140‑1 (MQ=255)                                         agTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGa                                                                    >  2:242832/1‑140 (MQ=255)                                               gTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGtt                                                              <  2:77247/140‑1 (MQ=255)                                                tttATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTc                                                             >  1:778064/1‑140 (MQ=255)                                                tttATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTc                                                             >  1:417241/1‑140 (MQ=255)                                                         aCCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGaaa                                                    >  1:281122/1‑140 (MQ=255)                                                         aCCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGaaa                                                    <  1:242832/140‑1 (MQ=255)                                                              tttCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATcg                                               >  2:185585/1‑140 (MQ=255)                                                               ttCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATcgc                                              >  1:654237/1‑140 (MQ=255)                                                                tctcGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATcgcc                                             >  2:694936/1‑140 (MQ=255)                                                                  tcGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATcgccgc                                           >  2:746135/1‑140 (MQ=255)                                                                         aGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATcgc                                              <  1:695746/130‑1 (MQ=255)                                                                         aGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATcgc                                              >  2:695746/1‑130 (MQ=255)                                                                            tttCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTg                                 <  1:95180/140‑1 (MQ=255)                                                                             ttCGCAGTTTGACACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGc                                <  2:139535/140‑1 (MQ=255)                                                                              tCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCt                               >  1:454556/1‑140 (MQ=255)                                                                              tCGCAGTTTGGCACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGACGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCt                               >  1:212515/1‑140 (MQ=255)                                                                                         cACCCAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTccacca                    <  2:654237/140‑1 (MQ=255)                                                                                            ccAACTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCt                 >  1:699402/1‑140 (MQ=255)                                                                                              aaCTCGCCCCACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCg               <  1:426529/140‑1 (MQ=255)                                                                                                    ccccACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAAc         <  2:821582/140‑1 (MQ=255)                                                                                                     cccACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTc                         <  2:127728/123‑1 (MQ=255)                                                                                                     cccACTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTc                         >  1:127728/1‑123 (MQ=255)                                                                                                        aCTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGa     >  2:271919/1‑140 (MQ=255)                                                                                                        aCTGTGGCTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGa     >  2:681935/1‑140 (MQ=255)                                                                                                           gtgGCTCCCGACGTCCATTATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGACCt  >  2:156723/1‑140 (MQ=255)                                                                                                               |                                                                                                                                        CCCTGCCCCACCCGCTGCTTCACCTGGTCAGTTTAGGTTTAGTCTCGTTTATCTTTACCCTTTTCTCGCTTGAGCTTTCGCAGTTTGGCACCCAACTCGCCCCACTGTGGTTCCCGACGTCCATCATGATGGTGGCGTTTTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGACCT  >  NC_000913/2143297‑2143542

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