Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,571,423 A→C 13.2% intergenic (+191/‑1588) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,571,4230AC13.2% 224.9 / 33.4 106intergenic (+191/‑1588)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (53/39);  new base C (8/6);  total (61/45)
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.45e-01

Only 100 of 130 total aligned reads displayed.

TGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGCTTTTCCATC  >  Exported/2571355‑2571488
                                                                    |                                                                 
tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                                                 <  1:1268388/71‑1 (MQ=255)
tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                                                 <  1:3810209/71‑1 (MQ=255)
tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                                                 <  1:4960781/71‑1 (MQ=255)
tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                                                  <  1:4984534/70‑1 (MQ=255)
  ctctCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                                                               <  1:151787/71‑1 (MQ=255)
    ctcATTAGATTTCTCATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatga                                                             >  1:1943869/1‑71 (MQ=255)
     tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             >  1:4314960/1‑70 (MQ=255)
      cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                                              <  1:3904393/68‑1 (MQ=255)
      cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                                              <  1:1937277/68‑1 (MQ=255)
      cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGata                                                           >  1:2811266/1‑71 (MQ=255)
      cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGata                                                           >  1:1951030/1‑71 (MQ=255)
       aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             <  1:552476/68‑1 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCActaatg                                                              <  1:4527043/67‑1 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                                                           <  1:2542554/70‑3 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                                                           <  1:3093088/70‑3 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                                                           <  1:4734283/70‑3 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                                          >  1:1715371/1‑71 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                                          >  1:3556466/1‑71 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                                          >  1:2753796/1‑71 (MQ=255)
        ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                                                  <  1:5173566/62‑1 (MQ=255)
         tGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                                                          >  1:5176006/1‑69 (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)
           gATTTCTTATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGaat                                                      >  1:3122013/1‑69 (MQ=255)
             ttttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATg                                                         <  1:231130/66‑1 (MQ=255)
             ttttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATg                                                         <  1:2994731/66‑1 (MQ=255)
              tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                                              >  1:3515198/1‑60 (MQ=255)
              tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                                              >  1:2172702/1‑60 (MQ=255)
               ttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgttgtt                                                  >  1:622813/1‑71 (MQ=255)
                tttATCGGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCaa                                                                   >  1:731348/1‑53 (MQ=255)
                tttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtc                                                 >  1:2823603/1‑71 (MQ=255)
                  tATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcc                                                <  1:2406726/70‑1 (MQ=255)
                    tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             <  1:5384567/55‑1 (MQ=255)
                    tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             <  1:3553576/55‑1 (MQ=255)
                    tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             <  1:1538276/55‑1 (MQ=255)
                      cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTcctcct                                            <  1:5342248/70‑1 (MQ=255)
                      cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                                           >  1:5491644/1‑71 (MQ=255)
                        tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                                           <  1:3229031/69‑1 (MQ=255)
                          ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                                                          >  1:4034113/1‑52 (MQ=255)
                          ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGAATGATGATATGAGTTGTTCCTCCTGCGACg                                      >  1:1619006/1‑71 (MQ=255)
                            gTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                                       <  1:1230706/69‑1 (MQ=255)
                             ttCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             >  1:3234548/1‑46 (MQ=255)
                             ttCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             <  1:4301107/46‑1 (MQ=255)
                              tCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCaa                                                                   >  1:2466197/1‑39 (MQ=255)
                               cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             <  1:1971366/44‑1 (MQ=255)
                               cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc                                  >  1:423132/1‑71 (MQ=255)
                               cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGct                                   >  1:3810355/1‑70 (MQ=255)
                                tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                                                  <  1:3935302/38‑1 (MQ=255)
                                tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                                                  <  1:64923/38‑1 (MQ=255)
                                tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACgg                                     >  1:5083692/1‑67 (MQ=255)
                                  ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                                >  1:1664118/1‑38 (MQ=255)
                                  ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                                >  1:2298830/1‑38 (MQ=255)
                                     cAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                            <  1:4786608/71‑1 (MQ=255)
                                      aGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                                       <  1:3721742/59‑1 (MQ=255)
                                        gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                                                     <  1:4446608/43‑1 (MQ=255)
                                        gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                                                     <  1:3708620/43‑1 (MQ=255)
                                        gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTa                                <  1:5339379/64‑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)
                                          tGAGACGATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTaa                               >  1:1459/1‑63 (MQ=255)
                                              aCTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaaaa                       >  1:4354855/1‑67 (MQ=255)
                                              aCTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAagag                    >  1:242859/1‑70 (MQ=255)
                                                tataCATCGTCGAGATCGCAATGATGTTATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAg                 >  1:3801726/1‑71 (MQ=255)
                                                tataCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaaaa                       >  1:171021/1‑65 (MQ=255)
                                                  taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCTACGGCTCTAATCAAAAAAGAGCAgg                <  1:3670887/70‑1 (MQ=255)
                                                  taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAAt                              <  1:5045480/56‑1 (MQ=255)
                                                  taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAgg                <  1:2745143/70‑1 (MQ=255)
                                                  taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAg                 >  1:4592541/1‑69 (MQ=255)
                                                  taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAg                 >  1:2680435/1‑69 (MQ=255)
                                                    cATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc                                  >  1:63521/1‑50 (MQ=255)
                                                    cATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                            >  1:4450503/1‑56 (MQ=255)
                                                    cATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                            >  1:3092649/1‑56 (MQ=255)
                                                    cATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGt             >  1:4219500/1‑71 (MQ=255)
                                                     aTCGTCGAGATCGCAATGATGATATGAGTTGTtcct                                               >  1:4080725/1‑36 (MQ=255)
                                                     aTCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAg                 <  1:1990394/66‑1 (MQ=255)
                                                      tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc           >  1:1351916/1‑71 (MQ=255)
                                                      tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc           >  1:4018276/1‑71 (MQ=255)
                                                      tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc           >  1:3781592/1‑71 (MQ=255)
                                                      tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc           >  1:224432/1‑71 (MQ=255)
                                                      tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc           >  1:4684665/1‑71 (MQ=255)
                                                           gagaTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAg              >  1:5419034/1‑63 (MQ=255)
                                                             gaTCGCAATGATGATATGAGTTGTTCCTCCTGCGAAGGCTCTaa                               >  1:1148926/1‑44 (MQ=255)
                                                              aTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa                         <  1:2593939/49‑1 (MQ=255)
                                                              aTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAag                      <  1:124382/52‑1 (MQ=255)
                                                              aTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAag                      <  1:4736522/52‑1 (MQ=255)
                                                              aTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAgg                <  1:84878/58‑1 (MQ=255)
                                                              aTCGCAATCATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa                         <  1:4850145/49‑1 (MQ=255)
                                                                cGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGCTTTTCCATc  >  1:1901843/1‑70 (MQ=255)
                                                                   aaTGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCTGGAGTGCtt         >  1:1967423/1‑60 (MQ=255)
                                                                   aaTGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGCtt         >  1:4438096/1‑60 (MQ=255)
                                                                    atgatgATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa                         <  1:3310542/43‑1 (MQ=255)
                                                                    atgatgATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                            <  1:5639258/40‑1 (MQ=255)
                                                                    |                                                                 
TGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGCTTTTCCATC  >  Exported/2571355‑2571488

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