Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,899,621 G→A V399I (GTT→ATT)  pdeD → putative c‑di‑GMP phosphodiesterase PdeD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,899,6210GA97.8% 146.0 / NA 46V399I (GTT→ATT) pdeDputative c‑di‑GMP phosphodiesterase PdeD
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (26/19);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATGTTGGCATTAACGTCGCACCCAGCCATTTTCGCCGTGGTGTGCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATGTTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTGAAACATTACGTCCTGACGTGCTGAAAATTGATAAGTCA  >  NC_000913/1899484‑1899752
                                                                                                                                         |                                                                                                                                   
ttCATGTTGGCATTAACGTCGCACCCAGCCATTTTCGCCGTGGTGTGCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATAtt                                                                                                                                   >  2:886894/1‑140 (MQ=255)
         gCATTAACGTCGCACCCAGCCATTTTCGCCGTGGTGTGCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCgg                                                                                                                          >  1:923898/1‑140 (MQ=255)
                   cGCACCCAGCCATTTTCGCCGTGGTGTGCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCcgcg                                                                                                                >  1:242150/1‑140 (MQ=255)
                    gCACCCAGCCATTTTCGCCGTGATGTGCTTATAAAAGGTCTCAATCAGTACTGGTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAg                                                                                                              <  1:442521/140‑1 (MQ=255)
                     cACCCAGCCATTTTCGCCGTGGTGTGCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAg                                                                                                              <  1:886894/140‑1 (MQ=255)
                     cACCCAGCCATTTTCGCCGTGGTGTGCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAg                                                                                                              <  1:782816/140‑1 (MQ=255)
                           gCCATTTTCGCCGTGGTGTGCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCAt                                                                                                        >  1:990633/1‑140 (MQ=255)
                                    gCCGTGGTGTGCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAAc                                                                                               >  1:563569/1‑140 (MQ=255)
                                      cGTGGTGTGCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCTGCGAGCTGCATCGTAAAAACGt                                                                                             >  2:670784/1‑140 (MQ=255)
                                        tggtgTGCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCa                                                                                           >  2:228363/1‑140 (MQ=255)
                                             tgCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTg                                                                                      >  2:58212/1‑140 (MQ=255)
                                              gCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTgg                                                                                     <  1:540419/140‑1 (MQ=255)
                                               cTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGc                                                                                    <  2:519429/140‑1 (MQ=255)
                                                  aTAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGAt                                                                                 >  2:577206/1‑140 (MQ=255)
                                                       aGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATtgatg                                                                            >  2:8965/1‑140 (MQ=255)
                                                            tcAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTc                                                                       >  2:805411/1‑140 (MQ=255)
                                                            tcAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTc                                                                       >  1:41601/1‑140 (MQ=255)
                                                            tcAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTc                                                                       >  2:636690/1‑140 (MQ=255)
                                                              aaTCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCgg                                                                     >  2:549252/1‑140 (MQ=255)
                                                                 cAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCTGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCAc                                                                  <  1:670784/140‑1 (MQ=255)
                                                                    tACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCgg                                                               >  2:845721/1‑140 (MQ=255)
                                                                     aCTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGc                                                              <  1:577206/140‑1 (MQ=255)
                                                                          tttAGCGCTCACCCGACCCAGCTACTGATCCACAAAATCTCCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAg                                                         <  2:447331/140‑1 (MQ=255)
                                                                              gcgcTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGGTTGTTGACTTCGGGACCGGCAACAGCTCg                                                     >  1:447331/1‑140 (MQ=255)
                                                                              gcgcTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCAcc                                                                 >  2:919300/1‑128 (MQ=255)
                                                                              gcgcTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCAcc                                                                 <  1:919300/128‑1 (MQ=255)
                                                                              gcgcTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCg                                                     >  2:821007/1‑140 (MQ=255)
                                                                              gcgcTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCg                                                     >  2:236305/1‑140 (MQ=255)
                                                                              gcgcTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCg                                                     >  1:26700/1‑140 (MQ=255)
                                                                                   cACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAGCGTCAAACAGGCTATTGATGACTTCGGCACCGGCAACAGCTCGTTTTc                                                >  2:383249/1‑140 (MQ=255)
                                                                                   cACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTc                                                >  2:658586/1‑140 (MQ=255)
                                                                                   cACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTc                                                >  2:999902/1‑140 (MQ=255)
                                                                                         aTTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCt                                          <  2:357944/140‑1 (MQ=255)
                                                                                         aTTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCt                                          >  2:740835/1‑140 (MQ=255)
                                                                                             aGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTGaa                                      <  1:58212/140‑1 (MQ=255)
                                                                                                  cTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCGCCGGCAACAGCTCGTTTTCCTGGCTTGAAACAtt                                 <  1:979981/140‑1 (MQ=255)
                                                                                                  cTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTGAAACAtt                                 <  1:805411/140‑1 (MQ=255)
                                                                                                  cTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTGAAACAtt                                 <  1:636690/140‑1 (MQ=255)
                                                                                                            aaaTCACCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTGAAACATTACGTCCTGAc                       <  2:242150/140‑1 (MQ=255)
                                                                                                                cACCGAACGCGATGCCTTACTGCATGTTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTGAAACATTACGTCCTGACGTGc                   <  1:383249/140‑1 (MQ=255)
                                                                                                                 aCCGAACGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTGAAACATTACGTCCTGACGTGCt                  >  2:84365/1‑140 (MQ=255)
                                                                                                                  ccGAACGCGATGCCTTACTGGACATTGATTATCGGATTGCCCGCGCGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTTAAACATTACGTCATGACGTGCTg                 >  2:870621/1‑140 (MQ=255)
                                                                                                                   cGAACGCGATGCGTTACTGGATGTTGATTATGGGACTGCCCGCGAGCTCCATCGTAAAAACGTCAAGCTTGCGATTGATAACTTCCGCACCGGCAACAGCTCGGTTTCCTGGCTTGAAACATTACGTCCTGACGTGCTGa                <  2:428095/140‑1 (MQ=255)
                                                                                                                      aCGCGATGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTGAAACATTACGTCCTGACGTGCTGaaaa             <  1:8965/140‑1 (MQ=255)
                                                                                                                           aTGCCTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTGAAACATTACGTCCTGACGTGCTGAAAATTGAt        <  1:821007/140‑1 (MQ=255)
                                                                                                                              ccTTACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTGAAACATTACGTCCTGACGTGCTGAAAATTGATAAg     <  2:984722/140‑1 (MQ=255)
                                                                                                                                 tACTGGATATTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTGAAACATTACGTCCTGACGTGCTGAAAATTGATAAGTCa  <  2:645835/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                   
TTCATGTTGGCATTAACGTCGCACCCAGCCATTTTCGCCGTGGTGTGCTGATAAAAGATCTCAATCAGTACTGGTTTAGCGCTCACCCGATTCAGCAACTGATCCTCGAAATCACCGAACGCGATGCCTTACTGGATGTTGATTATCGGATTGCCCGCGAGCTGCATCGTAAAAACGTCAAACTGGCGATTGATGACTTCGGCACCGGCAACAGCTCGTTTTCCTGGCTTGAAACATTACGTCCTGACGTGCTGAAAATTGATAAGTCA  >  NC_000913/1899484‑1899752

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: