Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,571,404 A→G 19.1% intergenic (+172/‑1607) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,571,4040AG19.1% 277.7 / 58.1 137intergenic (+172/‑1607)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (45/65);  new base G (15/11);  total (60/77)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 148 total aligned reads displayed.

TTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAG  >  Exported/2571334‑2571473
                                                                      |                                                                     
ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                                                                       <  1:1616543/71‑1 (MQ=255)
   ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                                     >  1:2060905/1‑70 (MQ=255)
          cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                              <  1:5004618/70‑1 (MQ=255)
          cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                              <  1:3511341/70‑1 (MQ=255)
          cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                              <  1:5097885/70‑1 (MQ=255)
           tGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                              <  1:774512/69‑1 (MQ=255)
             gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                              <  1:5436525/67‑1 (MQ=255)
             gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                            <  1:4143204/69‑1 (MQ=255)
             gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                            <  1:4358569/69‑1 (MQ=255)
             gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                            <  1:666272/69‑1 (MQ=255)
               tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa                                                                   <  1:940354/60‑1 (MQ=255)
               tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa                                                                   <  1:1392556/60‑1 (MQ=255)
               tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa                                                                   <  1:3179875/60‑1 (MQ=255)
                gAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                                       >  1:1537116/1‑71 (MQ=255)
                gAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                                       >  1:463623/1‑71 (MQ=255)
                 gctgtgATCTCATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTgagaa                                                         >  1:1551697/3‑67 (MQ=255)
                 aGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                                       <  1:1820505/70‑1 (MQ=255)
                  gtgtGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCa                                                     <  1:227182/71‑1 (MQ=255)
                     tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                                  <  1:3810209/71‑1 (MQ=255)
                     tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                                  <  1:4960781/71‑1 (MQ=255)
                     tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                                  <  1:1268388/71‑1 (MQ=255)
                     tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                                   <  1:4984534/70‑1 (MQ=255)
                      gCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag                                                           <  1:1624458/61‑1 (MQ=255)
                      gCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag                                                           <  1:417316/61‑1 (MQ=255)
                      gCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag                                                           <  1:4165716/61‑1 (MQ=255)
                       ctctCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                                                <  1:151787/71‑1 (MQ=255)
                       atctcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACg                                                       <  1:1855778/63‑1 (MQ=255)
                       atctcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACg                                                       <  1:4993226/63‑1 (MQ=255)
                        tctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                                     >  1:2847201/1‑49 (MQ=255)
                         ctcATTAGATTTCTCATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatga                                              >  1:1943869/1‑71 (MQ=255)
                         ctcATTAGATTTCTCATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatga                                              >  1:2317238/1‑71 (MQ=255)
                          tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                              >  1:4314960/1‑70 (MQ=255)
                          tcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCacta                                                  <  1:4315466/66‑4 (MQ=255)
                          tcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCacta                                                  <  1:1135986/66‑4 (MQ=255)
                           cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                                     >  1:1706652/1‑46 (MQ=255)
                           cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                               <  1:3904393/68‑1 (MQ=255)
                           cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                               <  1:1937277/68‑1 (MQ=255)
                           cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGata                                            >  1:2811266/1‑71 (MQ=255)
                            aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                              <  1:552476/68‑1 (MQ=255)
                            aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                              <  1:3581123/68‑1 (MQ=255)
                            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTgagaa                                                         >  1:1460628/1‑56 (MQ=255)
                            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTgagaa                                                         >  1:1372611/1‑56 (MQ=255)
                            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCActaatg                                               <  1:4527043/67‑1 (MQ=255)
                            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                                            <  1:3093088/70‑3 (MQ=255)
                            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                                            <  1:4734283/70‑3 (MQ=255)
                            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                           >  1:3556466/1‑71 (MQ=255)
                            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                           >  1:2753796/1‑71 (MQ=255)
                             ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCa                                                     <  1:1107040/60‑1 (MQ=255)
                             ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                                   <  1:5173566/62‑1 (MQ=255)
                               ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                                                <  1:5399116/63‑1 (MQ=255)
                               ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                 >  1:4112445/1‑62 (MQ=255)
                               ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAg                                        >  1:4664338/1‑71 (MQ=255)
                               ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAg                                        >  1:2577682/1‑71 (MQ=255)
                               ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                                         <  1:947627/70‑1 (MQ=255)
                               ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                                         <  1:1182047/70‑1 (MQ=255)
                                gATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgt                                       >  1:383292/1‑71 (MQ=255)
                                gATTTCTTATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGaat                                       >  1:3122013/1‑69 (MQ=255)
                                  ttttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATg                                          <  1:2994731/66‑1 (MQ=255)
                                  tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                                     >  1:1970997/1‑55 (MQ=255)
                                  tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                                     >  1:5134465/1‑55 (MQ=255)
                                   tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                               >  1:2172702/1‑60 (MQ=255)
                                   tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                               >  1:3515198/1‑60 (MQ=255)
                                    ttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgttgtt                                   >  1:622813/1‑71 (MQ=255)
                                       tATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcc                                 <  1:2406726/70‑1 (MQ=255)
                                         tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcctcc                              >  1:2529112/1‑71 (MQ=255)
                                           cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTcctcct                             <  1:5342248/70‑1 (MQ=255)
                                           cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                            >  1:3330444/1‑71 (MQ=255)
                                             tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                            <  1:5372862/69‑1 (MQ=255)
                                             tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                            <  1:3229031/69‑1 (MQ=255)
                                               ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                        >  1:1955615/1‑71 (MQ=255)
                                               ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                        >  1:1260997/1‑71 (MQ=255)
                                               ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGAATGATGATATGAGTTGTTCCTCCTGCGACg                       >  1:1625714/1‑71 (MQ=255)
                                                  ttCTGCCCCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                              >  1:4292406/1‑46 (MQ=255)
                                                   tCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCaa                                                    >  1:2466197/1‑39 (MQ=255)
                                                   tCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCa                                                     >  1:813073/1‑38 (MQ=255)
                                                    cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc                   >  1:5671987/1‑71 (MQ=255)
                                                    cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc                   >  1:423132/1‑71 (MQ=255)
                                                    cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGct                    >  1:3810355/1‑70 (MQ=255)
                                                     tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                                   <  1:3935302/38‑1 (MQ=255)
                                                     tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACgg                      >  1:5083692/1‑67 (MQ=255)
                                                       ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                 >  1:4033043/1‑38 (MQ=255)
                                                       ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                 >  1:2298830/1‑38 (MQ=255)
                                                       ccGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctct                  <  1:482278/69‑1 (MQ=255)
                                                       ccGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctct                  <  1:5208430/69‑1 (MQ=255)
                                                          cAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa             <  1:4786608/71‑1 (MQ=255)
                                                           aGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                        <  1:3721742/59‑1 (MQ=255)
                                                             gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                                      <  1:2142377/43‑1 (MQ=255)
                                                             gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                                      <  1:3708620/43‑1 (MQ=255)
                                                             gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                                      <  1:4446608/43‑1 (MQ=255)
                                                             gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa          <  1:5505732/71‑1 (MQ=255)
                                                             gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa          <  1:3290580/71‑1 (MQ=255)
                                                              aTGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                                      <  1:4218532/42‑1 (MQ=255)
                                                               tGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTaa                >  1:5210211/1‑63 (MQ=255)
                                                               tGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa          >  1:2546351/1‑69 (MQ=255)
                                                               tGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa             >  1:4208551/1‑66 (MQ=255)
                                                               tGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa             >  1:711638/1‑66 (MQ=255)
                                                               tGAGACGATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTaa                >  1:1459/1‑63 (MQ=255)
                                                                   aCTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaaaa        >  1:4354855/1‑67 (MQ=255)
                                                                     tataCATCGTCGAGATCGCAATGATGTTATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAg  >  1:3801726/1‑71 (MQ=255)
                                                                     tataCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaaaa        >  1:171021/1‑65 (MQ=255)
                                                                      |                                                                     
TTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAG  >  Exported/2571334‑2571473

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