Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,141,558 T→A 8.4% V174E (GTA→GAA)  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,5580TA 8.4% 158.7 / 2.1 55V174E (GTA→GAA) PP_3645L‑proline amide hydrolase
Reads supporting (aligned to +/- strand):  ref base T (22/28);  new base A (2/3);  total (24/31)
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

CCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCAT  >  NC_002947/4141411‑4141705
                                                                                                                                                    |                                                                                                                                                    
cccccGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTa                                                                                                                                                     <  4:171661/149‑1 (MQ=255)
 ccccGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTAc                                                                                                                                                    >  3:58068/1‑149 (MQ=255)
 ccccGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTAc                                                                                                                                                    >  1:309742/1‑149 (MQ=255)
   ccGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTg                                                                                                                                                  <  4:83908/149‑1 (MQ=255)
     gAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCAGTGAATACAAGTCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGc                                                                                                                                                >  5:263454/1‑149 (MQ=255)
       gCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGTACCTGGCGa                                                                                                                                              <  4:14728/149‑1 (MQ=255)
          aaGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGTACCTGGCGaatg                                                                                                                                           <  3:11575/149‑4 (MQ=255)
           aGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCaa                                                                                                                                          <  7:137344/149‑1 (MQ=255)
               ccGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGAc                                                                                                                                      >  4:268811/1‑149 (MQ=255)
                cGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGAcc                                                                                                                                     <  6:96798/149‑1 (MQ=255)
                 ggggTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGAccc                                                                                                                                    >  8:275008/1‑149 (MQ=255)
                    gTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCg                                                                                                                                   >  7:166403/1‑147 (MQ=255)
                      gATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAgg                                                                                                                               <  3:268811/149‑1 (MQ=255)
                       aTCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCTGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAggg                                                                                                                              <  6:42019/149‑1 (MQ=255)
                        tCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAGCGACCCGCAGGGc                                                                                                                             >  2:175490/1‑149 (MQ=255)
                          gCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑t                                                                                                                           <  2:222265/148‑1 (MQ=255)
                          gCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑t                                                                                                                           <  2:290442/148‑1 (MQ=255)
                            ccGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑Tac                                                                                                                         >  1:159804/1‑148 (MQ=255)
                                gcTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGGCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑Taccacac                                                                                                                    <  1:226206/149‑1 (MQ=255)
                                    cccAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑Taccacacca                                                                                                                  <  7:188938/147‑1 (MQ=255)
                                     ccACGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGTACCTGGCGaatg                                                                                                                                           <  6:69361/122‑4 (MQ=255)
                                     ccAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGTACCTGGCGaatg                                                                                                                                           >  5:69361/1‑119 (MQ=255)
                                      cAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTa                                                                                                              >  4:41183/1‑149 (MQ=255)
                                      cAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTa                                                                                                              >  2:257519/1‑149 (MQ=255)
                                       aaGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTAt                                                                                                             <  2:187037/149‑1 (MQ=255)
                                            gcggAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCt                                                                                                        <  7:61824/149‑1 (MQ=255)
                                              ggAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGt                                                                                                      <  7:275008/149‑1 (MQ=255)
                                                  gtGGGCCCAGACGCCAACGTCGATGCTGCACTGGACCGCTGGTTCAGCCGTGAATACAAGGATGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGc                                                                                                  <  6:263454/148‑1 (MQ=255)
                                                                cAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATg                                                                                     >  1:239386/1‑148 (MQ=255)
                                                                 aaCGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTa                                                                                   <  7:198332/149‑1 (MQ=255)
                                                                    gTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACcg                                                                                >  7:189954/1‑149 (MQ=255)
                                                                       gATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCgcgc                                                                             <  4:58068/149‑1 (MQ=255)
                                                                       gATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCgcgc                                                                             <  8:189954/149‑1 (MQ=255)
                                                                             gcCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCcgacga                                                                       >  5:275061/1‑149 (MQ=255)
                                                                                 cggACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTg                                                                   <  8:126443/148‑1 (MQ=255)
                                                                                   gACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTgg                                                                  <  6:275061/148‑1 (MQ=255)
                                                                                     ccGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCa                                                               <  2:159804/149‑1 (MQ=255)
                                                                                           gTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGCATACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATc                                                          <  1:257519/149‑1 (MQ=255)
                                                                                            ttCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TTCCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCa                                                        <  8:91942/149‑1 (MQ=255)
                                                                                             taaGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCTAGGAACTGGCGACCAACGACCCGCAGGGA‑AACCACACCAATAATTCCATGTTCGCCACTCAGGACATGTACCCCCCCGACGCCCTGGGCAGTaac                                                          >  1:80550/3‑145 (MQ=255)
                                                                                               aGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGt                                                     <  8:166403/149‑1 (MQ=255)
                                                                                                  cGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGcc                                                  >  1:315240/1‑149 (MQ=255)
                                                                                                       aTACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGc                                             >  5:120401/1‑149 (MQ=255)
                                                                                                             ggCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAGCGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCg                                       <  1:175490/149‑1 (MQ=255)
                                                                                                                 gCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCAc                                   >  8:234983/1‑149 (MQ=255)
                                                                                                                   cAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCg                                 >  3:223086/1‑149 (MQ=255)
                                                                                                                    aaCCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCgg                                >  5:151422/1‑149 (MQ=255)
                                                                                                                    aaCCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCgg                                >  3:36346/1‑149 (MQ=255)
                                                                                                                     aCCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGc                               >  8:26916/1‑149 (MQ=255)
                                                                                                                       cccGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGa                             <  3:41183/149‑1 (MQ=255)
                                                                                                                                    ccgccATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGc                >  8:249486/1‑149 (MQ=255)
                                                                                                                                      gccATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTc              <  6:120401/149‑1 (MQ=255)
                                                                                                                                               c‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGc     <  4:223086/149‑1 (MQ=255)
                                                                                                                                                  ggTACTGGCGAGCAACGACCCGCAGGGC‑TATCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCa   <  2:209521/149‑1 (MQ=255)
                                                                                                                                                   gTACTGGCGAGCAACGACCCGGAGGGC‑TACCACACCACTTATTCGCTGTTCTCCACTCTGGACATGTACAGCGCCGACGACCTGGGCAGTATACAGGTGCCGACGCTGATCGCCACCGGCGATCTCGACTCGGCCTCCACACCGGCCAt  >  5:14657/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCC‑AGGTACTGGCGAGCAACGACCCGCAGGGC‑TACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCAT  >  NC_002947/4141411‑4141705

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