Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,956 G→A 6.4% P31S (CCG→TCG)  allE ← S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913543,9560GA 6.4% 200.2 / 4.4 62P31S (CCG→TCG) allES‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base A (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

ACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCAAGTTGATAAAGCATGAGTGCCAACGTTTTGA  >  NC_000913/543816‑544103
                                                                                                                                            |                                                                                                                                                   
aCCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTAACAGACCATCCGGGGTTAAc                                                                                                                                              >  4:74589/1‑148 (MQ=255)
  cGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAg                                                                                                                                            <  4:113693/148‑1 (MQ=255)
   gAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTAACAGACCATCCGGGGTTAACAGt                                                                                                                                           <  3:74589/148‑1 (MQ=255)
           gttggttgCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCgg                                                                                                                                                     <  1:103732/130‑1 (MQ=255)
           gttggttgCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTGTTTACCAGACCATCCgg                                                                                                                                                     >  2:103732/1‑130 (MQ=255)
            ttggttgCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGaa                                                                                                                                      >  3:1105517/1‑144 (MQ=255)
            ttggttgCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGaa                                                                                                                                      <  4:1105517/144‑1 (MQ=255)
                     aCCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTa                                                                                                                         >  3:691593/1‑148 (MQ=255)
                     aCCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTa                                                                                                                         >  2:588044/1‑148 (MQ=255)
                      ccGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTa                                                                                                                                  <  3:664076/138‑1 (MQ=255)
                      ccGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGCCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTa                                                                                                                                  >  4:664076/1‑138 (MQ=255)
                        gTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACa                                                                                                                      <  3:3094/148‑1 (MQ=255)
                         ttttGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACaa                                                                                                                     <  4:691593/148‑1 (MQ=255)
                         ttttGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACaa                                                                                                                     <  1:588044/148‑1 (MQ=255)
                                    tgtgACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTAc                                                                                                                                 <  2:619514/125‑1 (MQ=255)
                                    tgtgACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTAc                                                                                                                                 >  1:1067480/1‑125 (MQ=255)
                                    tgtgACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTAc                                                                                                                                 <  2:1067480/125‑1 (MQ=255)
                                    tgtgACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTAc                                                                                                                                 >  1:619514/1‑125 (MQ=255)
                                             aTAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAg                                                                                                                                            <  3:896558/105‑1 (MQ=255)
                                             aTAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAg                                                                                                                                            >  4:896558/1‑105 (MQ=255)
                                                        aaGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTAc                                                                                                                                 >  3:404113/1‑105 (MQ=255)
                                                        aaGAGGCACCCAGCTTTGGCGGGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTAc                                                                                                                                 <  4:404113/105‑1 (MQ=255)
                                                            ggCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGAGGTTAACAGTGCGAAATTACCgtg                                                                                                                             >  2:502455/1‑105 (MQ=255)
                                                            ggCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGAGGTTAACAGTGCGAAATTACCgtg                                                                                                                             <  1:502455/105‑1 (MQ=255)
                                                                   aGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAg                                                                                                         <  4:842130/118‑1 (MQ=255)
                                                                   aGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAg                                                                                                         >  3:842130/1‑118 (MQ=255)
                                                                        tGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCTGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTAttt                                                                      <  4:474357/148‑1 (MQ=255)
                                                                                agagGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAAt                                                                                                                    <  1:944632/94‑1 (MQ=255)
                                                                                agagGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAAt                                                                                                                    >  2:944632/1‑94 (MQ=255)
                                                                                  agGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGt                                                                                <  2:261163/128‑1 (MQ=255)
                                                                                  agGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGt                                                                                >  1:261163/1‑128 (MQ=255)
                                                                                     aTTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTaa                                                                    >  3:71045/1‑137 (MQ=255)
                                                                                     aTTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTaa                                                                    <  4:71045/137‑1 (MQ=255)
                                                                                     aTTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGAGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCAttt                                                         >  4:980429/1‑148 (MQ=255)
                                                                                          cgcgTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGaa                                                                                                                                      <  1:496039/66‑1 (MQ=255)
                                                                                          cgcgTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGaa                                                                                                                                      >  2:496039/1‑66 (MQ=255)
                                                                                          cgcgTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAAt                                                                                                                    <  4:1101396/84‑1 (MQ=255)
                                                                                          cgcgTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAAt                                                                                                                    >  3:1101396/1‑84 (MQ=255)
                                                                                            cgTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGaa                                                                                                                                      >  4:184648/1‑64 (MQ=255)
                                                                                            cgTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGaa                                                                                                                                      <  3:184648/64‑1 (MQ=255)
                                                                                              tCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGAGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCt                                                <  3:980429/148‑1 (MQ=255)
                                                                                                 caATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTa                                                                                                                                  >  3:70616/1‑63 (MQ=255)
                                                                                                 caATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGtt                                                                                                           <  2:627763/86‑1 (MQ=255)
                                                                                                 caATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGtt                                                                                                           >  1:627763/1‑86 (MQ=255)
                                                                                                 caATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGtt                                                                          >  1:1128895/1‑119 (MQ=255)
                                                                                                 caATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGtt                                                                          <  2:1128895/119‑1 (MQ=255)
                                                                                                 caATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATACGGGGTTAACAGTGCGAAATTa                                                                                                                                  <  4:70616/63‑1 (MQ=255)
                                                                                                  aaTTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGt                                                                                <  4:71905/112‑1 (MQ=255)
                                                                                                  aaTTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGt                                                                                >  3:71905/1‑112 (MQ=255)
                                                                                                   aTTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTAttt                                           >  1:632061/1‑148 (MQ=255)
                                                                                                           tgCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACaa                                                                                                                     >  4:780616/2‑66 (MQ=255)
                                                                                                           aGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACaa                                                                                                                     <  3:780616/66‑1 (MQ=255)
                                                                                                           aGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTaaa                                                                   <  1:519368/116‑1 (MQ=255)
                                                                                                           aGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTaaa                                                                   >  2:519368/1‑116 (MQ=255)
                                                                                                                     aataTTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCAAGTTGATa                         >  1:1019621/1‑148 (MQ=255)
                                                                                                                        atTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAg                                                                                                     >  4:580567/1‑69 (MQ=255)
                                                                                                                        atTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAg                                                                                                     <  3:580567/69‑1 (MQ=255)
                                                                                                                        atTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATg                                                      >  3:723550/1‑116 (MQ=255)
                                                                                                                        atTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATg                                                      <  4:723550/116‑1 (MQ=255)
                                                                                                                         ttttttACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCAAGTTGATAAAGc                     <  2:632061/148‑1 (MQ=255)
                                                                                                                              tACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCAAGTTGATAAAGCATGAg                <  2:1019621/148‑1 (MQ=255)
                                                                                                                                    aCCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCAAGTTGa                           <  4:822155/131‑1 (MQ=255)
                                                                                                                                    aCCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCAAGTTGa                           >  3:822155/1‑131 (MQ=255)
                                                                                                                                         ccGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCAAGTTGATAAAGCATGAGTGCCAACGttt     >  4:579918/1‑148 (MQ=255)
                                                                                                                                            gggTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCAAGTTGATAAAGCATGAGTGCCAACGTTTTGa  >  3:529353/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
ACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCAAGTTGATAAAGCATGAGTGCCAACGTTTTGA  >  NC_000913/543816‑544103

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