Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 789,991 C→T 7.4% G347R (GGA→AGA)  galT ← galactose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913789,9910CT 7.4% 139.3 / 2.7 54G347R (GGA→AGA) galTgalactose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (2/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

AGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGATAAAAGTG  >  NC_000913/789849‑790138
                                                                                                                                              |                                                                                                                                                   
aGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTc                                                                                                                                               >  1:844463/1‑149 (MQ=255)
    cGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGa                                                                                                                                           >  1:189577/1‑149 (MQ=255)
      aaaCCGGCGTTGTAGGCGGTGTGCTCACCACTCAAATTCACGCGGCGCGGCGCCTGAATGGTGGGAGTGGCAGGGCAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGaaa                                                                                                                                         <  1:651380/149‑1 (MQ=255)
            tcgtTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGat                                                                                                                                   <  1:68603/149‑1 (MQ=255)
                    gTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGa                                                                                                                           >  1:820098/1‑149 (MQ=255)
                        gtgtgtTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTgc                                                                                                                       >  1:951610/1‑149 (MQ=255)
                             ttCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGCCTGCGCGCa                                                                                                                  >  2:60990/1‑149 (MQ=255)
                                aCCAATCAAATTCCCGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTCGCCAAATGCGTTGGCAAACAGAGAGTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGTAAATGGATATCGCTGACTGCGCGCAAAc                                                                                                               >  2:334825/1‑149 (MQ=255)
                                aCCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTATTATTTCAGACTCATTTCTTACACTCCGGAGTCGCGAAAATGGATAGCGCTGACTGCGCGCAAAc                                                                                                               >  2:678462/1‑149 (MQ=255)
                                aCCAATCAAATTCACGCGGCCAGGCGCCGGAAGGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAAc                                                                                                               <  1:1022913/149‑1 (MQ=255)
                                   aaTCAAATTCACGCGGCCCGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGct                                                                                                            >  1:192847/1‑149 (MQ=255)
                                     tCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCTGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGctct                                                                                                          >  2:467421/1‑149 (MQ=255)
                                     tCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCTGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGctct                                                                                                          >  2:118996/1‑149 (MQ=255)
                                      cAAATTCACGCGGACAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTctg                                                                                                         <  2:253203/149‑1 (MQ=255)
                                        aatgctcgcgGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGAAATCGCTGACTGCGCGCATACGCTctgct                                                                                                       <  2:58541/143‑1 (MQ=255)
                                             aCGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTg                                                                                                  >  2:631215/1‑149 (MQ=255)
                                               gcgGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGAAGCCAAATGCGTAGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTCCGCGCAAACGCTCTGCTGCCTGtt                                                                                                >  1:702973/1‑149 (MQ=255)
                                               gcgGCCAGGCGCCGGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGTTTGTGGTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGtt                                                                                                <  1:586640/149‑1 (MQ=255)
                                                 ggCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCG‑AAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTg                                                                                             <  2:890400/149‑1 (MQ=255)
                                                      ggagcCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCTGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGt                                                                                         <  1:467421/146‑1 (MQ=255)
                                                      ggCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCTGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGt                                                                                         <  1:118996/149‑1 (MQ=255)
                                                                  gtgtGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTg                                                                             <  1:147659/149‑1 (MQ=255)
                                                                   tgtgAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTgg                                                                            <  1:631215/149‑1 (MQ=255)
                                                                             aGGGTAGCCAAATGGGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCa                                                                  <  2:964802/149‑1 (MQ=255)
                                                                               ggTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGc                                                                >  2:764602/1‑149 (MQ=255)
                                                                                  aGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCAtt                                                             <  2:820098/149‑1 (MQ=255)
                                                                                  aGCCAAATGCGTTGGCAAACAGAGATTGGGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCAtt                                                             >  1:289494/1‑149 (MQ=255)
                                                                                     cAAATGAGATGGCAAACAGCGAATGTGTTTAATCTTTCAGACTCAATTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGCCAGGTCTCGCTGGGGCTCTGCCAGCATTTCa                                                          <  1:445515/149‑1 (MQ=255)
                                                                                      aaaTGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGAGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCAt                                                         >  2:926417/1‑149 (MQ=255)
                                                                                       aaTGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGCCTGCGCGCAAACGCTCTGCTGCCTGTGCGGCGGTCAGGTCTCGGTGGGTCGCGGCGAGCATTTCAtt                                                        >  2:611138/1‑148 (MQ=255)
                                                                                        aTGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGGCTCTGCCAGCATTTCATaa                                                       <  2:936524/149‑1 (MQ=255)
                                                                                                 aaaaCTGAGATTGTGTTTTTTCTTTAAGAATCATTTCTTACCCTCCGGATTAGCGAAAATGGATTTCGCTGACTGCGCGCAAACGCTCTTCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATa                                              <  2:734547/148‑1 (MQ=255)
                                                                                                  aaaCAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATaa                                             >  2:519460/1‑149 (MQ=255)
                                                                                                  aaaCAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATAGATGACTGCGCGCAAACGCTCTGCTGCCTGGTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCAGTTCATAACCAACCATaa                                             >  2:589588/1‑149 (MQ=255)
                                                                                                   aaCAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATaaa                                            <  2:624622/149‑1 (MQ=255)
                                                                                                       gagaTTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTta                                        >  2:384207/1‑149 (MQ=255)
                                                                                                        agaTTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTtac                                       >  1:839025/1‑149 (MQ=255)
                                                                                                          aTTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTtacgt                                     >  2:1005075/1‑149 (MQ=255)
                                                                                                            tgtgTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTtacgtac                                   >  1:988978/1‑149 (MQ=255)
                                                                                                             gtgtTTTTTCTTTCAGACTCATTTCTTTCACTCCGGATTCGCGAAAAAGGATATCGCTGACTGCGCGCAAACTCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGTGTCTCTGCCAGCATTTCATAAGCAACGATAAATTtacgtacg                                  <  2:127721/149‑1 (MQ=255)
                                                                                                             gtgtTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCATAAGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTtacgtacg                                  <  2:538775/149‑1 (MQ=255)
                                                                                                              tgtTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATTTCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACgg                                 <  2:143202/149‑1 (MQ=255)
                                                                                                              tgtTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGGCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACgg                                 <  2:844463/149‑1 (MQ=255)
                                                                                                                      atttcAGACTCATTTCTTACCCTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCTCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAAATTACGTACGGTGGCGGAg                         <  2:1029229/148‑1 (MQ=255)
                                                                                                                       tttcAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAgc                        >  2:72772/1‑149 (MQ=255)
                                                                                                                            gACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAGCgcagc                   >  1:932122/1‑149 (MQ=255)
                                                                                                                                     cTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCggt          >  2:189627/1‑148 (MQ=255)
                                                                                                                                     cTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGa          >  1:245246/1‑149 (MQ=255)
                                                                                                                                     cTTACACTCCGGATCCGCGAAAATGGATATCGCTGCCTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGa          <  1:851323/149‑1 (MQ=255)
                                                                                                                                      ttACACTCCGGAGTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGAt         <  1:524158/149‑1 (MQ=255)
                                                                                                                                        acacTCCGGATTAGCGAAAATGGATATCGCTTACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCGAGCATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGATaa       >  2:493229/1‑149 (MQ=255)
                                                                                                                                          acTCCGGATTCGCGAAAATGGCTATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGATaaaa     <  1:239920/149‑1 (MQ=255)
                                                                                                                                          acTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGATaaaa     <  1:764602/149‑1 (MQ=255)
                                                                                                                                            tCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGATAAAAGt   >  1:303598/1‑149 (MQ=255)
                                                                                                                                             ccGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGAATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGATAAAAGTg  >  1:647233/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
AGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGATAAAAGTG  >  NC_000913/789849‑790138

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