Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 545,990 T→C Q125R (CAG→CGG)  allD ← ureidoglycolate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913545,9900TC97.7% 148.2 / NA 44Q125R (CAG→CGG) allDureidoglycolate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (22/21);  minor base A (1/0);  total (23/21)
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

GGATCTCGTCGCCTTCTCCCGGCGCGGCAAAGGCCAGGGGGTTAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCTGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCATGTTCCATACCCATTTTCGCCGCGACCTGTCCGGCGGCATTGTC  >  NC_000913/545853‑546120
                                                                                                                                         |                                                                                                                                  
ggATCTCGTCGCCTTCTCCCGGCGCGGCAAAGGCCAGGGGGTTAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCtgccgc                                                                                                                                  >  2:576320/1‑140 (MQ=255)
 gATCTCGTCGCCTTCTCCCGGCGCGGCAAAGGCCAGGGGGTTAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGca                                                                                                                                 <  1:332973/140‑1 (MQ=255)
   tctcGTCGCCTTCTCCCGGCGCGGCAAAGGCCAGGGGGTTAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGcaca                                                                                                                               <  2:83045/140‑1 (MQ=255)
           ccTTCTCCCGGCGCGGCAAAGGCCAGGGGGTTAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAaga                                                                                                                       >  2:605942/1‑140 (MQ=255)
           ccTTCTCCCGGCGCGGCAAAGGCCAGGGGGTTAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAaga                                                                                                                       >  1:23793/1‑140 (MQ=255)
            cTTCTCCCGGCGCGGCAAAGGCCAGGGGGTTAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATAGAAATGCCCATGAATCGGGCGCGGGCTGCCGGCCGCACAAAATAagag                                                                                                                      >  2:261492/1‑140 (MQ=255)
             ttCTCCCGGCGCGGCAAAGGCCAGGGGGTTAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGAGCAGCCTGCAGCACAAAATAAgaga                                                                                                                     >  1:260287/1‑140 (MQ=255)
                                gCCAGGGGGTTAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCAt                                                                                                  <  1:646516/140‑1 (MQ=255)
                                gCCAGGGGGTTAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCAt                                                                                                  <  1:480960/140‑1 (MQ=255)
                                     ggggTTAGTACCGTAGTGAATTGCCGCGCCGCCAAACGGGACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGc                                                                                             <  2:260287/140‑1 (MQ=255)
                                     ggggTTAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGc                                                                                             <  1:36793/140‑1 (MQ=255)
                                          tAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCGTGCCGCACGAAATAAGAGATTGCGCAGCTGGGACCCATCCGGCTGATa                                                                                        >  2:227975/1‑140 (MQ=255)
                                          tAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATa                                                                                        >  2:47621/1‑140 (MQ=255)
                                                 gtagtaAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCTacca                                                                                 >  1:609352/1‑140 (MQ=255)
                                                  tagtaAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGaccac                                                                                >  2:123815/1‑140 (MQ=255)
                                                  tagtaAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGaccac                                                                                <  1:605942/140‑1 (MQ=255)
                                                  tagtaAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGaccac                                                                                >  1:165760/1‑140 (MQ=255)
                                                     taAATTTCCGCGCCGCCAAACGCCATCACTATTGGATCGGGCTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACcgccgc                                                                             <  1:194671/140‑1 (MQ=255)
                                                          ttCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAGATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGc                                                                        >  2:52560/1‑140 (MQ=255)
                                                            ccGCGCCGCCAAACGCCACCACCATTGGACCGGACTGCCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCa                                                                      <  1:188362/140‑1 (MQ=255)
                                                             cgcgCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTCTGACCCATCCGACTGATACCGACCACCGCAACGCCAt                                                                     >  1:328206/1‑140 (MQ=255)
                                                               cgccgcCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATCGCGCCGCTCTGATCCACCCGGCTGAAACCGGCCACCGCAACGCCAttt                                                                   >  2:208971/1‑140 (MQ=255)
                                                               cgccgcCAAACGGCACCACCATTGGATCGGACGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCAtttt                                                                  >  2:83979/1‑140 (MQ=255)
                                                                      aaaCGGCACCACCATTGGATCGGACTGGCACTTCGAAATTCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTggg                                                            <  2:240608/140‑1 (MQ=255)
                                                                      aaaCGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCTACCGCCGCAACGCCATTTTGCTggg                                                            <  2:609352/140‑1 (MQ=255)
                                                                       aaCGGCACCACCATTCGATCGGTCTCGCACATCGAAACGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATCTTGCTGGGc                                                           <  1:200445/140‑1 (MQ=255)
                                                                             accaccATTGGATCGGACGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTg                                                    <  1:83979/140‑1 (MQ=255)
                                                                              ccaccaTTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTg                                                    <  2:547353/140‑1 (MQ=255)
                                                                                 ccaTTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACGGACCACCGCAATGCGATTTTGATGGGCGGTTTTGATg                                                 >  1:169682/1‑140 (MQ=255)
                                                                                   aTTGGATCGGACTGGCATATCGAAATGCCAATGAATCCGTCGCGGGCTACCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGc                                               <  2:285617/140‑1 (MQ=255)
                                                                                   aTTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCATAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGc                                               <  1:620243/140‑1 (MQ=255)
                                                                                     tGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCAt                                             <  1:47621/140‑1 (MQ=255)
                                                                                     tGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCAt                                             >  2:357751/1‑140 (MQ=255)
                                                                                     tGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCAACGCGACGCCATTTTGCTGGGCGGTTTTGATGGCAt                                             >  2:229914/1‑140 (MQ=255)
                                                                                            gACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCATGTTCCAt                                      >  1:118895/1‑140 (MQ=255)
                                                                                            gACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCATGTTCCAt                                      >  1:450722/1‑140 (MQ=255)
                                                                                            gACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCATGTTCCAt                                      >  2:593684/1‑140 (MQ=255)
                                                                                                 gCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCATGTTCCATACCCa                                 >  2:390754/1‑140 (MQ=255)
                                                                                                              ccAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATCAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCATGTTCCATACCCATTTTCGCCGCGAc                    <  1:261492/140‑1 (MQ=255)
                                                                                                              ccAATGAATCCGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCATGTTCCATACCCATTTTCGCCGCGAc                    <  1:97323/140‑1 (MQ=255)
                                                                                                                        cGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCTGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCATGTTCCATACCCATTTTCGCCGCGACCTGTCcggcg          <  2:118895/140‑1 (MQ=255)
                                                                                                                        cGGCGCGGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCATGTTCCATACCCATTTTCGCCGCGACCTGTCcggcg          >  1:59373/1‑140 (MQ=255)
                                                                                                                             cgGGCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTGGATGGCATGTTCGATACCCATTTTCCCCGCAACCTGTCCGGCGGCAtt     >  2:222106/1‑140 (MQ=255)
                                                                                                                               ggCTGCCTGCCGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCATGTTCCATACCCATTTTCGCCGCGACCTGTCCGGCGGCATTGt   >  2:537216/1‑140 (MQ=255)
                                                                                                                                gCTGCCTGCCGCACACAATAAGGGACTGCCCCGCTGTAACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCATGTTCCATACCCATTTTCGCCGCGACCTGTCCGGCGGCATTGTc  <  1:247431/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                  
GGATCTCGTCGCCTTCTCCCGGCGCGGCAAAGGCCAGGGGGTTAGTACCGTAGTAAATTTCCGCGCCGCCAAACGGCACCACCATTGGATCGGACTGGCACATCGAAATGCCAATGAATCCGGCGCGGGCTGCCTGCTGCACAAAATAAGAGATTGCGCCGCTGTGACCCATCCGGCTGATACCGACCACCGCAACGCCATTTTGCTGGGCGGTTTTGATGGCATGTTCCATACCCATTTTCGCCGCGACCTGTCCGGCGGCATTGTC  >  NC_000913/545853‑546120

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