Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,881,464 G→A 8.9% V183I (GTA→ATA)  uacT → uric acid permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,881,4640GA 8.9% 126.6 / 2.2 45V183I (GTA→ATA) uacTuric acid permease
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base A (2/2);  total (22/23)
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.62e-01

CCCACCGCTGGTCACCGGTACGGTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCTACATCCTGGCCGGTTCCATCGGTATGGTCGACCTGTCGGGCCTGAATGAAGCACCGTGGTTGC  >  NC_002947/4881331‑4881611
                                                                                                                                     |                                                                                                                                                   
cccACCGCTGGTCACCGGTACGGTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATcc                                                                                                                                      >  6:180630/1‑149 (MQ=255)
cccACCGCTGGTCACCGGTACGGTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATc                                                                                                                                       <  5:140355/148‑1 (MQ=255)
   aCCGCTGGTCACCGGTACGGTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCctgc                                                                                                                                   <  6:102451/149‑1 (MQ=255)
     cGCTGGTCACCGGTACGGTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGa                                                                                                                                         <  1:139136/141‑1 (MQ=255)
     cGCTGGTCACCGGTACGGTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGa                                                                                                                                         >  2:139136/1‑141 (MQ=255)
          gTCACCGGTACGGTAATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCaa                                                                                                                            >  1:49203/1‑149 (MQ=255)
             aCCGGTACGGTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCg                                                                                                                         >  8:207951/1‑149 (MQ=255)
              ccGGTACGGTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGATACTGGCCGTGATCCTGCTGAt                                                                                                                               >  1:256940/1‑142 (MQ=255)
              ccGGTACGGTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGATACTGGCCGTGATCCTGCTGAt                                                                                                                               <  2:252950/142‑1 (MQ=255)
              ccGGTACGGTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGATACTGGCCGTGATCCTGCTGAt                                                                                                                               <  2:256940/142‑1 (MQ=255)
              ccGGTACGGTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGATACTGGCCGTGATCCTGCTGAt                                                                                                                               >  1:252950/1‑142 (MQ=255)
                      gTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTgg                                                                                                                                                     <  8:49661/112‑1 (MQ=255)
                      gTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTgg                                                                                                                                                     >  7:49661/1‑112 (MQ=255)
                         aTCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCgg                                                                                                             <  7:33874/149‑1 (MQ=255)
                         aTCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCgg                                                                                                             <  7:207951/149‑1 (MQ=255)
                            aCCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGttt                                                                                                          <  8:60565/149‑1 (MQ=255)
                             ccTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTc                                                                                                         >  8:146144/1‑149 (MQ=255)
                               tCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTg                                                                                                       >  6:262210/1‑149 (MQ=255)
                                 gATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTg                                                                                                       <  5:262210/147‑1 (MQ=255)
                                         tGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACgtg                                                                                             >  1:257533/1‑149 (MQ=255)
                                           tCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGa                                                                                                                                         >  7:15995/1‑103 (MQ=255)
                                           tCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGa                                                                                                                                         <  8:15995/103‑1 (MQ=255)
                                           tCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAAc                                                                                                >  7:82030/1‑144 (MQ=255)
                                           tCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAAc                                                                                                <  8:82030/144‑1 (MQ=255)
                                            cGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCg                                                                                          >  1:268356/1‑149 (MQ=255)
                                                    ccGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACtggt                                                                                  >  6:84418/1‑149 (MQ=255)
                                                     cGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACtggtg                                                                                 >  4:61526/1‑149 (MQ=255)
                                                      ggTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACtggtgg                                                                                <  2:268356/149‑1 (MQ=255)
                                                             gTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACgt                                                                                              <  3:2474/128‑1 (MQ=255)
                                                             gTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACgt                                                                                              >  4:2474/1‑128 (MQ=255)
                                                              tCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTgg                                                                                                      <  4:55910/119‑1 (MQ=255)
                                                              tCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATgg                                                                          <  5:84418/147‑1 (MQ=255)
                                                                     ggCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGc                                                                  >  2:136165/1‑148 (MQ=255)
                                                                     ggCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCt                                                                 <  7:146144/149‑1 (MQ=255)
                                                                          gTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCTACATc                                                            <  1:136165/149‑1 (MQ=255)
                                                                             gcggcCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCTACATCCt                                                          <  2:49203/148‑1 (MQ=255)
                                                                               ggcCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTg                                                                    <  6:58463/136‑1 (MQ=255)
                                                                               ggcCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAAAGTGTCGGTACTGGTGGGCATGGGCCTTg                                                                    >  5:58463/1‑136 (MQ=255)
                                                                                  cATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCTACATCCTGGCCgg                                                    >  8:17121/1‑149 (MQ=255)
                                                                                    tGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTCATCCTGCTGATCAACCGCTTCATGCGAGGTGTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCTACATCCTGGCCGGtt                                                  >  5:180209/1‑149 (MQ=255)
                                                                                             aaCGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCTACATCCTGGCCGGTTCCATCGGTa                                         <  1:3830/149‑1 (MQ=255)
                                                                                                    ggCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCTACATCCTGGCCGGTTCCATCGGTATGGTCGa                                  >  3:218190/1‑149 (MQ=255)
                                                                                                        cGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCTACATCCTGGCCGGTTCCATCGGTATGGTCGACCTg                              >  5:272920/1‑149 (MQ=255)
                                                                                                         gCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCTACATCCTGGCCGGTTCCATCGGTATGGTCGACCt                               <  3:61526/147‑1 (MQ=255)
                                                                                                         gCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCTACATCCTGGCCGGTTCCATCGGTATGGTCGACCTGt                             <  6:272920/149‑1 (MQ=255)
                                                                                                                              cggccTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCg                                                                                          >  7:204712/1‑67 (MQ=255)
                                                                                                                              cggccTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCg                                                                                          <  8:204712/67‑1 (MQ=255)
                                                                                                                                     gTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCTACATCCTGGCCGGTTCCATCGGTATGGTCGACCTGTCGGGCCTGAATGAAGCACCGTGGTTGc  <  5:180630/148‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                   
CCCACCGCTGGTCACCGGTACGGTGATCACCTCGATTGGCCTGTCGCTGTTCCCGGTGGCTGTCAACTGGGCCGGTGGCGGCCATGAAGCAGAAACGTTCGGCTCGCCGATCTACCTGCTGGTGGCCGGCCTGGTACTGGCCGTGATCCTGCTGATCAACCGCTTCATGCGCGGTTTCTGGGTAAACGTGTCGGTACTGGTGGGCATGGGCCTTGGCTACATCCTGGCCGGTTCCATCGGTATGGTCGACCTGTCGGGCCTGAATGAAGCACCGTGGTTGC  >  NC_002947/4881331‑4881611

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