Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,685,332 A→T L96Q (CTG→CAG)  pdeK ← putative c‑di‑GMP phosphodiesterase PdeK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,685,3320AT94.6% 120.9 / ‑0.1 37L96Q (CTG→CAG) pdeKputative c‑di‑GMP phosphodiesterase PdeK
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (16/19);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

CAGATAGGCAATTGGCTGCGGGTTTGCCGGACGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGA  >  NC_000913/3685202‑3685465
                                                                                                                                  |                                                                                                                                     
cAGATAGGCAATTGGCTGCGGGTTTGCCGGACGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGa                                                                                                                              <  1:294527/140‑1 (MQ=255)
cAGATAGGCAATTGGCTGCGGGTTTGCCGGACGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGa                                                                                                                              <  2:32792/140‑1 (MQ=255)
    tAGGCAATTGGCTGCGGGTTTGCCGGACGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTg                                                                                                                          >  1:318650/1‑140 (MQ=255)
     aGGCAATTGGCTGCGGGTTTGCCGGACGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTgg                                                                                                                         >  2:411814/1‑140 (MQ=255)
       gCAATTGGCTGCGGGTTTGCCGGACGTTCGAGCGAGAAATCGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCATCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGtt                                                                                                                       >  2:624403/1‑140 (MQ=255)
          aTTGGCTGCGGGTTTGCCGGACGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAgg                                                                                                                    <  2:92184/140‑1 (MQ=255)
               cTGCGGGTTTGCCGGACGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAgg                                                                                                                    <  1:45880/135‑1 (MQ=255)
               cTGCGGGTTTGCCGGACGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAgg                                                                                                                    >  2:45880/1‑135 (MQ=255)
                            ggACGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACg                                                                                                  >  2:188234/1‑140 (MQ=255)
                                         agTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCt                                                                                     >  2:36557/1‑140 (MQ=255)
                                          gTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCa                                                                                                    >  2:139346/1‑124 (MQ=255)
                                          gTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCa                                                                                                    <  1:139346/124‑1 (MQ=255)
                                                gCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGtt                                                                              <  1:384459/140‑1 (MQ=255)
                                                    aGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAAt                                                                          >  2:364722/1‑140 (MQ=255)
                                                         gaTTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTgg                                                                     <  1:233116/140‑1 (MQ=255)
                                                                 ccGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATa                                                             <  1:36557/140‑1 (MQ=255)
                                                                    gTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCt                                                          >  1:478080/1‑140 (MQ=255)
                                                                    gTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCt                                                          >  2:153616/1‑140 (MQ=255)
                                                                                    aGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAatat                                          <  2:318650/140‑1 (MQ=255)
                                                                                      tAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCg                                        <  2:565337/140‑1 (MQ=255)
                                                                                      tAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCg                                        <  1:153616/140‑1 (MQ=255)
                                                                                      tAACCATTACCGGCACCGGGCCCTCTGGCATAACACTTTTACGCTGCGCCTGGAACTGGTTAGCCAGCACTACATCGGCACGGCTCACCACGCCTGCCGGTTTAATGCCGGCAAGAATAGCTTCCGCCTCGCCAATATCg                                        >  2:153231/1‑140 (MQ=255)
                                                                                           aTTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCttt                                   >  1:537063/1‑140 (MQ=255)
                                                                                              aCCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCaa                                >  1:568413/1‑140 (MQ=255)
                                                                                               ccGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCaaa                               <  1:364722/140‑1 (MQ=255)
                                                                                                cGGCACCGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCaaaa                              >  2:437731/1‑140 (MQ=255)
                                                                                                    aCCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCg                                        >  1:232515/1‑126 (MQ=255)
                                                                                                    aCCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCg                                        <  2:232515/126‑1 (MQ=255)
                                                                                                     ccGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCt                         <  1:411814/140‑1 (MQ=255)
                                                                                                     ccGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCt                         <  1:369441/140‑1 (MQ=255)
                                                                                                     ccGGGCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCt                         <  2:568413/140‑1 (MQ=255)
                                                                                                       gggCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGa                       <  1:330424/140‑1 (MQ=255)
                                                                                                        ggCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGCATATCGCCTTTCAAAATAGCTGaa                      >  1:569557/1‑140 (MQ=255)
                                                                                                        ggCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCt                         >  1:56173/1‑137 (MQ=255)
                                                                                                        ggCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCt                         <  2:56173/137‑1 (MQ=255)
                                                                                                        ggCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGaa                      >  2:60280/1‑140 (MQ=255)
                                                                                                        ggCGCTCTGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGaa                      >  1:276746/1‑140 (MQ=255)
                                                                                                              ctGGAATAAAACTTTTACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAg                <  1:392698/140‑1 (MQ=255)
                                                                                                                            ttACGCTGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCGCGGCTGACCGCGCCTGCCGGTTTAATGCTGGCAAGGATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACaga  <  1:399561/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                     
CAGATAGGCAATTGGCTGCGGGTTTGCCGGACGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGA  >  NC_000913/3685202‑3685465

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