Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,172,648 G→A 6.9% A94A (GCG→GCA)  PP_3671 → aldo/keto reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,172,6480GA 6.9% 169.6 / 2.5 58A94A (GCG→GCA) PP_3671aldo/keto reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base A (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

CCTACGAAAACGAAGCAGCCGTTGGCCAGGCGCTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGACGAAGGCAAGGCGCGCCACATTGGCGTGGCCAACTTCCCATTGGG  >  NC_002947/4172506‑4172788
                                                                                                                                              |                                                                                                                                            
ccTACGAAAACGAAGCAGCCGTTGGCCAGGCGCTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTggg                                                                                                                                         >  3:231281/1‑148 (MQ=255)
ccTACGAAAACGAAGCAGCCGTTGGCCAGGCGCTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCAggg                                                                                                                                         >  1:208442/1‑148 (MQ=255)
      aaaaCGAAGCAGCCGTTGGCCAGGCGCTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGa                                                                                                                                   >  3:57597/1‑148 (MQ=255)
              gcagcCGTTGGCCAGGCGCTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGaa                                                                                                                                  >  6:329289/1‑141 (MQ=255)
              gcagcCGTTGGCCAGGCGCTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGaa                                                                                                                                  <  5:329289/141‑1 (MQ=255)
                 gcCGTTGGCCAGGCGCTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCt                                                                                                                       >  2:205238/1‑149 (MQ=255)
                       ggCCAGGCGCTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCt                                                                                                                                            >  1:40719/1‑122 (MQ=255)
                       ggCCAGGCGCTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCt                                                                                                                                            <  2:40719/122‑1 (MQ=255)
                        gCCAGGCGCTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCca                                                                                                                 <  3:200126/148‑1 (MQ=255)
                         ccAGGCGCTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCcac                                                                                                                >  4:320384/1‑148 (MQ=255)
                                cTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTg                                                                                                        <  4:57597/149‑1 (MQ=255)
                                    gcgACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGc                                                                                                      >  1:114727/1‑147 (MQ=255)
                                     cgACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGTGTCAGTCGCTGGAAGACAGCCTGCGCGCACTGGGTACCGaa                                                                                                                                  >  2:288087/1‑118 (MQ=255)
                                     cgACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGTGTCAGTCGCTGGAAGACAGCCTGCGCGCACTGGGTACCGaa                                                                                                                                  <  1:288087/118‑1 (MQ=255)
                                     cgACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAg                                                                                                   <  2:114727/149‑1 (MQ=255)
                                      gACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAgg                                                                                                  >  8:87977/1‑149 (MQ=255)
                                        tgtcGATGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCa                                                                                                <  3:347658/146‑1 (MQ=255)
                                        cAGCGACGTGACACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCATGGAGACTAAAGACATGCGTCAGTAGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAgg                                                                                                  <  4:70818/147‑1 (MQ=255)
                                         aGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAggtcc                                                                                               >  7:77187/1‑146 (MQ=255)
                                             aCGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGAc                                                                                           <  1:205238/149‑1 (MQ=255)
                                                 gCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAgg                                                                                                                              >  5:151552/1‑110 (MQ=255)
                                                 gCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAgg                                                                                                                              <  6:151552/110‑1 (MQ=255)
                                                 gCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTgg                                                                                        >  4:82825/1‑148 (MQ=255)
                                                    aCGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGAcc                                                                                    >  4:86625/1‑149 (MQ=255)
                                                       cgAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTgg                                                                                 <  7:87977/149‑1 (MQ=255)
                                                        gAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTTGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGGCAGGCACCGACTGGGACCTGGc                                                                                <  6:66909/149‑1 (MQ=255)
                                                        gAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGc                                                                                >  3:58947/1‑149 (MQ=255)
                                                         aGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGcc                                                                               <  3:320384/149‑1 (MQ=255)
                                                          gCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGccc                                                                              >  8:208014/1‑149 (MQ=255)
                                                                  aCCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCACT‑GGTTCCCCACCGGTAGCACTGTTGCCCCTCCCCTGGCCAGGCCCCGGCTGGGCCCAGGCCCgccggatag                                                                     >  3:113544/1‑142 (MQ=255)
                                                                     tGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGa                                                                    >  2:14775/1‑148 (MQ=255)
                                                                         caccaAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACAccc                                                               <  3:86625/149‑1 (MQ=255)
                                                                           ccaAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTg                                                             <  7:208014/149‑1 (MQ=255)
                                                                            caAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTgg                                                            >  2:164385/1‑149 (MQ=255)
                                                                               ggTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAg                                                         >  8:100728/1‑149 (MQ=255)
                                                                                    gttggGACCGCCTGTAGCGTAATGCCATGCTTCAGTGGCTGTATGACTCCCTGCGCGCAC‑GGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGc                                                   <  4:113544/147‑1 (MQ=255)
                                                                                     gtggGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTgg                                                            <  1:216877/140‑1 (MQ=255)
                                                                                     gtggGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTgg                                                            >  2:216877/1‑140 (MQ=255)
                                                                                     gtggGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTAGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGc                                                   >  6:62870/1‑149 (MQ=255)
                                                                                     gtggGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTAGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGc                                                   <  5:62870/149‑1 (MQ=255)
                                                                                       gggACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCACCGACAcc                                                                <  2:191581/134‑1 (MQ=255)
                                                                                       gggACCGCCTGGAGCCTAAAGCCATGCATCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCACCGACAcc                                                                >  1:191581/1‑134 (MQ=255)
                                                                                         gACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGa                                               <  3:82825/149‑1 (MQ=255)
                                                                                          aCCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGAc                                              >  1:313520/1‑149 (MQ=255)
                                                                                             gCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGACGaa                                           <  7:195006/149‑1 (MQ=255)
                                                                                                       tAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACAccc                                                               <  4:19437/119‑1 (MQ=255)
                                                                                                       tAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGAc                                              >  5:94697/1‑136 (MQ=255)
                                                                                                       tAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGAc                                              <  6:94697/136‑1 (MQ=255)
                                                                                                                           tGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGt                                                           >  4:243220/1‑103 (MQ=255)
                                                                                                                           tGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGt                                                           <  3:243220/103‑1 (MQ=255)
                                                                                                                              aaGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGACGAAGGCAAGGCGCGCCACATTGGCGTGGCCAACTTc          <  1:137425/149‑1 (MQ=255)
                                                                                                                                     ccTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGACGAAGGCAAGGCGCGCCACATTGGCGTGGCCAACTTCCCATTgg   <  7:100728/149‑1 (MQ=255)
                                                                                                                                     ccTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGACGAAGGCAAGGCGCGCCACATTGGCGTGGCCAACTTCCCATTgg   <  6:241877/149‑1 (MQ=255)
                                                                                                                                     ccTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGACGAAGGCAAGGCGCGCCACATTGGCGTGGCCAACTTCCCATTgg   <  4:116628/149‑1 (MQ=255)
                                                                                                                                     ccTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGACGAAGGCAAGGCGCGCCACATTGGCGTGGCCAACTTCCCATTgg   >  3:116628/1‑149 (MQ=255)
                                                                                                                                      cTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGACGAAGGCAAGGCGCGCCACATTGGCGTGGCCAACTTCCCATTggg  <  1:14775/149‑1 (MQ=255)
                                                                                                                                       tGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGACGaaggcaaggc                                   <  7:93872/115‑1 (MQ=255)
                                                                                                                                       tGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGACGaaggcaaggc                                   >  8:93872/1‑115 (MQ=255)
                                                                                                                                              |                                                                                                                                            
CCTACGAAAACGAAGCAGCCGTTGGCCAGGCGCTGCGCGACAGCGACGTGCCACGCGAGCAGATCCACCTGACCACCAAGGTCTGGTGGGACCGCCTGGAGCCTAAAGCCATGCGTCAGTCGCTGGAAGACAGCCTGCGCGCGCTGGGTACCGAACAGGTCGACCTGTTCCACATCCACTGGCCAGGCACCGACTGGGACCTGGCCCGCAGCATCGACACCCTGGTAGCCCTGCGTGACGAAGGCAAGGCGCGCCACATTGGCGTGGCCAACTTCCCATTGGG  >  NC_002947/4172506‑4172788

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