Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,699,544 C→A 8.7% P63H (CCC→CAC)  PP_1494 → GGDEF domain‑containing response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,699,5440CA 8.7% 141.6 / 2.9 46P63H (CCC→CAC) PP_1494GGDEF domain‑containing response regulator
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base A (2/2);  total (23/23)
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.63e-01

CAACGAAAACTCGGCGATGGTGCTGCTGGTCGACGATCAGGCCATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAGGAAGACCCGCTGGTGAAGAGTGCGGCATTTGCC  >  NC_002947/1699395‑1699692
                                                                                                                                                     |                                                                                                                                                    
caACGAAAACTCGGCGATGGTGCTGCTGGTCGACGATCAGGCCATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGcc                                                                                                                                                      >  1:672309/1‑150 (MQ=255)
caACGAAAACTCGGCGATGGTGCTGCTGGTCGACGATCAGGCCATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGcc                                                                                                                                                      >  1:501746/1‑150 (MQ=255)
    gAAAACTCGGCGATGGTGCTGCTGGTCGACGATCAGTCCATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACg                                                                                                                                                  <  1:970907/150‑1 (MQ=255)
          tCGGCGATGGTGCTGCTGGTCGACGATCAGGCCATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCAtc                                                                                                                                            <  2:882896/150‑2 (MQ=255)
                             tCGACGATCAGGCCATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGc                                                                                                                         <  1:517854/150‑1 (MQ=255)
                                acgaTCAGGCCATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTg                                                                                                                      <  1:165888/150‑1 (MQ=255)
                                acgaTCAGGCCATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTg                                                                                                                      <  1:166379/150‑1 (MQ=255)
                                  gaTCAGGCCATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTggg                                                                                                                    >  2:900811/1‑150 (MQ=255)
                                  gaTCAGGCCATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTggg                                                                                                                    <  1:565396/150‑1 (MQ=255)
                                        gCCATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGAc                                                                                                              >  2:917877/1‑150 (MQ=255)
                                          cATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACgg                                                                                                            >  1:599681/1‑150 (MQ=255)
                                              aTCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCtt                                                                                                         <  1:1245683/149‑1 (MQ=255)
                                              aTCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCtt                                                                                                         >  2:1245683/1‑149 (MQ=255)
                                               tCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGa                                                                                                       <  1:459173/150‑1 (MQ=255)
                                                                    gCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGca                                                                                  >  2:608852/1‑150 (MQ=255)
                                                                     cTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGcaa                                                                                 <  1:900811/150‑1 (MQ=255)
                                                                      ttGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGc                                                                                                                         <  1:1079017/109‑1 (MQ=255)
                                                                      ttGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGc                                                                                                                         >  2:1079017/1‑109 (MQ=255)
                                                                                 aGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGAc                                                                     >  1:456143/1‑150 (MQ=255)
                                                                                     aacaTCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAAcc                                                                            >  1:1222441/1‑139 (MQ=255)
                                                                                     aacaTCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAAcc                                                                            <  2:1222441/139‑1 (MQ=255)
                                                                                            aCTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATcc                                                          >  1:763822/1‑150 (MQ=255)
                                                                                               tCCACTTCTGTGCCGACCCGCACCAGGCTGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGa                                                       <  1:949573/150‑1 (MQ=255)
                                                                                                ccACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGat                                                      <  2:1164248/150‑1 (MQ=255)
                                                                                                  aCTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGatca                                                    >  1:386561/1‑150 (MQ=255)
                                                                                                  aCTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGatca                                                    >  2:827121/1‑150 (MQ=255)
                                                                                                       tgtgCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTc                                               <  2:386561/150‑1 (MQ=255)
                                                                                                           ccGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCACACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGatcatc                                                  <  1:295835/143‑1 (MQ=255)
                                                                                                           ccGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCACACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGatcatc                                                  <  1:295766/143‑1 (MQ=255)
                                                                                                           ccGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCACACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGatcatc                                                  >  2:295766/1‑143 (MQ=255)
                                                                                                           ccGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCACACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGatcatc                                                  >  2:295835/1‑143 (MQ=255)
                                                                                                              aCCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGa                                        >  2:221156/1‑150 (MQ=255)
                                                                                                              aCCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGa                                        <  1:917877/150‑1 (MQ=255)
                                                                                                                ccGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGAcc                                      >  2:1230631/1‑150 (MQ=255)
                                                                                                                         gCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAGGAAGAc                             <  1:221156/150‑1 (MQ=255)
                                                                                                                                     gCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCa                                                                                                                            <  2:1233414/43‑1 (MQ=255)
                                                                                                                                     gCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCa                                                                                                                            >  1:1233414/1‑43 (MQ=255)
                                                                                                                                       cATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAGGAAGAc                             <  2:761682/136‑1 (MQ=255)
                                                                                                                                       cATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAGGAAGAc                             >  1:761682/1‑136 (MQ=255)
                                                                                                                                       cATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAGGAAGACCCGCTGGTGAagag               >  1:47372/1‑150 (MQ=255)
                                                                                                                                        tcactcATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGAAATCACGATCATCGACCTGTCGACCAAGGAAGACCCGCTGGTGAAGAGt              >  1:44283/6‑150 (MQ=255)
                                                                                                                                         tGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAgg                                  >  2:113416/1‑129 (MQ=255)
                                                                                                                                         tGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAgg                                  <  1:113416/129‑1 (MQ=255)
                                                                                                                                           cgcATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAGGAAGACCCGCTGGTGAAGAGTGCg           <  1:1230631/150‑1 (MQ=255)
                                                                                                                                             cATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAGGAAGACCCGCTGGTGAAGAGTGCGGc         >  2:258208/1‑150 (MQ=255)
                                                                                                                                                  aGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAGGAAGACCCGCTGGTGaa                  <  2:264077/136‑1 (MQ=255)
                                                                                                                                                  aGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAGGAAGACCCGCTGGTGaa                  >  1:264077/1‑136 (MQ=255)
                                                                                                                                                   gCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAGGAAGACCCGCTGGTGAAGAGTGCGGCATTTgc   >  2:409256/1‑150 (MQ=255)
                                                                                                                                                    cccACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAGGAAGACCCGCTGGTGAAGAGTGCGGCATTTgcc  <  2:763822/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAACGAAAACTCGGCGATGGTGCTGCTGGTCGACGATCAGGCCATGATCGGCGAAGCGGTGCGGCGTGGCTTGGCCAATGAAGACAACATCGACTTCCACTTCTGTGCCGACCCGCACCAGGCGGTGGCCCAGGCCATGCGCATCAAGCCCACGGTCATTCTTCAAGACCTGATCATGCCTGGGCTGGACGGCTTGACCCTGGTGCGCGAGTACCGCAACAACCCGGTGACCCAGGACATCCCGATCATCGTCCTGTCGACCAAGGAAGACCCGCTGGTGAAGAGTGCGGCATTTGCC  >  NC_002947/1699395‑1699692

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