Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,571,425 G→A 9.7% intergenic (+193/‑1586) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,571,4250GA 9.7% 252.2 / 12.1 93intergenic (+193/‑1586)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (52/31);  new base A (6/3);  total (58/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

TGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGCTTTTCCATCTATTGC  >  Exported/2571355‑2571494
                                                                      |                                                                     
tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                                                       <  1:3810209/71‑1 (MQ=255)
tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                                                       <  1:1268388/71‑1 (MQ=255)
tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                                                       <  1:4960781/71‑1 (MQ=255)
    ctcATTAGATTTCTCATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatga                                                                   >  1:2317238/1‑71 (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)
       aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                                   <  1:3581123/68‑1 (MQ=255)
       aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                                                                <  1:3135749/71‑1 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCActaatg                                                                    <  1:4527043/67‑1 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                                                                 <  1:4734283/70‑3 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                                                >  1:1508719/1‑71 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                                                >  1:2753796/1‑71 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                                                >  1:1715371/1‑71 (MQ=255)
       aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                                                >  1:3556466/1‑71 (MQ=255)
         tGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                                                                >  1:5176006/1‑69 (MQ=255)
          ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                                                                     <  1:5399116/63‑1 (MQ=255)
          ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAg                                                             >  1:4664338/1‑71 (MQ=255)
          ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAg                                                             >  1:2577682/1‑71 (MQ=255)
          ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTTTGCATTGTTGAGAACGCACTAATGACATGa                                                              <  1:5040939/70‑1 (MQ=255)
          ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                                                              <  1:1182047/70‑1 (MQ=255)
          ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                                                              <  1:947627/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)
              tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                                                    >  1:3515198/1‑60 (MQ=255)
              tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                                                    >  1:2172702/1‑60 (MQ=255)
               ttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgttgtt                                                        >  1:622813/1‑71 (MQ=255)
                tttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtc                                                       >  1:2823603/1‑71 (MQ=255)
                  tATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcc                                                      <  1:2406726/70‑1 (MQ=255)
                    tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                                   <  1:1538276/55‑1 (MQ=255)
                    tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                                   <  1:5384567/55‑1 (MQ=255)
                    tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcctcc                                                   >  1:2529112/1‑71 (MQ=255)
                      cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTcctcct                                                  <  1:5342248/70‑1 (MQ=255)
                      cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                                                 >  1:5491644/1‑71 (MQ=255)
                      cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                                                 >  1:3330444/1‑71 (MQ=255)
                        tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                                                 <  1:3229031/69‑1 (MQ=255)
                        tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                                                 <  1:5372862/69‑1 (MQ=255)
                          ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                                                                >  1:4034113/1‑52 (MQ=255)
                          ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                                             >  1:1955615/1‑71 (MQ=255)
                          ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGAATGATGATATGAGTTGTTCCTCCTGCGACg                                            >  1:1619006/1‑71 (MQ=255)
                          ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGAATGATGATATGAGTTGTTCCTCCTGCGACg                                            >  1:1625714/1‑71 (MQ=255)
                            gTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                                             <  1:1230706/69‑1 (MQ=255)
                             ttCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                                   <  1:4301107/46‑1 (MQ=255)
                             ttCTGCCCCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                                   >  1:4292406/1‑46 (MQ=255)
                               cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                                   <  1:1971366/44‑1 (MQ=255)
                               cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc                                        >  1:4362226/1‑71 (MQ=255)
                               cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc                                        >  1:5671987/1‑71 (MQ=255)
                               cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGct                                         >  1:3810355/1‑70 (MQ=255)
                                  ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                                      >  1:4033043/1‑38 (MQ=255)
                                  ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                                      >  1:2298830/1‑38 (MQ=255)
                                  ccGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctct                                       <  1:5208430/69‑1 (MQ=255)
                                   cGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCg                                               >  1:3029933/1‑60 (MQ=255)
                                     cAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                                  <  1:4786608/71‑1 (MQ=255)
                                      aGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                                             <  1:3721742/59‑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)
                                        gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa                               <  1:258067/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)
                                            agaCTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAag                            <  1:3733606/70‑1 (MQ=255)
                                              aCTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAagag                          >  1:242859/1‑70 (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)
                                                    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:224432/1‑71 (MQ=255)
                                                      tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc                 >  1:3781592/1‑71 (MQ=255)
                                                      tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc                 >  1:1351916/1‑71 (MQ=255)
                                                             gaTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTaa                                     >  1:4148712/1‑44 (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)
                                                              aTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAag                            <  1:4729998/52‑1 (MQ=255)
                                                              aTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAgg                      <  1:84878/58‑1 (MQ=255)
                                                              aTCGCAATCATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa                               <  1:4850145/49‑1 (MQ=255)
                                                                cGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAg                    >  1:5447684/1‑58 (MQ=255)
                                                                cGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGCTTTTCCATc        >  1:2793358/1‑70 (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)
                                                                     tgatgaTATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAg                       <  1:1920484/50‑1 (MQ=255)
                                                                     tgatgaTATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGCTTTTCCATCTATTGc  <  1:1911131/71‑1 (MQ=255)
                                                                     tgatgaTATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGCTTTTCCATCTATTGc  >  1:716569/1‑71 (MQ=255)
                                                                      |                                                                     
TGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGCTTTTCCATCTATTGC  >  Exported/2571355‑2571494

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