Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,141,168 T→A 8.0% L44Q (CTG→CAG)  PP_3645 → L‑proline amide hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,141,1680TA 8.0% 139.9 / 3.3 50L44Q (CTG→CAG) PP_3645L‑proline amide hydrolase
Reads supporting (aligned to +/- strand):  ref base T (20/26);  new base A (2/2);  total (22/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 = 6.06e-01

GGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCG  >  NC_002947/4141026‑4141308
                                                                                                                                              |                                                                                                                                            
ggAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCaa                                                                                                                                        <  6:243249/149‑1 (MQ=255)
    cgcCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGAc                                                                                                                                    <  1:135987/149‑1 (MQ=255)
            aTGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGt                                                                                                                            >  1:305731/1‑149 (MQ=255)
             tGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTc                                                                                                                           <  3:299772/149‑1 (MQ=255)
                 tCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAG‑TCGTTGGCCAGGCCAACGaactga                                                                                                                                <  4:310893/139‑6 (MQ=255)
                 tCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAG‑TCGTTGGCCAGGCCAACGaactga                                                                                                                                >  3:310893/1‑134 (MQ=255)
                     ccTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTa                                                                                                                   <  4:113496/149‑1 (MQ=255)
                           gCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACAt                                                                                                             <  6:344572/149‑1 (MQ=255)
                           gCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACAt                                                                                                             >  6:139621/1‑149 (MQ=255)
                                           ccGGCAGCTGCTTTCTGGATGTCTGCGAGGGCCAACCCGTGGTGCTGATCCACGGCGTTGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTc                                                                                               <  3:63868/147‑1 (MQ=255)
                                                 ccAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc                                                                                       <  4:93624/149‑1 (MQ=255)
                                                  cAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgcg                                                                                      >  5:273438/1‑149 (MQ=255)
                                                      tATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTAcc                                                                                  >  2:255100/1‑149 (MQ=255)
                                                         cTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGgc                                                                               >  1:78511/1‑149 (MQ=255)
                                                                  ggCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGacac                                                                      >  3:12643/1‑149 (MQ=255)
                                                                  ggCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGaca                                                                       >  1:313317/1‑148 (MQ=255)
                                                                  ggCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGaca                                                                       >  8:206173/1‑148 (MQ=255)
                                                                   gCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCCCGTACCGGCAGCCGacaca                                                                     >  7:47895/1‑149 (MQ=255)
                                                                        ggCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCaaa                                                                                            <  6:8203/121‑1 (MQ=255)
                                                                        ggCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCaaa                                                                                            >  5:8203/1‑121 (MQ=255)
                                                                             aCCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAg                                                           >  3:218557/1‑149 (MQ=255)
                                                                                 gtggtgCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTa                                                       >  7:321314/1‑149 (MQ=255)
                                                                                    gtgCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATg                                                     <  2:313317/148‑1 (MQ=255)
                                                                                    gtgCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGc                                                    <  2:305731/149‑1 (MQ=255)
                                                                                          aTCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGat                                                <  4:218557/147‑2 (MQ=255)
                                                                                          aTCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGAc                                                >  6:162712/1‑147 (MQ=255)
                                                                                          aTCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGAc                                                <  1:255100/147‑1 (MQ=255)
                                                                                          aTCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGAc                                                <  5:139621/147‑1 (MQ=255)
                                                                                            ccACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGc                                            >  5:232841/1‑149 (MQ=255)
                                                                                             cACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCt                                           >  5:344287/1‑149 (MQ=255)
                                                                                                       gCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTg                                 >  6:312365/1‑149 (MQ=255)
                                                                                                         cTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCt                               >  2:105765/1‑149 (MQ=255)
                                                                                                            aaCAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCtt                              <  4:12643/147‑1 (MQ=255)
                                                                                                            aaCAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTg                             <  3:136268/148‑1 (MQ=255)
                                                                                                            aaCAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGa                            <  5:312365/149‑1 (MQ=255)
                                                                                                            aaCAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCGCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGa                            <  7:238234/149‑1 (MQ=255)
                                                                                                              cAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGacc                          <  1:105765/149‑1 (MQ=255)
                                                                                                                aGGAAATGTGGGGAGGTCAGTTCGTTGGCCTGGCAAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACTGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCGGACCAGCTCGCCGAACTGCTTGaccac                        <  8:165546/149‑1 (MQ=255)
                                                                                                                  gAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCt                      <  6:273438/149‑1 (MQ=255)
                                                                                                                           ggCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGacc                          >  5:187741/1‑136 (MQ=255)
                                                                                                                           ggCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGacc                          <  6:187741/136‑1 (MQ=255)
                                                                                                                           ggCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGc             <  6:232841/149‑1 (MQ=255)
                                                                                                                               gTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTTTGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAATTTGcaggt         >  1:280903/1‑146 (MQ=255)
                                                                                                                               gTCAGTGCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACaa         <  8:292289/149‑1 (MQ=255)
                                                                                                                                tCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAg        <  2:78511/149‑1 (MQ=255)
                                                                                                                                  aGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGc       >  8:27907/1‑148 (MQ=255)
                                                                                                                                      cGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGCCCACGGGCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCAcc   <  8:321314/148‑1 (MQ=255)
                                                                                                                                       gTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCg  <  7:117324/148‑1 (MQ=255)
                                                                                                                                       gTTGGCCAGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcg                                                                                        <  4:174159/62‑1 (MQ=255)
                                                                                                                                       gTTGGCCAGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcg                                                                                        >  3:174159/1‑62 (MQ=255)
                                                                                                                                           gcctggccAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCtt                              >  1:327116/1‑116 (MQ=255)
                                                                                                                                           gcctggccAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCtt                              <  2:327116/116‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCG  >  NC_002947/4141026‑4141308

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