Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,272,277 T→C W518R (TGG→CGG)  pdeN → putative c‑di‑GMP phosphodiesterase PdeN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,272,2770TC96.0% 74.4 / ‑4.0 25W518R (TGG→CGG) pdeNputative c‑di‑GMP phosphodiesterase PdeN
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (8/16);  total (8/17)
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.70e-01

TGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGTGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCGTTGCAGATGATTGTGCGCATACTGCTGCTGTTTATCGCTCTGTTCACCTTTG  >  NC_000913/2272138‑2272409
                                                                                                                                           |                                                                                                                                    
tGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGAGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGAATTTGTTCTCTGCCTCAAGAACCCGTATACGCCGCAGt                                                                                                                                      >  1:328923/1‑140 (MQ=255)
      aGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAg                                                                                                                                <  1:605600/140‑1 (MQ=255)
         ggTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGgtg                                                                                                                             <  1:362773/140‑1 (MQ=255)
         ggTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGATGGCTAAGGAATCCATATACGTCGCAGCGGTAAGgtg                                                                                                                             >  1:264585/1‑140 (MQ=255)
                aCGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTAc                                                                                                                      <  2:461776/140‑1 (MQ=255)
                   ccGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGTGGTAAGGTGTGCTTACGTc                                                                                                                   <  2:124466/140‑1 (MQ=255)
                   ccGGAACAGGCGCGATGCCTAAGCGAACGCGTCGTGAATCTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTc                                                                                                                   <  2:251995/140‑1 (MQ=255)
                    cGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTcc                                                                                                                  >  2:559225/1‑140 (MQ=255)
                        aCAGGCGCGATGGCTAAGCGAACGTGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCtta                                                                                                              <  2:31681/140‑1 (MQ=255)
                         cAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCttat                                                                                                             <  2:106517/140‑1 (MQ=255)
                              gcgATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTc                                                                                                        <  1:571650/140‑1 (MQ=255)
                                       aaGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGaaa                                                                                               <  1:572290/140‑1 (MQ=255)
                                              cgcgGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGcc                                                                                        >  1:43444/1‑140 (MQ=255)
                                                  gcgTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCGTACGTCCCTTAGTAATCATAGTGAAAGCATGCCGGAt                                                                                    >  1:324583/1‑140 (MQ=255)
                                                                         tGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAgg                                                             <  2:359441/140‑1 (MQ=255)
                                                                                gccgccCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCg                                                      <  1:611446/140‑1 (MQ=255)
                                                                                gccgccCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCAGAGTGAAAGGATGGCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCg                                                      >  2:229804/1‑140 (MQ=255)
                                                                                           ccGCTGGATGATTTTGTACGCTGGCTAAAGCGACCGTATGCGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCGTTGCAgatgat                                           <  2:324583/140‑1 (MQ=255)
                                                                                           ccGCTGGACGATTTTGTTCGCTTGCTAAAGATATCGTATGCCCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCGTTGCAgatgat                                           <  2:303694/140‑1 (MQ=255)
                                                                                              cTGGACGATTTTGTTCGCTGGCTAATGAATCCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCGTTGCAGATGATtgt                                        <  1:229804/140‑1 (MQ=255)
                                                                                                 gACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCGTTGCAGATGATTGTgcg                                     >  1:642476/1‑140 (MQ=255)
                                                                                                   cGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCGTTGCAGATGATTGTGCGCa                                   <  2:642476/140‑1 (MQ=255)
                                                                                                                tGGCTAAAGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCGTTGCAGTTGATTGTGCACATACTGCTGCTGtt                      >  1:286771/1‑140 (MQ=255)
                                                                                                                     aaaGAAACCGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCGTTGCAGATGATTGTGCGCATACTGCTGCTGTTTATCg                 <  1:459676/140‑1 (MQ=255)
                                                                                                                            ccGTATACGCCGCAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCGTTGCAGATGATTGTGCGCATACTGCTGCTGTTTATCGCTCTGtt          <  1:559225/140‑1 (MQ=255)
                                                                                                                                    gccgcAGCGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCGTTGCAGATGATTGTGCGCATACTGCTGCTGTTTATCGCTCTGTTCACCTTTg  >  2:61921/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                    
TGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGTGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTGCGGCTAATGATGAGTAAAAGGAAATCCGTTGCAGATGATTGTGCGCATACTGCTGCTGTTTATCGCTCTGTTCACCTTTG  >  NC_000913/2272138‑2272409

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