Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,299,267:1 +A 11.1% intergenic (+22/‑1915) ychE → / → oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,299,2671.A11.1% 194.1 / 21.7 54intergenic (+22/‑1915)ychE/oppAUPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base . (18/30);  new base A (3/3);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTT  >  NC_000913/1299126‑1299413
                                                                                                                                              |                                                                                                                                                  
gTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATg                                                                                                                                             >  2:537890/1‑149 (MQ=255)
    tACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTcc                                                                                                                                         <  4:140264/149‑1 (MQ=255)
      cGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGcatggc                                                                                                                                                >  4:427402/1‑136 (MQ=255)
      cGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGcatggc                                                                                                                                                <  3:427402/141‑6 (MQ=255)
       gCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCATGGCAGCTCCGtt                                                                                                                                      <  4:202013/149‑1 (MQ=255)
       gCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGtt                                                                                                                                      <  2:500472/149‑1 (MQ=255)
          aGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTAc                                                                                                                                   <  4:261042/149‑1 (MQ=255)
          aGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTAc                                                                                                                                   <  4:225960/149‑1 (MQ=255)
               ggCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGt                                                                                                                              <  3:248397/149‑1 (MQ=255)
                     aaCGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCATGGCAGCTCCGtttcatttga                                                                                                                              <  4:408585/143‑8 (MQ=255)
                     aaCGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCATGGCAGCTCCGtttcatttga                                                                                                                              >  3:408585/1‑136 (MQ=255)
                     aaCGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCattat                                                                                                                        >  2:401304/1‑149 (MQ=255)
                     aaCGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCattat                                                                                                                        >  1:6123/1‑149 (MQ=255)
                      aCGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCattatt                                                                                                                       >  2:481421/1‑149 (MQ=255)
                      aCGTGATGACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCattatt                                                                                                                       >  2:399133/1‑149 (MQ=255)
                          gATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGACATGCTCCGTTTACTTCGTCATTATTTTTa                                                                                                                   <  1:440355/149‑1 (MQ=255)
                            ttACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACt                                                                                                                 <  1:481421/149‑1 (MQ=255)
                             tACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACtt                                                                                                                <  4:245490/149‑1 (MQ=255)
                                gcgTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGttt                                                                                                                                     <  1:164358/125‑1 (MQ=255)
                                gcgTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGttt                                                                                                                                     >  2:164358/1‑125 (MQ=255)
                                    attatGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTccc                                                                                                         <  3:276388/149‑1 (MQ=255)
                                                             tGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGaa                                                                                <  4:265834/149‑1 (MQ=255)
                                                             tGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGaa                                                                                <  3:415446/149‑1 (MQ=255)
                                                              ggggATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAAt                                                                               >  1:295293/1‑149 (MQ=255)
                                                                ggATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAg                                                                             <  3:9192/149‑1 (MQ=255)
                                                                      aaTTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTaaaa                                                                       <  1:537890/149‑1 (MQ=255)
                                                                         ttATTGTTACTGGTATTAGGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGATCCGTTTAAATCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAACGGAATAGGTAAAAata                                                                    <  3:410011/149‑1 (MQ=255)
                                                                          tATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAatat                                                                   <  4:23742/149‑1 (MQ=255)
                                                                               ttACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAttcttc                                                              <  2:305084/149‑1 (MQ=255)
                                                                                        ttAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAtt                                                                  <  2:200475/136‑1 (MQ=255)
                                                                                        ttAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAtt                                                                  >  1:200475/1‑136 (MQ=255)
                                                                                            ggggATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagtt                                                 >  2:578452/1‑149 (MQ=255)
                                                                                             gggATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagtta                                                >  3:442058/1‑149 (MQ=255)
                                                                                              ggATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagttag                                               >  2:586310/1‑149 (MQ=255)
                                                                                              ggATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagttag                                               <  4:3792/149‑1 (MQ=255)
                                                                                               gATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagttagt                                              >  1:163473/1‑149 (MQ=255)
                                                                                                 ttttCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACa                                                    >  3:130120/1‑141 (MQ=255)
                                                                                                 ttttCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATACTCTTCTCAAATTAcc                                                    <  4:130120/141‑2 (MQ=255)
                                                                                                    tCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATaa                                         <  3:311283/149‑1 (MQ=255)
                                                                                                            cTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTcc                                 <  3:407722/149‑1 (MQ=255)
                                                                                                                ttaattaaTTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTaa                             <  1:190085/149‑1 (MQ=255)
                                                                                                                 taattaatTCCTCCCAAATGATACGGATCTGCCACGCTCCGGTTACTTCATCATCAATATTACATTGTTCCCGCGTAGATAGCAAAAGCAAAAGGAATAGGTAAAAATATTGTTCCCAAATCACAGTTAGTTATAAGGATTTCCTCAAc                            <  4:331697/149‑1 (MQ=255)
                                                                                                                   attaattCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagttagtt                                             <  1:426897/130‑1 (MQ=255)
                                                                                                                   attaattCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagttagtt                                             >  2:426897/1‑130 (MQ=255)
                                                                                                                       attCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTtc                      <  3:10594/149‑1 (MQ=255)
                                                                                                                       attCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTtc                      <  1:399133/149‑1 (MQ=255)
                                                                                                                           cTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTctcct                  >  2:298952/1‑149 (MQ=255)
                                                                                                                            tttCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTctcctc                 >  1:40747/1‑149 (MQ=255)
                                                                                                                                  agtaaaCGGAGCATGGCAGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACcatc           >  3:212504/4‑149 (MQ=255)
                                                                                                                                   tGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACcatca          >  4:344056/1‑149 (MQ=255)
                                                                                                                                    gAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAATTTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACcatcat         <  3:564993/149‑1 (MQ=255)
                                                                                                                                     aaaCGGAGCATGGCAGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAttc                                                                 <  1:231321/92‑1 (MQ=38)
                                                                                                                                     aaaCGGAGCATGGCAGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAttc                                                                 >  2:231321/1‑92 (MQ=38)
                                                                                                                                         ggAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTAt    >  2:23491/1‑149 (MQ=255)
                                                                                                                                         ggAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTAt    <  3:446198/149‑1 (MQ=255)
                                                                                                                                           aGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTAttt  <  3:344056/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTT  >  NC_000913/1299126‑1299413

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