Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,707,463 G→A 8.5% A267V (GCG→GTG)  PP_4164 ← alpha/beta hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,707,4630GA 8.5% 131.6 / 2.4 47A267V (GCG→GTG) PP_4164alpha/beta hydrolase
Reads supporting (aligned to +/- strand):  ref base G (23/20);  new base A (2/2);  total (25/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 = 6.73e-01

CCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCGGGCAACTCACCCGCCGCACATAGGC  >  NC_002947/4707323‑4707605
                                                                                                                                            |                                                                                                                                              
ccGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTc                                                                                                                                        <  5:139676/149‑1 (MQ=255)
    aTTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCACATCGTTcaggt                                                                                                                                    >  2:50252/1‑149 (MQ=255)
         aGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCGAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTcaggtgcagg                                                                                                                               <  1:147352/149‑1 (MQ=255)
         aGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTcaggtgcagg                                                                                                                               <  4:304735/149‑1 (MQ=255)
             ccAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCg                                                                                                                                           <  5:242786/133‑1 (MQ=255)
             ccAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCg                                                                                                                                           >  6:242786/1‑133 (MQ=255)
                         gCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCggg                                                                                                                 >  6:175000/1‑147 (MQ=255)
                         gCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCggg                                                                                                                 >  2:221177/1‑147 (MQ=255)
                                   ggTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTg                                                                                                     >  2:273069/1‑149 (MQ=255)
                                     ttGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGt                                                                                                    >  1:359288/1‑148 (MQ=255)
                                          tcCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAg                                                                                              >  2:159307/1‑149 (MQ=255)
                                           cccTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAgg                                                                                             >  7:166679/1‑149 (MQ=255)
                                            ccTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAg                                                                                              <  3:22762/147‑1 (MQ=255)
                                            ccTGTTCAACCCGCCAAGCCTCGACACTAGTCCGGTGCAAGTCAAGCAAAGCAAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGAGATGGCCGCCGGGCAGCGTCACCTGTTCAAAggg                                                                                            >  1:201686/1‑149 (MQ=255)
                                                ttCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAg                                                                                         <  1:159307/148‑1 (MQ=255)
                                                   aaCCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAgctgc                                                                                     <  5:175000/149‑1 (MQ=255)
                                                    aCCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAgctgct                                                                                    >  3:187984/1‑149 (MQ=255)
                                                     cccGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTc                                                                                   <  8:166679/149‑1 (MQ=255)
                                                      ccGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTcc                                                                                  <  2:359288/149‑1 (MQ=255)
                                                        gCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAgg                                                                                             <  5:286035/136‑1 (MQ=255)
                                                        gCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAgg                                                                                             >  6:286035/1‑136 (MQ=255)
                                                             gCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAAcc                                                                           >  7:214615/1‑149 (MQ=255)
                                                               cTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCg                                                                         <  4:120704/149‑1 (MQ=255)
                                                                tCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCgg                                                                        >  1:297861/1‑149 (MQ=255)
                                                                       ttGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCACATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGgcg                                                                 >  6:9648/1‑149 (MQ=255)
                                                                              gTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAg                                                                                         >  5:66099/1‑118 (MQ=255)
                                                                              gTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAg                                                                                         <  6:66099/118‑1 (MQ=255)
                                                                                       aaGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCACATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTc                                                  <  5:9648/148‑1 (MQ=255)
                                                                                         gCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGCCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGcc                                                     <  5:149071/143‑1 (MQ=255)
                                                                                         gCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGCCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGcc                                                     >  6:149071/1‑143 (MQ=255)
                                                                                                 aaaGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGccacca                                       <  1:273069/149‑1 (MQ=255)
                                                                                                 aaaGAAGCGAATGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCTGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGccacca                                       <  5:259379/149‑1 (MQ=255)
                                                                                                          aTTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAgcg                              >  1:302471/1‑149 (MQ=255)
                                                                                                           ttGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGtt                                                                                                   >  6:61969/1‑79 (MQ=255)
                                                                                                           ttGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGtt                                                                                                   <  5:61969/79‑1 (MQ=255)
                                                                                                             gATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGcc                            <  4:187984/148‑1 (MQ=255)
                                                                                                              aTACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCgg                          >  1:221224/1‑149 (MQ=255)
                                                                                                              aTACAGTCTGCAACAAGGGCTGCGCCCTGCACATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCg                           <  1:50252/148‑1 (MQ=255)
                                                                                                                      tGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCGGGCAACTCa                  >  5:287069/1‑149 (MQ=255)
                                                                                                                        caacaaGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCGGGCAACTCAcc                >  5:153978/1‑149 (MQ=255)
                                                                                                                          acaaGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCGGGCAACTCACccg              >  2:128127/1‑149 (MQ=255)
                                                                                                                            aaGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCACCCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCGGGCAACTCACccgcc            >  1:164169/1‑149 (MQ=255)
                                                                                                                             aGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCGGGCAACTCACccgccg           <  2:221224/149‑1 (MQ=255)
                                                                                                                              gggCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCGGGCAACTCACCcgccgc          <  8:214615/149‑1 (MQ=255)
                                                                                                                                gCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCGGGCAACTCACCCGCCGcac        <  1:128127/149‑1 (MQ=255)
                                                                                                                                  tGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCGGGCAACTCACCCGCCGcaca       >  3:264685/1‑148 (MQ=255)
                                                                                                                                       ccTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCGGGCAACTCACCCGCCGCACATAGGc  >  6:123303/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                              
CCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATCGTTCAGGTGCAGGTGATGGCCGCCGGGCAGCGTCACCTGTTCAAAGGGTAGCTGCTCCAGCAACCCGGTGTGGCGCGCCAGCATGCCGTCGGCGGCCACCACCAAGAGCGCCGGGCAACTCACCCGCCGCACATAGGC  >  NC_002947/4707323‑4707605

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