Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,852,397 T→A 9.6% Q312L (CAG→CTG)  PP_2505 ← GAF domain/GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,852,3970TA 9.6% 139.4 / 2.0 49Q312L (CAG→CTG) PP_2505GAF domain/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (25/19);  new base A (3/2);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACG  >  NC_002947/2852249‑2852540
                                                                                                                                                    |                                                                                                                                               
ccGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCt                                                                                                                                                 <  7:46101/149‑1 (MQ=255)
    aTGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGt                                                                                                                                             <  8:134163/149‑1 (MQ=255)
    aTGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGt                                                                                                                                             <  3:141206/149‑1 (MQ=255)
     tgtgATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTg                                                                                                                                            <  3:11770/149‑1 (MQ=255)
            cGTACCCTGCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCAGGGTGAGCAga                                                                                                                                      <  8:119329/148‑2 (MQ=255)
             gTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATTACGCTGGCACATAGTCAACTTGGTCACGTGCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGgc                                                                                                                                    <  8:93151/149‑1 (MQ=255)
                 ccTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCAc                                                                                                                                >  5:153396/1‑149 (MQ=255)
                   tCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCAcc                                                                                                                               >  4:292371/1‑148 (MQ=255)
                            cGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCAt                                                                                                                     <  6:153396/149‑1 (MQ=255)
                              gagaTCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGCGCCCTGCTCACCCa                                                                                                                             >  6:163722/1‑139 (MQ=255)
                              gagaTCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGCGCCCTGCTCACCCa                                                                                                                             <  5:163722/139‑1 (MQ=255)
                               agaTCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGt                                                                                                                  >  4:229878/1‑149 (MQ=255)
                                   cTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTc                                                                                                                 <  7:113689/146‑1 (MQ=255)
                                   cTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTc                                                                                                                 >  8:113689/1‑146 (MQ=255)
                                   cTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACAATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATc                                                                                                              >  7:122474/1‑149 (MQ=255)
                                    ttAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAACCTGCTCACCCAGGGTGAGCAgat                                                                                                                                     >  7:274554/1‑123 (MQ=255)
                                    ttAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAACCTGCTCACCCAGGGTGAGCAgat                                                                                                                                     <  8:274554/125‑3 (MQ=255)
                                         ttCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGCGCCCTGCTCACCCagg                                                                                                                           <  3:9155/130‑3 (MQ=255)
                                         ttCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACAGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGCGCCCTGCTCACCCagg                                                                                                                           >  4:9155/1‑128 (MQ=255)
                                           ccTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCa                                                                                                                                 <  8:142449/122‑1 (MQ=255)
                                           ccTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCa                                                                                                                                 >  7:142449/1‑122 (MQ=255)
                                           ccTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGa                                                                                                        <  3:292371/147‑1 (MQ=255)
                                                     aCGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGAATAGAGACCGCCGTCAGAGAGCTTTATTGCATCATAAATCTGCTCACCCTGCGTGAGCATTGCTCACACAATGCTCATGGTCATCTTCCGATGGCCGCCggt                                                                                             >  8:332402/1‑147 (MQ=255)
                                                                   aCAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGgcgc                                                                              >  4:230690/1‑149 (MQ=255)
                                                                     aaCTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGgcgcgc                                                                            >  8:146064/1‑149 (MQ=255)
                                                                     aaCTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTAATGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACCGCACCGgcgcgc                                                                            >  8:290894/1‑149 (MQ=255)
                                                                          ggTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGc                                                                         <  3:229878/147‑1 (MQ=255)
                                                                               cGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCCCCCTGGGGGAGCAGGGCGCACCCAATCCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATtga                                                                  <  4:293300/149‑1 (MQ=255)
                                                                                  ccaccaTGTTTGGCCGCATAGAGACAGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGTAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATtgact                                                                >  6:359336/1‑148 (MQ=255)
                                                                                                 gCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCg                                                 >  4:156711/1‑148 (MQ=255)
                                                                                                  cATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGAGCGCTGCCGATTGACTGACGCAAGCGCTCCGc                                                >  7:339486/1‑148 (MQ=255)
                                                                                                      gagaCCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCtt                                                                                                            <  8:45001/84‑1 (MQ=255)
                                                                                                      gagaCCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCtt                                                                                                            >  7:45001/1‑84 (MQ=255)
                                                                                                        gaCCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTgc                                         <  3:230690/149‑1 (MQ=255)
                                                                                                         aCCGCCGTCAGCGCGATTTATTGCAGGATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAAGAGTAACTGCGCGCTGCCGATTGACTGACGCAAGCGCTTCGCAATCTgc                                         <  8:122474/148‑1 (MQ=255)
                                                                                                          ccgccgTCAGCGCGCTTTATTTCATCATCAATCTGCTCACCCTGGGTGAGCAGTGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTgcg                                        <  7:332402/148‑1 (MQ=255)
                                                                                                             ccgTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCaa                                    >  3:71831/1‑149 (MQ=255)
                                                                                                                     gcgcTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGc                            >  4:295387/1‑149 (MQ=255)
                                                                                                                            aTTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCtt                      <  3:156711/148‑1 (MQ=255)
                                                                                                                                 atcatcAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAAccc                >  2:326178/1‑149 (MQ=255)
                                                                                                                                   catcaATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTgc                                         <  4:33670/122‑1 (MQ=255)
                                                                                                                                   catcaATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTgc                                         >  3:33670/1‑122 (MQ=255)
                                                                                                                                   catcaATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCg               >  1:64504/1‑148 (MQ=255)
                                                                                                                                    atcaATCTGCTCACCCTGGGTGAGCAGGGCGCTCCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCg               >  8:155221/1‑147 (MQ=255)
                                                                                                                                       aaTCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGTCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGc           >  7:82535/1‑148 (MQ=255)
                                                                                                                                        cccTGCTCACCCAGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGc           >  6:88184/3‑147 (MQ=255)
                                                                                                                                          cTGCTCACCCAGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACTCGTAGccttg       >  8:295912/1‑145 (MQ=255)
                                                                                                                                            gCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGc      >  6:119061/1‑148 (MQ=255)
                                                                                                                                                aCCCTGGGTGAGTAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGTCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTAc   <  5:359336/147‑1 (MQ=255)
                                                                                                                                                aCCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACg  >  1:168124/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACG  >  NC_002947/2852249‑2852540

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: ‑