Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,213,182 A→G 10.7% Q110Q (CAA→CAG patA → putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,213,1820AG10.7% 92.4 / 2.0 35Q110Q (CAA→CAGpatAputrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (14/16);  new base G (2/2);  total (16/19)
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.92e-01

GAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAGGGCCAGGAGTTTATCGACTGCCTGGGAGGTTTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTAT  >  CP009273/3212918‑3213454
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                
gaAGCGAACGCTGGATCATGAGGAGATGAAAGTACTTATCCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGAGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAGGGCCAGGAGTTTATCGACTGCCTGGGAGGTTTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAaac                                                                                                                                                                                                                                                                        <  1:768391/275‑1 (MQ=255)
          cTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGTGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAGGGCCAGGAGTTTAGCGACTGCCTGGGAGGTTTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGca                                                                                                                                                                                                                                                              <  1:59575/275‑1 (MQ=255)
          cTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCAGGAAGTTTAAATACGCTTGTCGACACCCAGGGCCAGGAGTTTATCGACTGCCTGGGAGGTTTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAAGAACCGCTGca                                                                                                                                                                                                                                                              <  1:798344/275‑1 (MQ=255)
                                                                    aTGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAGGGCCAGGAGTTTCTCGACTGCCTGGGAGGTTTTGGAATTTTCAACGTGGGGCCCCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACTGGCGATGTTTGCGAAAACCCTTTCTTCGCt                                                                                                                                                                                                     >  1:301954/1‑274 (MQ=255)
                                                                    aTGTCAATCCGGGGTTTTTAGAGTATAGCAAATCGGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAGGGCCAGGAGTTTATCGACTGCCTGGGAGGTTTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGNCGTACAGANTCAACTTGCGAAACAACCGc                                                                                                                                                                                                                                                                  >  1:550289/1‑213 (MQ=255)
                                                                                                                                   tAGAGTGGCAAGCGGGAAGATTAAATACGCTTGTCGACACCCAGGGGCAGGAGTTCACCGACTGCCTGGGAGGTTGTGGAATTTTCAACGTGGGGAACCGTAATCCAGTTGTGGTTTCCGCCGTACAGCATCAACATGCGGAACAACCGCTGcac                                                                                                                                                                                                                                                             >  2:198088/1‑155 (MQ=255)
                                                                                                                                       gTGGCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAGGGCCAGGAGTTTATCGACTGCCTGGGAGGTTTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGcc                                                                                                                                                                                                                                                         <  1:198091/155‑1 (MQ=255)
                                                                                                                                          gCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAGGGCCAGGAGTTTATCGACTGCCTGGGAGGTTTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCAATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGAGGACCCGAGTCCGTTGAAGACGCACTGAAg                                                                                                                               >  1:476273/1‑274 (MQ=255)
                                                                                                                                                 ggAAGTTTGAATACGCTTGTCGACACCCAGGGCCAGGAGGTTATCGACTGCCTGGTAGGTTTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGc                                                                                                                                                                                                                                                               <  2:316023/139‑1 (MQ=255)
                                                                                                                                                 ggAAGTTTAAATACGCTTGTCGACACCCAGGGCCAGGAGTTTATCGTCTGCCGGGGAGGTTTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGa                                                                                                                        <  1:345947/274‑1 (MQ=255)
                                                                                                                                                         aaaTACGCTTGTCGACACCCAGGGTCAGGAGTTTATCGACTGCCTGGGAGGTTTTGGTATTTTCAACGTGGGGCACCGTAATGCAGTTGTGGTTTACGCCGTACAGAATCAAATTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCACGGTAAACTGAAATACAGCTTCGTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTAc                                                                                                                <  1:274033/274‑1 (MQ=255)
                                                                                                                                                                                    ggAGTTTATCGACTGCCTGGGAGGTTTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTa                                                                                    >  1:350642/1‑275 (MQ=255)
                                                                                                                                                                                                             tttGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGc                                                                                                                                                                                                                                                                  >  2:455308/1‑76 (MQ=255)
                                                                                                                                                                                                               tGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAAcc                                                                                                                                                                                                                >  2:15714/1‑124 (MQ=255)
                                                                                                                                                                                                                 gAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGca                                                                                                                                                                                                                                                              <  1:455311/76‑1 (MQ=255)
                                                                                                                                                                                                                   aTTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTtg                                                                                                                                                                                                            <  1:15714/124‑1 (MQ=255)
                                                                                                                                                                                                                      ttCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAGCTTGCGAAACAACCGc                                                                                                                                                                                                                                                                  >  2:337234/1‑67 (MQ=255)
                                                                                                                                                                                                                          aCGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGca                                                                                                                                                                                                                                                              <  1:337237/67‑1 (MQ=255)
                                                                                                                                                                                                                           cggggggCACCGTAATCCAGTTGTGGTTTGCGCCGGACAGAATCAGCATGCTCAACATCCGCTGCACAGCCAGGAGCGTCTAGTTCCGTTACGGGCGATGTTGGAGAAAGCCCGTGCTGCGCTAACGCCCGGTAAACTGGAATACAGCTTGTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAGATCACTTGgcg                                             <  2:350639/272‑1 (MQ=255)
                                                                                                                                                                                                                           ccgggggCACCGTCATCCAGTTGTGCTTTCGGCCGTACAGAATCAACTTGCGAAACAACCGAGGCCCAGCCAGGTGCTGTTCGGTCCGTTACGGTCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATCCAGCTTCTTCTGTAATAGCGGCACCGAGTCCGGTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGgcg                                             <  2:803410/272‑1 (MQ=255)
                                                                                                                                                                                                                                      gTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTa                                                                                                                                                                                                                                     >  1:434548/1‑80 (MQ=255)
                                                                                                                                                                                                                                       tAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGcac                                                                                                                                                                                                                                                             >  1:604188/1‑55 (MQ=255)
                                                                                                                                                                                                                                       tAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGt                                                                                                                                                   >  1:639171/1‑161 (MQ=255)
                                                                                                                                                                                                                                          tCCAGTTGTGGTTTCCGCCGTACAGAATCAGCTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTgg                                                                                                                                                                                                                        >  1:763015/1‑89 (MQ=255)
                                                                                                                                                                                                                                          tCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACggg                                                                                                                                                                                                                                 <  2:434545/80‑1 (MQ=255)
                                                                                                                                                                                                                                           ccAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGcc                                                                                                                                                                                                                                                         <  2:604182/55‑1 (MQ=255)
                                                                                                                                                                                                                                           ccAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGt                                                                                                                                               <  2:639165/161‑1 (MQ=255)
                                                                                                                                                                                                                                              gTTGTGGTTTCCGCCGTACAGAATCAGCTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGaa                                                                                                                                                                                                                    <  2:763009/89‑1 (MQ=255)
                                                                                                                                                                                                                                                    gTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAgctgc                                                                                                                                                                                                                                                >  1:246889/1‑55 (MQ=255)
                                                                                                                                                                                                                                                    gTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGc                                                                                                                                                             >  2:792987/1‑138 (MQ=255)
                                                                                                                                                                                                                                                        ccgccgTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGCGCTGCTCGa                                                                                                                                                                                                                                            <  2:246886/55‑1 (MQ=255)
                                                                                                                                                                                                                                                        ccgccgTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCa                                                                                                                                                         <  1:792993/138‑1 (MQ=255)
                                                                                                                                                                                                                                                            cgTAAAAAATCACCTGGCCAAACACCCGCTGCACAGCCAGGTGCTCCGCGTTCTGTTTCGGGCTATGTTGCCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATGCAGCTTCTTCTGTAATAGCGGCACCTAGGCCGTTGAATCAGCGCGTAAGCTTGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACGTTTATTGCCAGCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCa            <  2:651328/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGccc                                                                                                                                                                                              >  2:239697/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                      cAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTCCCAGTCACCGCGCGGCAAGNTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTAt  >  1:651334/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                        aCTTGCGAAACAACCGGTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGAGGTTGGCGAAAACCCTTGCTGCGCTAACGcgc                                                                                                                                                                                              >  1:717396/1‑83 (MQ=255)
                                                                                                                                                                                                                                                                        aCTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAg                                                                                                                                                                             >  2:455999/1‑102 (MQ=255)
                                                                                                                                                                                                                                                                        aCTTGCGAAACAACCGCTGCACAGCCAGCCGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGCGTCCgct                                                                                                                                              >  1:649077/1‑131 (MQ=255)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                
GAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAGGGCCAGGAGTTTATCGACTGCCTGGGAGGTTTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTAT  >  CP009273/3212918‑3213454

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: