Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,742,721	C→T	G113G (GGC→GGT)	dgcN →	diguanylate cyclase DgcN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,742,721	0	C	T	100.0%	191.3 / NA	56	G113G (GGC→GGT)	dgcN	diguanylate cyclase DgcN
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/30);  total (26/30)

TTTTTGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTT  >  NC_000913/2742588‑2742855                                                                                                                                      |                                                                                                                                       tttttGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACt                                                                                                                                  <  1:490863/140‑1 (MQ=255)   tttGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACttt                                                                                                                                <  2:819894/140‑1 (MQ=255)   tttGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACttt                                                                                                                                <  2:300566/140‑1 (MQ=255)           tGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAAtt                                                                                                                        <  2:50208/140‑1 (MQ=255)                cTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGt                                                                                                                            >  2:366028/1‑131 (MQ=255)                cTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGt                                                                                                                            <  1:366028/131‑1 (MQ=255)                     gcaACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCa                                                                                                              <  2:1056637/140‑1 (MQ=255)                      caACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCAc                                                                                                             <  1:882466/140‑1 (MQ=255)                        aCTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTg                                                                                                           >  1:824115/1‑140 (MQ=255)                        aCTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTg                                                                                                           >  2:642417/1‑140 (MQ=255)                          tGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGc                                                                                                         <  2:1448466/140‑1 (MQ=255)                                acTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGctc                                                                                                       <  2:938053/136‑1 (MQ=255)                                acTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGctc                                                                                                       >  1:938053/1‑136 (MQ=255)                                  tGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCt                                                                                                 <  1:713259/140‑1 (MQ=255)                                   ggCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTg                                                                                                >  2:710713/1‑140 (MQ=255)                                    gCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGc                                                                                               >  1:903034/1‑140 (MQ=255)                                    gCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGc                                                                                               <  1:642417/140‑1 (MQ=255)                                          cTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGct                                                                                                        <  2:1089670/125‑1 (MQ=255)                                          cTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGct                                                                                                        >  1:1089670/1‑125 (MQ=255)                                             ggCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCcatcat                                                                                      >  2:406727/1‑140 (MQ=255)                                              gCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCAtt                                                                                     <  1:1190197/140‑1 (MQ=255)                                              gCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCAtt                                                                                     >  1:114872/1‑140 (MQ=255)                                                   cAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGccc                                                                                             <  2:294730/127‑1 (MQ=255)                                                   cAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGccc                                                                                             >  1:294730/1‑127 (MQ=255)                                                     aGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGa                                                                              <  1:710713/140‑1 (MQ=255)                                                       ggCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGAtt                                                                            <  2:421736/140‑1 (MQ=255)                                                                              cGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACaa                                                                    <  1:1023618/125‑1 (MQ=255)                                                                              cGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTgg                                                     >  2:796250/1‑140 (MQ=255)                                                                              cGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTgg                                                     >  2:914353/1‑140 (MQ=255)                                                                              cGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAAGTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACaa                                                                    >  2:1023618/1‑125 (MQ=255)                                                                                tGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCg                                                   <  1:759737/140‑1 (MQ=255)                                                                                tGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCACCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCg                                                   <  1:551220/140‑1 (MQ=255)                                                                                  aTAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGaa                                                 <  2:58593/140‑1 (MQ=255)                                                                                   tAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAg                                                >  2:1292545/1‑140 (MQ=255)                                                                                        agcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCg                                                   <  2:1218899/132‑1 (MQ=255)                                                                                        agcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCg                                                   >  1:1218899/1‑132 (MQ=255)                                                                                          caAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTgg                                                     >  2:945094/1‑128 (MQ=255)                                                                                          caAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTgg                                                     <  1:945094/128‑1 (MQ=255)                                                                                            aaaTATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTaa                                       <  2:745523/140‑1 (MQ=255)                                                                                               tatTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCg                                    <  2:1136151/140‑1 (MQ=255)                                                                                                 ttCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCt                                  <  2:824115/140‑1 (MQ=255)                                                                                                     ggCATCCTGGCATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTcgcg                              <  1:1392586/140‑1 (MQ=255)                                                                                                               cATTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAAt                    <  1:1292545/140‑1 (MQ=255)                                                                                                                aTTACACCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATc                   >  2:1269151/1‑140 (MQ=255)                                                                                                                    cacCCGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGcc               >  1:1321029/1‑140 (MQ=255)                                                                                                                        cGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGa                             >  1:826143/1‑122 (MQ=255)                                                                                                                        cGTAAGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGa                             <  2:826143/122‑1 (MQ=255)                                                                                                                           aaGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTAt        <  1:1269151/140‑1 (MQ=255)                                                                                                                           aaGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTAt        >  2:769528/1‑140 (MQ=255)                                                                                                                            aGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGcc               >  1:735611/1‑132 (MQ=255)                                                                                                                            aGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGcc               <  2:735611/132‑1 (MQ=255)                                                                                                                            aGGATCCAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATc       >  1:140804/1‑140 (MQ=255)                                                                                                                                 ccAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGtt  >  1:1431091/1‑140 (MQ=255)                                                                                                                                 ccAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGtt  >  2:1448032/1‑140 (MQ=255)                                                                                                                                 ccAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGtt  >  2:1466912/1‑140 (MQ=255)                                                                                                                                 ccAGGTGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGtt  >  2:287453/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                       TTTTTGCCGATGGCCCTGCAGCAACTGAAACACTGGCAGCGCTGGGCCAGCAAGGGCAATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTT  >  NC_000913/2742588‑2742855

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