Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,508,231 G→A 11.5% intergenic (+85/‑165) cycA → / → urtA D‑alanine/beta‑alanine/D‑serine/glycine permease/urea ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,508,2310GA11.5% 136.0 / 7.9 52intergenic (+85/‑165)cycA/urtAD‑alanine/beta‑alanine/D‑serine/glycine permease/urea ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base A (3/3);  total (27/25)
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.25e-01

CCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGCCCAACACCATAA  >  NC_002947/5508095‑5508374
                                                                                                                                        |                                                                                                                                               
ccAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCtgt                                                                                                                                     >  3:286692/1‑149 (MQ=255)
    cTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGca                                                                                                                                  >  5:269593/1‑148 (MQ=255)
     tCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGcaca                                                                                                                                >  7:47807/1‑149 (MQ=255)
       gTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGcaca                                                                                                                                <  3:229594/147‑1 (MQ=255)
        tACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACaa                                                                                                                             <  6:269593/149‑1 (MQ=255)
            cagcagGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATcc                                                                                                                         <  1:178429/149‑1 (MQ=255)
             agcagGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATcc                                                                                                                         >  8:266884/1‑148 (MQ=255)
                agGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATc                                                                                                                                         >  2:181748/1‑129 (MQ=255)
                agGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATc                                                                                                                                         <  1:181748/129‑1 (MQ=255)
                 ggCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAg                                                                                                                    >  2:208894/1‑149 (MQ=255)
                             gCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTgg                                                                                                        <  7:266884/149‑1 (MQ=255)
                              cctatctgTGCGTAACGGCTGATGGCCTCGGGGGCTGTTGCCACGAAGGCACTGTGTACGCCGTCCTATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTgg                                                                                                        <  3:260838/141‑1 (MQ=39)
                               aCGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGAt                                                                                                      <  4:15038/149‑1 (MQ=255)
                               aCGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGAt                                                                                                      <  4:286692/149‑1 (MQ=255)
                                    cGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATccaac                                                                                                 <  5:269119/149‑1 (MQ=255)
                                    cGTGCGTAACGGCTGATGCCCTCGGGGGCGGTTCCCACGAACGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATTGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATccaac                                                                                                 <  6:43095/149‑1 (MQ=255)
                                      tGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCATGTGCACACAAATCAATCAGGCACCAACCTAGaagatctg                                                                                                >  4:260838/1‑141 (MQ=39)
                                        cGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCcaaccaac                                                                                             <  2:295675/149‑1 (MQ=255)
                                               gCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCaa                                                                                      >  8:165560/1‑149 (MQ=255)
                                                 tGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCa                                                                                                                     >  3:239821/1‑116 (MQ=255)
                                                 tGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCa                                                                                                                     <  4:239821/116‑1 (MQ=255)
                                                 tGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACAGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGa                                                                                    >  7:223120/1‑149 (MQ=255)
                                                      gCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACTAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCg                                                                               <  8:47807/149‑1 (MQ=255)
                                                        cTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGCTCCAACCCACGCCCCAAGACCTCGcc                                                                             >  8:114919/1‑148 (MQ=255)
                                                        cTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGca                                                                             <  1:322919/149‑1 (MQ=255)
                                                              ggCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAgg                                                                                                                   <  3:161782/105‑1 (MQ=255)
                                                              ggCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAgg                                                                                                                   >  4:161782/1‑105 (MQ=255)
                                                              ggCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGACCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGccc                                                                       >  6:73670/1‑149 (MQ=255)
                                                              ggCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGACCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGccc                                                                       >  6:69042/1‑149 (MQ=255)
                                                                      ccACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTaaa                                                                <  7:165560/148‑1 (MQ=255)
                                                                      ccACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTaaa                                                                <  1:208894/148‑1 (MQ=255)
                                                                       cACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCACCCAACGCCCCATGACCTCGCCCAGCCCCTCGAAAcc                                                              >  4:86887/1‑149 (MQ=255)
                                                                           aaGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTAc                                                          >  7:141659/1‑149 (MQ=255)
                                                                                aCGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCCACCTGGGTCCCACCCACGCCCCAATCCCTCGCACCGCCCCTCTATACCGTACCTATc                                                     >  8:117993/1‑149 (MQ=38)
                                                                                         gCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCACCCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAgg                                            >  4:176817/1‑149 (MQ=255)
                                                                                                gATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTATCGCCACAAGGCACAGcc                                     >  1:294881/1‑149 (MQ=255)
                                                                                                 aTGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACAGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGccc                                    >  1:215354/1‑149 (MQ=255)
                                                                                                   gtgGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCtt                                  >  4:93621/1‑149 (MQ=255)
                                                                                                        aGCGCCCTTGTGTCGCGATTGGCCTCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAAc                             <  8:223120/149‑1 (MQ=255)
                                                                                                          cgcCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGGACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCtt                                  <  7:265595/142‑1 (MQ=255)
                                                                                                          cgcCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCAGGATCCAACCAACGCCCCAAGACCTCGCACACCCCCTATAAACCGTACCTTACGCCACTACCAACAGACCtt                                  >  8:265595/1‑142 (MQ=255)
                                                                                                           gcCCTTGTGTCGCGATGGGCCGCCCCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGcc                          >  6:248737/1‑149 (MQ=255)
                                                                                                           gcCCTTGTGTCGCGATGGGCCGCACCGCGACCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGcc                          <  5:69042/149‑1 (MQ=255)
                                                                                                           gcCCTTGTGTCGCGATGGGCCGCACCGCGACCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGcc                          <  5:73670/149‑1 (MQ=255)
                                                                                                                 gtgtCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTccc                     <  7:114919/148‑1 (MQ=255)
                                                                                                                   gtCGCGATGGGCCGCACCGCGACCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAAc                                                               >  7:227934/1‑104 (MQ=255)
                                                                                                                   gtCGCGATGGGCCGCACCGCGACCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAAc                                                               <  8:227934/104‑1 (MQ=255)
                                                                                                                      gcgATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCtt                <  2:215354/148‑1 (MQ=255)
                                                                                                                           gggCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGCCCaa          >  1:166695/1‑149 (MQ=255)
                                                                                                                            ggCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGc              >  7:147671/1‑144 (MQ=255)
                                                                                                                            ggCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGc              <  8:147671/144‑1 (MQ=255)
                                                                                                                                    cgcgGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCCAGACCTCGCACAGCCCCTCTAAACCGTACATTTCGCCACTAGGCACAGCCCTTGCAACCCCCTCCACGCTTGCCCCACACCATaa  >  4:197539/1‑148 (MQ=255)
                                                                                                                                        |                                                                                                                                               
CCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGCCCAACACCATAA  >  NC_002947/5508095‑5508374

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