Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,266,853 G→A 8.1% R30H (CGC→CAC)  PP_4642 → cyclic diguanosine monophosphate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,266,8530GA 8.1% 140.3 / 2.4 49R30H (CGC→CAC) PP_4642cyclic diguanosine monophosphate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (19/26);  new base A (2/2);  total (21/28)
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.86e-01

GGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAACGAGGTCCGTGTGCAGATGCAG  >  NC_002947/5266707‑5267001
                                                                                                                                                   |                                                                                                                                                    
ggCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGTTCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAATGCATCGACTTCGATGCCCCCACCGAGCTGCGcc                                                                                                                                                    <  1:47425/149‑1 (MQ=255)
ggCAAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCTGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTgcgc                                                                                                                                                     >  1:251026/1‑148 (MQ=255)
   cAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAgg                                                                                                                                                 <  1:99792/149‑1 (MQ=255)
      aTTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCTACTTCGATGCCCCCACCGAGCTGCGCCAGGGGc                                                                                                                                              <  8:10549/149‑1 (MQ=255)
          cAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGACTCGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTg                                                                                                                                            <  1:98166/148‑1 (MQ=255)
           aGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAggg                                                                                                                                                <  8:198247/142‑1 (MQ=255)
           aGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAggg                                                                                                                                                >  7:198247/1‑142 (MQ=255)
             tACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACg                                                                                                                                       >  6:174170/1‑149 (MQ=255)
                   aTGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGc                                                                                                                                  <  5:174170/148‑1 (MQ=255)
                   aTGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGc                                                                                                                                  <  2:251026/148‑1 (MQ=255)
                   aTGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCAATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGc                                                                                                                                  <  6:82842/148‑1 (MQ=255)
                          tGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCGGGCCTGCaaaa                                                                                                                          <  4:312469/149‑1 (MQ=255)
                           gCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAAc                                                                                                                         >  2:276209/1‑149 (MQ=255)
                               ctctAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAActgct                                                                                                                     <  1:276209/149‑1 (MQ=255)
                                    aCTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTgg                                                                                                                   <  5:136563/146‑1 (MQ=255)
                                    aCTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGAcc                                                                                                                <  6:266308/149‑1 (MQ=255)
                                        taGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCttt                                                                                                             <  6:239783/148‑1 (MQ=255)
                                             aGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAgg                                                                                                                                                 >  2:57140/1‑107 (MQ=255)
                                             aGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAgg                                                                                                                                                 <  1:57140/107‑1 (MQ=255)
                                                 gtcatccTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAgg                                                                                                   >  1:22722/6‑149 (MQ=255)
                                                              ccTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTgg                                                                                                                                   <  4:323579/104‑1 (MQ=255)
                                                              ccTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTgg                                                                                                                                   >  3:323579/1‑104 (MQ=255)
                                                                                       cGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCaaa                                                                                                     <  4:302849/109‑1 (MQ=255)
                                                                                       cGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCaaa                                                                                                     >  3:302849/1‑109 (MQ=255)
                                                                                       cGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGAcc                                                             <  3:216920/149‑1 (MQ=255)
                                                                                            aTGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCACCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACt                                                        >  7:200747/1‑149 (MQ=255)
                                                                                            aTGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCACCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACt                                                        >  7:195404/1‑149 (MQ=255)
                                                                                                cGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATc                                                                                         >  8:314145/1‑112 (MQ=255)
                                                                                                cGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATc                                                                                         <  7:314145/112‑1 (MQ=255)
                                                                                                     cgccgccgcTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGAttt                                                >  3:56477/1‑149 (MQ=255)
                                                                                                      gccgccgcTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTc                                               <  7:193782/149‑1 (MQ=255)
                                                                                                        cgccgcTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTc                                                                                                            >  2:58013/1‑86 (MQ=255)
                                                                                                        cgccgcTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGa                                             <  2:245708/149‑1 (MQ=255)
                                                                                                          ccgcTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACg                                           >  6:11174/1‑149 (MQ=255)
                                                                                                            gcTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGc                                          >  7:332456/1‑148 (MQ=255)
                                                                                                             cTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCca                                        >  7:139296/1‑149 (MQ=255)
                                                                                                              tttCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCca                                        >  2:93269/1‑148 (MQ=255)
                                                                                                                  aaCGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCgggc                                                             <  4:102274/123‑4 (MQ=255)
                                                                                                                  aaCGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCgggc                                                             >  3:102274/1‑120 (MQ=255)
                                                                                                                       aTCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCt                              >  1:352204/1‑149 (MQ=255)
                                                                                                                               cGATGCCCCCACCGAGCTGCACCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAAc                       <  8:195404/148‑1 (MQ=255)
                                                                                                                               cGATGCCCCCACCGAGCTGCACCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAAc                       <  8:200747/148‑1 (MQ=255)
                                                                                                                                gATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAAc                       >  7:232336/1‑147 (MQ=255)
                                                                                                                                gATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAACg                      <  4:56477/148‑1 (MQ=255)
                                                                                                                                 aTGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAACGAg                    <  2:22722/149‑1 (MQ=255)
                                                                                                                                  tGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAACGAgg                   <  5:11174/149‑1 (MQ=255)
                                                                                                                                     ccccACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGATTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAACGAGGTcc                >  7:44850/1‑149 (MQ=255)
                                                                                                                                     ccccACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAACGAGGTcc                >  2:199261/1‑149 (MQ=255)
                                                                                                                                     ccccACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAACGAGGTcc                <  8:232336/149‑1 (MQ=255)
                                                                                                                                              gCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAACGAGGTCCGTGTGCAGa       <  8:139296/149‑1 (MQ=255)
                                                                                                                                                   gcCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAACGAGGTCCGTGTGCAGATGCAg  <  2:352204/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCATGACGA‑ACGCCGCCGCTTTCAACGCATCGACTTCGATGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGACAACGAGGTCCGTGTGCAGATGCAG  >  NC_002947/5266707‑5267001

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