Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,402,027 A→G 17.4% L9P (CTG→CCG)  ydaM ← diguanylate cyclase, csgD regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,402,0270AG17.4% 47.1 / 3.8 23L9P (CTG→CCG) ydaMdiguanylate cyclase, csgD regulator
Reads supporting (aligned to +/- strand):  ref base A (11/8);  new base G (2/2);  total (13/10)
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.29e-01

TTCTTTTTGTGAGGAGACGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAA  >  CP009273/1401882‑1402175
                                                                                                                                                 |                                                                                                                                                    
ttCTTTTTGTGAGGAGACGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAggg                                                                                                                                                   >  1:2513/1‑149 (MQ=255)
                  ggAAGAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTAGCGTCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTaa                                                                                                                                 <  1:143569/149‑1 (MQ=255)
                                           cTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTAtt                                                                                                              >  2:269496/1‑143 (MQ=255)
                                                     ggTCCTGCAGCCACAGCTTCGCCGCACTATCTGCACCAGCTAACTGTGCCTCAAAGTGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATCGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAAtt                                                                                                    <  2:145028/143‑1 (MQ=255)
                                                      gTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTAATCAAttt                                                                                                   >  2:20236/1‑143 (MQ=255)
                                                        cTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAAttttt                                                                                                 <  2:2086/143‑1 (MQ=255)
                                                                aaCAGTTTCCCTGCGCCATTCGCAGAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGGAATTCCAGGGTATTGAAGTTGTGCTTAATCATTGAGATCCCGTTATTTTGATCAATTTTTGTTGCTATCCGATc                                                                                   <  1:22561/149‑1 (MQ=255)
                                                                     tttCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGAt                                                                                    >  2:261573/1‑143 (MQ=255)
                                                                      ttCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATc                                                                                   <  2:83467/143‑1 (MQ=255)
                                                                                     gCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCGGGGTATTGCCGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATg                                                              >  1:260516/1‑149 (MQ=255)
                                                                                     gCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCGGGGTATTGAAGTTGTGCGTCATCATCGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCTTATg                                                              >  1:96156/1‑149 (MQ=255)
                                                                                          aaTTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCatat                                                               >  2:93881/1‑143 (MQ=255)
                                                                                                              gCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGa                                     >  1:71010/1‑149 (MQ=255)
                                                                                                                      aTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTAttt                             <  1:56473/149‑1 (MQ=255)
                                                                                                                      aTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTAttt                             >  1:216327/1‑149 (MQ=255)
                                                                                                                      aTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTGGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACg                                   >  2:136915/1‑143 (MQ=255)
                                                                                                                         cAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACgtgg                                <  2:214262/143‑4 (MQ=255)
                                                                                                                                        aGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCacac                 >  2:24601/1‑143 (MQ=255)
                                                                                                                                        aGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCacac                 >  2:153534/1‑143 (MQ=255)
                                                                                                                                        aGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTg           >  1:21727/1‑149 (MQ=255)
                                                                                                                                         gTAAGTCCGGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACg                <  2:85029/143‑1 (MQ=255)
                                                                                                                                         gTAAGTCCGGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATAATAGACGCTATTTAACACTTCACACg                <  2:233145/143‑1 (MQ=255)
                                                                                                                                         gTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGt          <  1:180063/149‑1 (MQ=255)
                                                                                                                                                 aGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCCCGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGaa  <  1:106329/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTCTTTTTGTGAGGAGACGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAA  >  CP009273/1401882‑1402175

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑