Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,683,736 Δ1 bp coding (1883/1989 nt) yhjK ← cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,683,7350T.97.5% 178.5 / ‑0.8 40coding (1884/1989 nt)yhjKcyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (18/21);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

AAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCA  >  NC_000913/3683589‑3683874
                                                                                                                                                  |                                                                                                                                           
aaaCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGtgtgc                                                                                                                                         >  1:101422/1‑147 (MQ=255)
aaaCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGtgtgc                                                                                                                                         >  2:84139/1‑147 (MQ=255)
 aaCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCt                                                                                                                                       >  2:357047/1‑151 (MQ=255)
    cAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTc                                                                                                                                     >  2:52555/1‑150 (MQ=255)
    cAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCa                                                                                                                                    <  2:25736/151‑1 (MQ=255)
     aGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTc                                                                                                                                     >  2:108205/1‑149 (MQ=255)
             ttAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCAc                                                                                                                            <  2:120649/150‑1 (MQ=255)
              tAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGc                                                                                                                          >  1:368744/1‑151 (MQ=255)
                   ttttGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCg                                                                                                                     >  1:58927/1‑151 (MQ=255)
                      tGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGc                                                                                                                   <  2:103520/150‑1 (MQ=255)
                         aaTCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATc                                                                                                               >  2:423391/1‑151 (MQ=255)
                               tttGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTg                                                                                                          <  2:474564/150‑1 (MQ=255)
                                      tAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTc                                                                                                   >  2:200874/1‑150 (MQ=255)
                                          tACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGc                                                                                               >  1:138069/1‑150 (MQ=255)
                                              tttCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTgc                                                                                          >  2:68711/1‑151 (MQ=255)
                                                tcttcCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTgcgc                                                                                        >  1:37198/1‑151 (MQ=255)
                                                  ttcCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGcc                                                                                       >  1:166382/1‑150 (MQ=255)
                                                   tcCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCa                                                                                      <  1:50549/150‑1 (MQ=255)
                                                    ccAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAg                                                                                     <  2:190387/150‑1 (MQ=255)
                                                    ccAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGc                                                                                    >  1:151916/1‑151 (MQ=255)
                                                      aGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCAt                                                                                  <  1:503151/151‑1 (MQ=255)
                                                        gTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATg                                                                                 >  1:452164/1‑150 (MQ=255)
                                                                   cGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTgcagc                                                                      >  1:144429/1‑150 (MQ=255)
                                                                           tCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCAtgc                                                              <  2:28849/150‑1 (MQ=255)
                                                                             aaTAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCAtgctgc                                                           >  1:495007/1‑151 (MQ=255)
                                                                                aGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTa                                                         <  1:236209/150‑1 (MQ=255)
                                                                                          cGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCaa                                              <  2:383411/151‑1 (MQ=255)
                                                                                              cAAAAAGGAAGCCCTGGGCAATACCAACGCCCACTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCt                                          >  1:50946/1‑151 (MQ=255)
                                                                                               aaaaaGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCtt                                         <  2:226660/151‑1 (MQ=255)
                                                                                                 aaaGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTca                                       <  2:38953/151‑1 (MQ=255)
                                                                                                     gAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTcaaca                                    >  2:121066/1‑150 (MQ=255)
                                                                                                             gggCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTg                           <  2:577847/151‑1 (MQ=255)
                                                                                                              ggCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGt                          <  1:503129/151‑1 (MQ=255)
                                                                                                                cAATACCAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCg                        <  2:217745/151‑1 (MQ=255)
                                                                                                                     ccAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGAtttt                   <  1:356981/151‑1 (MQ=255)
                                                                                                                     ccAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGAtttt                   >  2:555326/1‑151 (MQ=255)
                                                                                                                      cAACGCCCGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTc                  <  2:149191/151‑1 (MQ=255)
                                                                                                                           cccGCTTTTGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTAc             <  1:499098/151‑1 (MQ=255)
                                                                                                                              gCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGt           <  2:527930/151‑1 (MQ=255)
                                                                                                                                 tttGCCAGCCAGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCg         <  2:101436/149‑1 (MQ=255)
                                                                                                                                    gccagccaGTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTg     <  2:397296/150‑1 (MQ=255)
                                                                                                                                      cagccagTCGCG‑TGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCa  <  2:253501/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
AAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCA  >  NC_000913/3683589‑3683874

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑