Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,508,074 T→A 11.4% Q48L (CAG→CTG)  PP_2201 ← alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,508,0740TA11.4% 124.6 / 7.3 44Q48L (CAG→CTG) PP_2201alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base A (2/3);  total (21/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.11e-01

GCGCGATGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCA  >  NC_002947/2507925‑2508212
                                                                                                                                                     |                                                                                                                                          
gcgcGATGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCt                                                                                                                                            <  2:251775/150‑1 (MQ=255)
gcgcGATGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCAATGTAGCCCATCTGCAGGGCCt                                                                                                                                            >  1:72604/1‑150 (MQ=255)
        aCCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCAGGCCCTgc                                                                                                                                    <  2:137754/150‑2 (MQ=255)
         cccGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTACCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCAGGCCCTgca                                                                                                                                   <  2:703108/150‑3 (MQ=255)
          ccGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTg                                                                                                                                     <  1:1268/147‑1 (MQ=255)
          ccGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTg                                                                                                                                     >  2:1268/1‑147 (MQ=255)
            ggTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGACAAATGCATGTAGCCCATCTGCAGGGCCTGCCCCTGGGGCt                                                                                                                                >  1:44881/1‑150 (MQ=255)
                     ccTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTGCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGa                                                                                                                       >  1:813848/1‑150 (MQ=255)
                     ccTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCATGGCCTGGCCCTGGGACTTGAAGTCGa                                                                                                                       >  2:234212/1‑150 (MQ=255)
                     ccTTGCTCAGGGCGTCGATAGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGCGCATCAGCACCCCGTTGTGGCCGTTGGGCTGGCCCTGGGCGGGCCCGTCCATGTGGCCCATCTGCAGGGCCTGGCCCTGCGACTTAACGTCGa                                                                                                                       >  2:78798/1‑150 (MQ=255)
                          cTCAGGGCGTCGATGTTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGc                                                                                                                  <  2:429036/150‑1 (MQ=255)
                          cTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAACGc                                                                                                                  <  1:851198/150‑1 (MQ=255)
                                     gATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGCCGAAATGCTTGAGGGGGTg                                                                                                       >  1:87149/1‑150 (MQ=255)
                                         gtggtTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCggg                                                                                                   >  2:497060/1‑150 (MQ=255)
                                                       ggcggcGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACtt                                                                                                                               >  1:435854/1‑108 (MQ=255)
                                                       ggcggcGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACtt                                                                                                                               <  2:435854/108‑1 (MQ=255)
                                                       ggcggcGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTc                                                                                     <  2:764981/150‑1 (MQ=255)
                                                              cAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTgg                                                                              <  1:395020/150‑1 (MQ=255)
                                                                gAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTgggg                                                                            <  2:840689/150‑1 (MQ=255)
                                                                    ttCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCg                                                                        >  2:596702/1‑150 (MQ=255)
                                                                     tCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCgt                                                                       <  1:38826/150‑1 (MQ=255)
                                                                       ttGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCgtag                                                                     <  1:508616/150‑1 (MQ=255)
                                                                                        cacGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCtt                                                                                      >  2:200212/1‑116 (MQ=255)
                                                                                        cacGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCtt                                                                                      <  1:200212/116‑1 (MQ=255)
                                                                                          cGTTGTGGCCGTGGGCCTGGCCTTGGGCGGGCAGGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGa                                                  <  1:78798/150‑1 (MQ=255)
                                                                                               tgGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGcc                                             >  1:638657/1‑150 (MQ=255)
                                                                                                  ccGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCATGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGaa                                          <  1:234212/150‑1 (MQ=255)
                                                                                                    gTTGGCCTGGCCTTGGGCGGGCAAGTTATAAAAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAAct                                        >  1:450730/1‑150 (MQ=255)
                                                                                                        gcctggcctTGGGCGGGAACGTCCATCTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAActgctg                                    <  2:412418/150‑1 (MQ=255)
                                                                                                          ctggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGc                                  >  2:743879/1‑150 (MQ=255)
                                                                                                          ctggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGc                                  >  1:690533/1‑150 (MQ=255)
                                                                                                           tggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGaa                                                                                           >  1:636738/1‑92 (MQ=255)
                                                                                                           tggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGaa                                                                                           <  2:636738/92‑1 (MQ=255)
                                                                                                             gcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGTTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAag                               <  2:72604/150‑1 (MQ=255)
                                                                                                               ctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTggggg                                                                           >  2:843725/1‑104 (MQ=255)
                                                                                                               ctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTggggg                                                                           <  1:843725/104‑1 (MQ=255)
                                                                                                                ttGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAaggca                            <  2:638657/150‑1 (MQ=255)
                                                                                                                            gTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGCTGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAg                >  1:79217/1‑150 (MQ=255)
                                                                                                                               cATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGc             <  2:690533/150‑1 (MQ=255)
                                                                                                                                  gTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAg                                                                                  <  2:50851/78‑1 (MQ=255)
                                                                                                                                  gTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAg                                                                                  >  1:50851/1‑78 (MQ=255)
                                                                                                                                         agtcccAGGGCCAGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgcgc   >  1:6305/6‑150 (MQ=255)
                                                                                                                                             ccAGGGCCAGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCgca  <  2:146509/146‑1 (MQ=255)
                                                                                                                                             ccAGGGCCAGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCgca  >  1:146509/2‑147 (MQ=255)
                                                                                                                                               aGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAActgct                                     >  2:643199/1‑110 (MQ=255)
                                                                                                                                               aGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAActgct                                     <  1:643199/110‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
GCGCGATGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCA  >  NC_002947/2507925‑2508212

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