Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,093,922	A→G	I349M (ATA→ATG)	dgcT →	putative diguanylate cyclase DgcT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,093,922	0	A	G	100.0%	188.1 / NA	52	I349M (ATA→ATG)	dgcT	putative diguanylate cyclase DgcT
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/26);  total (26/26)

AGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATC  >  NC_000913/1093793‑1094055                                                                                                                                  |                                                                                                                                      agcCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATg                                                                                                                             <  2:328823/140‑1 (MQ=255)  gcCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCg                                                                                                                                     >  2:958966/1‑131 (MQ=255)  gcCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGa                                                                                                                            >  1:956195/1‑140 (MQ=255)   cccAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGAc                                                                                                                           <  2:956195/140‑1 (MQ=255)   cccAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGAc                                                                                                                           >  2:499401/1‑140 (MQ=255)   cccAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCAT‑CGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGa                                                                                                                               <  2:639152/135‑1 (MQ=255)   cccAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCAT‑CGGTTGGCGATCAGGTGATAAAAACAGTGGGGAATATCATTGGGAAAAGCATGCGACCAGa                                                                                                                               >  1:639152/1‑135 (MQ=255)       aaaaaCGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCttt                                                                                                                       >  2:874241/1‑140 (MQ=255)              ccTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGa                                                                                                                               >  2:739832/1‑125 (MQ=255)              ccTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGa                                                                                                                               <  1:739832/125‑1 (MQ=255)                      cGTCATGATTATGGATATCGACCACTACAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTcggcg                                                                                                        <  2:740939/140‑1 (MQ=255)                         cATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGgcggcg                                                                                                     <  1:450749/140‑1 (MQ=255)                          aTGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAGAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGa                                                                                                    <  2:901000/140‑1 (MQ=255)                                 tGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGTTCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGttt                                                                                             >  1:645298/1‑140 (MQ=255)                                       tCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTc                                                                                       >  1:563591/1‑140 (MQ=255)                                        cGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCt                                                                                      >  1:169248/1‑140 (MQ=255)                                        cGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCt                                                                                      >  1:453612/1‑140 (MQ=255)                                          aCCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTg                                                                                    >  2:256825/1‑140 (MQ=255)                                              cTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGa                                                                                >  2:462630/1‑140 (MQ=255)                                                tCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACg                                                                              <  2:583057/140‑1 (MQ=255)                                                     aaaGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACAt                                                                         <  2:453612/140‑1 (MQ=255)                                                     aaaGCCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACAt                                                                         <  2:906233/140‑1 (MQ=255)                                                        gTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGa                                                                      <  1:462630/140‑1 (MQ=255)                                                               acacCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGaa                                                               <  2:717385/140‑1 (MQ=255)                                                                   cTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTg                                                           >  1:783774/1‑140 (MQ=255)                                                                     gggggCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCg                                                         <  2:563591/140‑1 (MQ=255)                                                                     gggggCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCg                                                         >  2:426050/1‑140 (MQ=255)                                                                          cATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGa                                                        <  2:555220/136‑1 (MQ=255)                                                                          cATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGa                                                        >  1:555220/1‑136 (MQ=255)                                                                          cATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGc                                                    >  1:355315/1‑140 (MQ=255)                                                                             ccGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCttt                                                 >  2:420082/1‑140 (MQ=255)                                                                                   ggCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGa                                           <  2:783774/140‑1 (MQ=255)                                                                                       aTCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAgg                                       >  1:107483/1‑140 (MQ=255)                                                                                               aTAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGa                               <  2:169248/140‑1 (MQ=255)                                                                                                tAAAAACAGTGGTGAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGaa                              <  1:256825/140‑1 (MQ=255)                                                                                                         tggtgAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGtt                        <  2:209363/137‑1 (MQ=255)                                                                                                         tggtgAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGtt                        >  1:209363/1‑137 (MQ=255)                                                                                                         tggtgAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAg                     >  2:803319/1‑140 (MQ=255)                                                                                                         tggtgAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAg                     >  1:936039/1‑140 (MQ=255)                                                                                                         tggtgAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAg                     >  1:909618/1‑140 (MQ=255)                                                                                                          ggtgAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGc                    >  2:427949/1‑140 (MQ=255)                                                                                                            tgAATATCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGt                  >  1:41411/1‑140 (MQ=255)                                                                                                                tatCATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTaa              <  2:909618/140‑1 (MQ=255)                                                                                                                   cATTGGGAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTg           >  2:596012/1‑140 (MQ=255)                                                                                                                       gggAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGa       <  1:759703/140‑1 (MQ=255)                                                                                                                       gggAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGa       <  2:936039/140‑1 (MQ=255)                                                                                                                       gggAAAAGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGa       <  1:686402/140‑1 (MQ=255)                                                                                                                          aaaaGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATaa    >  1:774545/1‑140 (MQ=255)                                                                                                                           aaaGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAAt   <  1:874241/140‑1 (MQ=255)                                                                                                                           aaaGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAAt   <  2:107483/140‑1 (MQ=255)                                                                                                                            aaGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATc  <  1:420082/140‑1 (MQ=255)                                                                                                                            aaGCATGCGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATc  <  1:596012/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                      AGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATC  >  NC_000913/1093793‑1094055

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