Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,010,122 G→A G32D (GGC→GAC)  hyuA → phenylhydantoinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,010,1220GA91.7% 64.5 / NA 24G32D (GGC→GAC) hyuAphenylhydantoinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/13);  minor base C (2/0);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAAAACCTGGAAATCTCACCTGGTGACCTGTGCCGCGATGTACGCTGCAATCCCGTTAAGCTGGCTGGCAACAGAGTAATTACTCTTCACTTATTCCCCTCCGGTTCGCCGGAGGGTTTTTGGAGTTTGCTATGCGCGTATTGATCAAAAACGGCACTGTCGTTAACGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGGCAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTCCCTGGCGGCGTGGATGTCCATACGCATTTCAATATTGATGTCGGCATCGCGCGCAGTTGTGATGATTTTTTTACCGGTACCCGCGCAGCTGCGTGTGGCGGTACAACAACCATTATTGACCATATGGGATTTGGCCCAAACGG  >  NC_000913/3009895‑3010326
                                                                                                                                                                                                                                   |                                                                                                                                                                                                            
gAAAAAACCTGGAAATCTCACCTGGTGACCTGTGCCGCGATGTACGCTGCAATCCCGTTAAGCTGGCTGGCAACAGAGTAATTACTCTTCACTTATTCCCCTCCGGTTCGCCGGAGGGTTTTTGGAGTTTGCTATGCGCGTATTGATCAAAAACGGCACTGTCGTTAACGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTnc                                                                                                                                                                                       <  1:347157/251‑1 (MQ=255)
          ggAAATCTCACCTGGTGACCTGTACCGCGATGTACGCTGAACTCCCGTTAAGCTGCCTGACAAAAGACAAATTACTCTTCACTTCTTCCCCTGCGATTCACCGGATGGCTTTTGGAGCTTGCTATCCGCGTATTGATCAAAAAAGACATTGTCGTTAACGCAGATGGACAAGCAAAACAGGAGTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGCCAACAATATTTCTCCGCAGCTCCCGTATga                                                                                                                                                                                >  1:434336/1‑248 (MQ=255)
                    ccTGGTGACCTGTGCCGCGATGTACGCTGCAATCCCGTTAAGCTGGCTGGCAACCGAGTCATTACTCTTCAATTATTCCCCTCCGGTTCGCCGGAGGGTTTTTGGAGTTTGATATGCGCGTATTGATCAAAAACGGCACTGTCGTTAACGCAGATGGACACGCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGcc                                                                                                                                                                   >  1:314413/1‑251 (MQ=255)
                            ccTGTGCCGCGATGTACGCTGCAATCCCGTTAAGCTGGCTGGCAACAGAGTAATTACTCTTCACGTATTCCCCTCCGGTTCGCCGGAGGGTTTTTGGAGTTTGCTATGCGCGTATTGATCAAAAACGGCACTGTCGTTAACGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGAGTGTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCtn                                                                                                                                                           <  1:129440/251‑2 (MQ=255)
                                           aCGCTGCAATCCCGTTAAGCTGGCTGGCAACAGAGTAATTACTCTTCACTTATTCCCCTCCGGTTCGCCGGAGGGTTTTTGGAGTTTGCTATGCGCGTATTGATCAAAAACGGCACTGTCGTTAACGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGTACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTCCCTgg                                                                                                                                            <  1:335478/251‑1 (MQ=255)
                                                        gTTAAGCTGGCTGGCAACAGAGTAATTACTCTTCACTTATTCCCCTCCGGTTCGCCGGAGGGTTTTTGGAGTTTGCTATGCGCGTATTGATCAAAAACGGCACTGTCGTTAACGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTCCCTGGCGGCGTGGATGTc                                                                                                                               <  1:238507/251‑1 (MQ=255)
                                                                ggctggcAACAGAGGAATTACTCTTCACTTATTCCCCCCCGGTTCGCCGGAGGGTTTTTGGAGTTTGCTATGCGCGTATTGATCAAAAACGGCACTGTCGGTAACGCAGATGGACAAGCCAACCAGCATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGACAACAATATTGCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGATTTCCCTGGCGGCGTGGATGTCCATACGc                                                                                                                        <  1:262706/250‑1 (MQ=255)
                                                                                    ctctTCACTTATTCCCCTCCGGTTCGCCGGAGGGTTTTTGGAGTTTGCTATGCGCGTATTGATCAAAAACGGCACTGTCGTTAACGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGAGTGTTCGCCAGTTGGATAACAATATTTcgcc                                                                                                                                                                                               <  1:124954/159‑1 (MQ=255)
                                                                                                                           ggAGTTTGCTATGCGCGTATTGATCAAAAACGGCACTGTCGTTAACGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTCCCTGGCGGCGTGGATGTCCATACGCATTTCAATATTGATGTCGGCATCGCGCGCAGTTGTGATGATTTTTTTACCGGTACCcgcg                                                            >  1:177851/1‑251 (MQ=255)
                                                                                                                                       gcgcGTATTGATCAAAAACGGCACTGTCGTTAACGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTCCCTGGCGGCGTGGATGTCCATACGCATTTCAATATTGATGTCGGCATCGCGCGCAGTTGTGATGATTTTTTTACCGGTACCCGCGCAGCTGCGtgtg                                                <  1:2814/251‑1 (MQ=255)
                                                                                                                                                    aaaaaCGGCACTGTCGTAAACGCAGATGGGCAAGCCAAACAGGATTTGCCGATTGAAAGCGGGATTGTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGGTTTACCTTGCGGCGTGGATGTCCATACGCCATTCAATATTGATGTCGGCATCGCGCGGAGTTGTGATGATTTTTTTACCGGTAACCGCGCAGCTGCGTGGGGCGGTAcaacaa                                    >  1:227943/1‑250 (MQ=255)
                                                                                                                                                            cACTGTCGTTAACGCAGATGGACAAGCCCAACAGGATTTGCTGATTGAAAGCGGCATTGACCGCCAGTTGGACAACAATATTTCGCCGCACCTCCCCTATGAAGAAATTGATGCAACTGGCTATTAAGTTTTTCCTGGCAGCGTGGATGTCCATACGCATTTCAATATTTATGTCGGCATCGCGCGCAGTTGTGATGAtttt                                                                            >  1:68461/1‑202 (MQ=255)
                                                                                                                                                             aCTGTCGTTAACGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGACAACACTATTTcgccgc                                                                                                                                                                                             >  2:31018/1‑88 (MQ=255)
                                                                                                                                                             aCTGTCGTTAACGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGACAACAATATTTcgccgc                                                                                                                                                                                             <  1:31018/88‑1 (MQ=255)
                                                                                                                                                                   gTTAACGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGATGGTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCtn                                                                                                                                                           <  1:78738/116‑2 (MQ=255)
                                                                                                                                                                        cGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTCCCTGGCGGCGTGGATGTCCATACGCATTTCAATATTGATGTCGGCAGCGCGCGCAGTtgtg                                                                                    <  1:124368/182‑1 (MQ=255)
                                                                                                                                                                          cAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCACCAGTTAGCCAACAATATTTCGCAGCAGCTCCAGCATGAAGAAATTGCTGACACTGGCTGTTACGTTTTCCCTGGCGGCGTGGATGTCCATACGCATTTCAATATTGATGTCGGCATCGCACGCAGTTGTGATGATTTTTTTACAGGTACCCGCGCAGCTGCGTGTGGCGATACAACAACCATTATTGACCATATGGGAttt             >  1:300789/1‑251 (MQ=255)
                                                                                                                                                                                     gCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTCCCTGGCGGCGTGGATGTCCATACGCATTTCAATATTGATGTCGGCATCGCGCGCAGTTGTGATGATTTTTTTACCGGTACCCGCGCAGCTGCGTGTGGCGGTACAACAACCATTATTGACCATATGGGATTTGGCCCAAACgg  >  1:354293/1‑251 (MQ=255)
                                                                                                                                                                                           cAGGATTTGCTGATTGCAAGCGGGATTGTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTCCCTGGCGGCGTGGATGTCCATAc                                                                                                                          >  1:155859/1‑125 (MQ=255)
                                                                                                                                                                                               aTTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTCCCTGGCGGCGTGGAGGTCCATACGCATTTCAATATTGATGTCGGCATCGCGCGCAGTTGTGATGATTTTTTTACCGGTAccc                                                               <  1:243869/180‑1 (MQ=255)
                                                                                                                                                                                                                      nttCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTccc                                                                                                                                               >  1:267479/2‑77 (MQ=255)
                                                                                                                                                                                                                      gTTCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTccc                                                                                                                                               <  2:267479/77‑1 (MQ=255)
                                                                                                                                                                                                                       ttCGCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTCCCTGGCGGCGTGGATGTCCATACGCATTTCAATATTGATGTCGGCATCGCGCGCAGTTGTGATGATTTTTTTACCGGTACCCGCGCAGCTGCGTGTGGCGGTACAACAACCATTATTGACCATATGGGATTTGGc          <  1:14551/209‑1 (MQ=255)
                                                                                                                                                                                                                          gCCAGTTGGACAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTCCCTGGCGGCGTGGATGTCCAtcc                                                                                                                          >  1:372995/1‑92 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                            
GAAAAAACCTGGAAATCTCACCTGGTGACCTGTGCCGCGATGTACGCTGCAATCCCGTTAAGCTGGCTGGCAACAGAGTAATTACTCTTCACTTATTCCCCTCCGGTTCGCCGGAGGGTTTTTGGAGTTTGCTATGCGCGTATTGATCAAAAACGGCACTGTCGTTAACGCAGATGGACAAGCCAAACAGGATTTGCTGATTGAAAGCGGGATTGTTCGCCAGTTGGGCAACAATATTTCGCCGCAGCTCCCGTATGAAGAAATTGATGCCACTGGCTGTTACGTTTTCCCTGGCGGCGTGGATGTCCATACGCATTTCAATATTGATGTCGGCATCGCGCGCAGTTGTGATGATTTTTTTACCGGTACCCGCGCAGCTGCGTGTGGCGGTACAACAACCATTATTGACCATATGGGATTTGGCCCAAACGG  >  NC_000913/3009895‑3010326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: