Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,143,684 (G)6→7 coding (419/3318 nt) dgcE → putative diguanylate cyclase DgcE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,143,6781.G93.5% 98.0 / 0.3 31coding (413/3318 nt)dgcEputative diguanylate cyclase DgcE
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (19/10);  total (20/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATC  >  NC_000913/2143545‑2143809
                                                                                                                                       |                                                                                                                                   
acgacCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCAGTCGGCAGCGCCATTGTTCCGCCTCTGTTgggggg                                                                                                                               <  2:358189/140‑1 (MQ=255)
   acCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTgggggggtg                                                                                                                            >  2:748944/1‑137 (MQ=255)
   acCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTgggggggtg                                                                                                                            >  2:626920/1‑137 (MQ=255)
     cATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTgggggggtgtt                                                                                                                          <  2:126615/140‑6 (MQ=255)
         aaTATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTgg                                                                                                                       >  2:62204/1‑140 (MQ=255)
         aaTATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTgg                                                                                                                       >  2:446522/1‑140 (MQ=255)
               gttgAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTcc                                                                                                                 >  1:314032/1‑140 (MQ=255)
               gttgAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTGCTGGTTGTcc                                                                                                                 >  1:554905/1‑140 (MQ=255)
                    aGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCctgctg                                                                                                            >  2:445890/1‑140 (MQ=255)
                     gCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGa                                                                                                           >  2:408249/1‑140 (MQ=255)
                     gCCGTGCTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTTCACAATCTGCCTGTCTGGCTGCCTCTCGCACTCGGCAGCGCCGTTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGa                                                                                                           >  2:30626/1‑140 (MQ=255)
                            tCGGGGCGGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGga                                                                                                   <  2:801041/140‑1 (MQ=255)
                                           tACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAg                                                                                     <  1:438864/140‑1 (MQ=255)
                                            aCGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAgg                                                                                    <  1:748944/140‑1 (MQ=255)
                                             cGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAggg                                                                                   <  1:62204/140‑1 (MQ=255)
                                               tAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCa                                                                                 >  1:254863/1‑140 (MQ=255)
                                                      ctgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGa                                                                          >  2:430490/1‑140 (MQ=255)
                                                       tgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGat                                                                         >  1:523249/1‑140 (MQ=255)
                                                        gcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGTCCGGAGACGATCCTCTCAGGGCATTTTTGata                                                                        <  2:554905/140‑1 (MQ=255)
                                                                        ttGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAAt                                                        >  1:58646/1‑140 (MQ=255)
                                                                         tGCAAAATCTGGCTGACTGGCTGCGTCTGGTACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                                       >  1:485203/1‑140 (MQ=255)
                                                                         tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                                       >  2:154311/1‑140 (MQ=255)
                                                                         tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                                       >  1:782225/1‑140 (MQ=255)
                                                                         tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                                       >  1:345888/1‑140 (MQ=255)
                                                                         tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                                       >  2:353041/1‑140 (MQ=255)
                                                                               aTCTGGCTGACTGGCTGCGTCTGGTACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCgg                                                 <  2:485203/140‑1 (MQ=255)
                                                                               aTCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCgg                                                 <  1:445890/140‑1 (MQ=255)
                                                                                          tGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCAc                                      >  2:40152/1‑140 (MQ=255)
                                                                                                  tCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCg                              >  2:714940/1‑140 (MQ=255)
                                                                                                   cgggCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGc                             <  2:340333/138‑1 (MQ=255)
                                                                                                   cTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGc                             <  2:314032/140‑1 (MQ=255)
                                                                                                          tCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGAtt                      <  1:446522/140‑1 (MQ=255)
                                                                                                                    cATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTaaa            <  1:408249/140‑1 (MQ=255)
                                                                                                                    cATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTaaa            >  1:127826/1‑140 (MQ=255)
                                                                                                                              ccTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATc  >  2:420871/1‑140 (MQ=255)
                                                                                                                              ccTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATc  >  1:185570/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                   
ACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATC  >  NC_000913/2143545‑2143809

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