Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,853,344 A→G 78.8% intergenic (‑13/‑305) PP_2505 ← / → PP_2506 GAF domain/GGDEF domain‑containing protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,853,3440AG78.8% 61.0 / 11.1 33intergenic (‑13/‑305)PP_2505/PP_2506GAF domain/GGDEF domain‑containing protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base G (12/14);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCTGGGCTAGTTCTTGGGTGAGTCGTACCAAGGTATCAAGATAATGCTCGGCCGGCGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCACTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAGCAAACCGGTTGGCTCTATCCCGGATCCTCCAG  >  NC_002947/2853200‑2853486
                                                                                                                                                |                                                                                                                                              
cACCTGGGCTAGTTCTTGGGTGAGTCGTACCAAGGTATCAAGATAATGCTCGGCCGGCGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGTACAAGCTGGCATGTCTGTTactccactcc                                                                                                                                            >  7:125298/1‑149 (MQ=255)
         tAGTTCTTGGGTGAGTCGTACCAAGGTATCAAGATAATGCTCGGCCGGCGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCtt                                                                                                                                   >  1:23819/1‑149 (MQ=255)
           gTTCTTGGGTGAGTCGTACCAAGGTATCAAGATAATGCTCGGCCGGCGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCttt                                                                                                                                  <  6:171853/148‑1 (MQ=255)
            ttCTTGGGTGAGTCGTACCAAGGTATCAAGATAATGCGCGGCCGGCGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCttttt                                                                                                                                <  8:91819/149‑1 (MQ=255)
                          gTACCAAGGTATCAAGATAATGCTCGGCCGGCGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGGTCTTTTTTTCACAAATGCGGa                                                                                                                  <  5:133536/149‑1 (MQ=255)
                                        aGATAATGCTCGGCCGGCGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCaa                                                                                                    <  4:68744/149‑1 (MQ=255)
                                           tAATGCTCGGCCGGCGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCa                                                                                                 >  8:123698/1‑149 (MQ=255)
                                                cTCGGCCGGCGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGATTCATCAGCGGGTAGGGGACAAGCAGGCATGTCTGTTACTCCACTACTTTGTTCTTTTTTTGACAAATGCGGATAGAAAATTCTAAAGCATCtt                                                                                             >  7:151857/1‑148 (MQ=255)
                                                  cggccggaGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCaa                                                                                                    <  5:142177/139‑1 (MQ=255)
                                                  cggccggaGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCaa                                                                                                    >  6:142177/1‑139 (MQ=255)
                                                        gcGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGc                                                                                    <  1:71400/149‑1 (MQ=255)
                                                        gcGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGc                                                                                    <  3:130810/149‑1 (MQ=255)
                                                                    aGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCg                                                                        >  8:216297/1‑149 (MQ=255)
                                                                      tCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCACTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGgc                                                                      >  5:176329/1‑149 (MQ=255)
                                                                         atatCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGcc                                                                   <  4:179161/149‑1 (MQ=255)
                                                                         atatCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCACTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGcc                                                                   >  3:192968/1‑149 (MQ=255)
                                                                              atcaAGTGTGCGCTGACGAAGTCCTTCATCACCGGGTACCGGACAAGCTGCCATGTCAGTTACTCCGCTCCTTGGATCTATTTTTCACAAAAGCGGATAGAAAATTCTCAAGCAACTTGACCAATGCCGCACAACCTCGCCCCCTTcat                                                              >  5:307966/1‑148 (MQ=255)
                                                                               tcaAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTcagg                                                             >  4:128491/1‑149 (MQ=255)
                                                                                                aaGTTCTTCATCAGCTGGTAGGGGAGAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCTCACGCTCGGCGCCTTCAGGCCAGGCCAGCCTGAATa                                            <  8:151857/149‑1 (MQ=255)
                                                                                                   tggttgatcaGCGGGTAGGTGTGAATCTGGCTTGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGAGCAATGCCGCACACGCTCGGCGCGTTCAGGCCAGGCCAGCCGGAATAggg                                         <  6:307966/143‑1 (MQ=255)
                                                                                                     cTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCAAGCCGGAATAgg                                          >  6:58215/1‑146 (MQ=255)
                                                                                                     cTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCACTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAg                                       >  4:199105/1‑149 (MQ=255)
                                                                                                          tcaGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAg                                  >  1:226489/1‑149 (MQ=255)
                                                                                                                ggTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAGCAAAcc                            >  7:289690/1‑149 (MQ=255)
                                                                                                                ggTAGGGGACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAGCAAAcc                            >  8:138451/1‑149 (MQ=255)
                                                                                                                    ggggACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAGCAAACCGGtt                        <  7:138451/149‑1 (MQ=255)
                                                                                                                      ggACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAGCAAACCGGTTgg                      <  7:123698/149‑1 (MQ=255)
                                                                                                                       gACAAGCTGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAGCAAACCGGTTGGc                     >  7:189250/1‑149 (MQ=255)
                                                                                                                       gACAAGCTGGCATGTCTGTTACTCCACTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAGCAAACCGGTTGGc                     <  6:176329/149‑1 (MQ=255)
                                                                                                                              tGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGa                                        >  7:104971/1‑123 (MQ=255)
                                                                                                                              tGGCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGa                                        <  8:104971/123‑1 (MQ=255)
                                                                                                                                gCATGTCTGTTACTCCGCTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAGCAAAc                             <  8:66028/132‑1 (MQ=255)
                                                                                                                                gCATGTCTGTTACTCCACTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAGCAAACCGGTTGGCTCTATccc             <  3:199105/148‑1 (MQ=255)
                                                                                                                                   tGTCTGTTACTACACTCCATGGTTCTTTTTTTCACAAATGCTGATAGAAAATTCTCAAGCATCTTGACCATTGCCGCACACGCTCTGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAGCAAACGGGTTGGCTCTATCCCGGa          <  8:125298/148‑1 (MQ=255)
                                                                                                                                           actccgctccTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAGCAAACCGGTTGGCTCTATCCCGGATCCTCCAg  <  8:289690/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CACCTGGGCTAGTTCTTGGGTGAGTCGTACCAAGGTATCAAGATAATGCTCGGCCGGCGTATCGACCAAGTCCATATCATCAAGTGTGCGCTGACGAAGTGCTTCATCAGCGGGTAGGGGACAAGCTGGCATGTCTGTTACTCCACTCCTTGGTTCTTTTTTTCACAAATGCGGATAGAAAATTCTCAAGCATCTTGACCAATGCCGCACACGCTCGGCGCCTTCAGGCCAGGCCAGCCGGAATAGGGAGCATAGCAAACCGGTTGGCTCTATCCCGGATCCTCCAG  >  NC_002947/2853200‑2853486

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: