Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,502,968 Δ1 bp 13.5% coding (553/1380 nt) xanP → xanthine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,502,9680A.13.5% 92.7 / 12.0 37coding (553/1380 nt)xanPxanthine permease
Reads supporting (aligned to +/- strand):  ref base A (16/14);  new base . (3/2);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

GAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGATACCTG  >  NC_002947/4502827‑4503117
                                                                                                                                              |                                                                                                                                                     
gagaTGGGCGTGTCGCGCATCCTGTCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATTGGCCATGGCAcc                                                                                                                                               <  2:974910/150‑1 (MQ=255)
    tGGGCGTGTCGCGCATCCTGCCCTTCGTCACACGGATGTTGACCCCGCTTGTGTCTGGCATCGTGGTGTTGATGAGTGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTg                                                                                                                                           <  1:247841/150‑1 (MQ=255)
           gtCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATgg                                                                                                                                    <  1:924806/150‑1 (MQ=255)
             cgcgCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCg                                                                                                                                  <  2:139336/150‑1 (MQ=255)
                      tGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgc                                                                                                                         <  2:306036/150‑1 (MQ=255)
                      tGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTAACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGACAACctgc                                                                                                                         >  2:360685/1‑150 (MQ=255)
                               tCAAACGCCTGCTCACCCCGCTGTTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAACGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCgg                                                                                                                <  1:228303/150‑1 (MQ=255)
                                     gCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTAGGTGCCATTGGCCATGGCAccaaa                                                                                                                                            <  1:1071449/116‑4 (MQ=255)
                                     gCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTAGGTGCCATTGGCCATGGCAccaaa                                                                                                                                            >  2:1071449/1‑113 (MQ=255)
                                           tCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTgg                                                                                                    <  2:247578/150‑1 (MQ=255)
                                               cccGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATTGGCCATGGCAccaaa                                                                                                                                            >  1:448107/1‑103 (MQ=255)
                                               cccGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATTGGCCATGGCAccaaa                                                                                                                                            <  2:448107/106‑4 (MQ=255)
                                                                       gtTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgc                                                                                                                         >  1:102878/1‑101 (MQ=255)
                                                                       gtTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACctgc                                                                                                                         <  2:102878/101‑1 (MQ=255)
                                                                          gatgatCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggt                                                                     >  1:714684/1‑150 (MQ=255)
                                                                              atCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGgtggtg                                                                 >  2:691160/1‑150 (MQ=255)
                                                                                   ccTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGAt                                                            >  1:681068/1‑150 (MQ=255)
                                                                                         gCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGGGGGGTGGATGCGCAg                                                      >  1:1076129/1‑150 (MQ=255)
                                                                                             aTCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATAGCAGGTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTgc                                                  >  1:850595/1‑150 (MQ=255)
                                                                                              tCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGATGGTGTGGATGCGCAGCTgcg                                                 >  1:714565/1‑150 (MQ=255)
                                                                                                      ggCCTGATCAGCATGGGCGGTGGCTTGGGTGCCATGTCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTTTGCATGCGCAGCTGCGCCATCGTc                                         >  2:564684/1‑150 (MQ=255)
                                                                                                         cTGATCCGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAACCCTGCTGATGTCGGGCGTGGTACTGGCCATACCCCGAGCCCCCACCCGTATCCCGGTGGCGTGGCTCCCCAGCTGCGCCATCGTCATc                                      >  2:309215/1‑150 (MQ=255)
                                                                                                         cTGATAAGCATGGGCGGTGGCTTTTGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATc                                      <  1:279018/150‑1 (MQ=255)
                                                                                                                          tGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATc                                            >  1:573458/1‑127 (MQ=255)
                                                                                                                          tGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATc                                            <  2:573458/127‑1 (MQ=255)
                                                                                                                            gCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCg                                                                       <  1:760838/98‑1 (MQ=255)
                                                                                                                            gCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCg                                                                       >  2:760838/1‑98 (MQ=255)
                                                                                                                              tttGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCAt                                       <  1:809406/128‑1 (MQ=255)
                                                                                                                              tttGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCAt                                       >  2:809406/1‑128 (MQ=255)
                                                                                                                               ttGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATgcgc                <  1:763265/150‑1 (MQ=255)
                                                                                                                                 ggTGCCATTGGCCATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggtggt                                                                  <  1:390076/98‑1 (MQ=255)
                                                                                                                                 ggTGCCATTGGCCATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCggtggt                                                                  >  2:390076/1‑98 (MQ=255)
                                                                                                                                    gCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGcc           <  1:546536/150‑1 (MQ=255)
                                                                                                                                     ccATGGCCAATGGCACCTTTGCCCGTCCCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCg          >  1:705118/1‑150 (MQ=255)
                                                                                                                                     ccATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTa                                                                            <  2:922915/84‑1 (MQ=255)
                                                                                                                                     ccATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGAACTGGCGATCATCGTCGTTCTCAACCGTa                                                                            >  1:922915/1‑84 (MQ=255)
                                                                                                                                       aTGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGa        >  1:115576/1‑150 (MQ=255)
                                                                                                                                              aaTGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGATACCTg  >  1:628528/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
GAGATGGGCGTGTCGCGCATCCTGCCCTTCGTCAAACGCCTGATCACCCCGCTGGTGACTGGCATCGTGGTGTTGATGATCGGCCTGACGCTGATCAAGGTCGGCCTGATCAGCATGGGCGGTGGCTTTGGTGCCATGGCCAATGGCACCTTTGCCAATGGCGAAAACCTGCTGATGTCGGGCGTGGTACTGGCGATCATCGTCGTTCTCAACCGTATCCCGGTGGTGTGGATGCGCAGCTGCGCCATCGTCATCGCCCTGGCGGTCGGCTATGCGCTGGCCGGATACCTG  >  NC_002947/4502827‑4503117

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