Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,572,795 A→C 11.3% intergenic (+1563/‑216) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,572,7950AC11.3% 139.4 / 13.6 62intergenic (+1563/‑216)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/25);  new base C (4/3);  total (34/28)
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

GCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCTTCTGTGTCTCTGTATCTGGTTA  >  Exported/2572727‑2572857
                                                                    |                                                              
tccTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc                                                              <  1:2827269/70‑1 (MQ=255)
tccTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc                                                              <  1:5675327/70‑1 (MQ=255)
 ccTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatca                                                             >  1:1458927/1‑71 (MQ=255)
  cTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                                            <  1:1215386/71‑1 (MQ=255)
   tGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                                            >  1:4255053/1‑70 (MQ=255)
   tGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctcat                                                            >  1:4848809/1‑70 (MQ=255)
    gggtgaTATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACtcctcatca                                                          <  1:3320751/69‑1 (MQ=255)
     ggtgaTATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACTcctcatcat                                                         <  1:3500579/70‑1 (MQ=255)
     aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                                          <  1:4024081/70‑1 (MQ=255)
     aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                                          >  1:4364238/1‑70 (MQ=255)
     aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTcatcatcat                                                         <  1:9140/71‑1 (MQ=255)
      gtgaTATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACtcctcat                                                            <  1:4888930/67‑1 (MQ=255)
       taatataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctcat                                                            <  1:1756528/63‑1 (MQ=255)
        gtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                                            <  1:3995348/65‑1 (MQ=255)
        gtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                                            <  1:3750460/65‑1 (MQ=255)
        aatataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATTCCTCATCAtgct                                                      >  1:3625482/3‑71 (MQ=255)
         ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTcatcatcat                                                         <  1:1729571/67‑1 (MQ=255)
         ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                                      <  1:704595/70‑1 (MQ=255)
         ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                                      <  1:1963142/70‑1 (MQ=255)
         atataCCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACTCCTCATCAtgctg                                                     <  1:1260426/70‑2 (MQ=255)
         atataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc                                                              >  1:2377207/2‑62 (MQ=255)
         atataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATTCCTCATCAtgctg                                                     >  1:5647503/2‑70 (MQ=255)
           ataCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATagtact                                                   >  1:311533/1‑71 (MQ=255)
            taCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATTCCTCATCATGctgctg                                                  >  1:1084867/1‑71 (MQ=255)
             aCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                                            <  1:3216827/60‑1 (MQ=255)
              ccGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTActgct                                                <  1:3679169/71‑1 (MQ=255)
                 ccccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCg                                             >  1:1675088/1‑71 (MQ=255)
                   ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                                      <  1:4925133/60‑1 (MQ=255)
                   ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATactactact                                                <  1:2410494/66‑1 (MQ=255)
                   ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCgg                                            >  1:1158939/1‑70 (MQ=255)
                   ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGc                                           <  1:527883/71‑1 (MQ=255)
                    ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                                          <  1:3793395/55‑1 (MQ=255)
                    ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                                          <  1:878015/55‑1 (MQ=255)
                    ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                                          <  1:2079703/55‑1 (MQ=255)
                    ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                                          <  1:1131489/55‑1 (MQ=255)
                    ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCa                                          <  1:3059692/71‑1 (MQ=255)
                     tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                                          <  1:3481780/54‑1 (MQ=255)
                     tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                                          <  1:1512071/54‑1 (MQ=255)
                     tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                                          <  1:303116/54‑1 (MQ=255)
                     tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                                          <  1:4503456/54‑1 (MQ=255)
                     tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                                          <  1:4509804/54‑1 (MQ=255)
                     tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                                          <  1:105533/54‑1 (MQ=255)
                     tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAg                                         >  1:3132992/1‑71 (MQ=255)
                      gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                                            <  1:946158/51‑1 (MQ=255)
                      gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                                            <  1:1960161/51‑1 (MQ=255)
                      gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGc                                           >  1:3048195/1‑68 (MQ=255)
                       ctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGt                                        <  1:4689544/70‑1 (MQ=255)
                       ctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTc                                       >  1:67307/1‑71 (MQ=255)
                         gcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGc                                           >  1:5064551/1‑65 (MQ=255)
                          cGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtg                                     <  1:4764864/70‑1 (MQ=255)
                            gTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtgatg                                  <  1:3982140/71‑1 (MQ=255)
                             ttAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtgat                                   <  1:5478436/69‑1 (MQ=255)
                             ttAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtgat                                   <  1:2148182/69‑1 (MQ=255)
                             ttAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtgat                                   <  1:3604419/69‑1 (MQ=255)
                              tAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGaa                                >  1:587680/1‑71 (MQ=255)
                              tAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGaa                                >  1:610744/1‑71 (MQ=255)
                               aaTGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATactact                                                   >  1:2881159/1‑51 (MQ=255)
                               aaTGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATactact                                                   >  1:5012677/1‑51 (MQ=255)
                               aaTGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGaa                                >  1:1630218/1‑70 (MQ=255)
                                 tGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtga                                    >  1:1852637/1‑64 (MQ=255)
                                 tGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtg                                     <  1:1380358/63‑1 (MQ=255)
                                    tACCGATAATAGTTATAGTACAGCCCTTATTCATCATCAta                                                        >  1:5017773/1‑41 (MQ=255)
                                    tACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta                                                        >  1:379287/1‑41 (MQ=255)
                                    tACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta                                                        >  1:2378101/1‑41 (MQ=255)
                                    tACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta                                                        >  1:143625/1‑41 (MQ=255)
                                    tACCGATAATAGTTATAGTACAGCCCTGATTAATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTg                          >  1:3958912/1‑71 (MQ=255)
                                     aCCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtga                                    >  1:86042/1‑60 (MQ=255)
                                      ccGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                                      <  1:4970106/41‑1 (MQ=255)
                                      ccGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTa                             >  1:1186802/1‑66 (MQ=255)
                                      ccGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCt                        >  1:678455/1‑71 (MQ=255)
                                        gATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCttt                      <  1:2099569/71‑2 (MQ=255)
                                          taataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCtttta                    >  1:548902/1‑68 (MQ=255)
                                          taataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCtttta                    >  1:5567465/1‑68 (MQ=255)
                                          taataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCtttta                    <  1:3806569/71‑4 (MQ=255)
                                            ataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCg                                             >  1:2202934/1‑44 (MQ=255)
                                            ataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCg                                             >  1:2367906/1‑44 (MQ=255)
                                            ataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAg                                         <  1:4014853/48‑1 (MQ=255)
                                            ataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAg                                         <  1:1650515/48‑1 (MQ=255)
                                            ataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAg                                         <  1:2247157/48‑1 (MQ=255)
                                              ggTTATAGTACAGCCCTGATTCCTCATCATgctgct                                                   >  1:4151177/2‑36 (MQ=255)
                                                  ataGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCTTCTGTGTCTCTg            >  1:3680414/1‑71 (MQ=255)
                                                     atACAGGCCGGACTCCTCATCATACTGCTGACCCGGCAgg                                        <  1:2486317/39‑2 (MQ=255)
                                                      tACAGGCCGGACTCCTCATCATACTGCTGCCCCGGCAGCCggata                                  <  1:5636144/45‑2 (MQ=255)
                                                      tACAGCCCTGATTCATCATCATACTACTGCTGCGGCAg                                         <  1:3016099/38‑1 (MQ=255)
                                                      tACAGCCCTGATTCATCATCATACTACTGCTGCGGCAg                                         <  1:920544/38‑1 (MQ=255)
                                                       aCAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCtttt                     <  1:4272338/57‑3 (MQ=255)
                                                       aCAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCtttt                     <  1:1686976/57‑3 (MQ=255)
                                                       aCAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCtttt                     <  1:2640343/57‑3 (MQ=255)
                                                        cAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTAtt                           >  1:344311/1‑50 (MQ=255)
                                                        cAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTAtt                           >  1:2445121/1‑50 (MQ=255)
                                                        cAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCTTCTGTGTCTCTGTATCTg      >  1:2969893/1‑71 (MQ=255)
                                                             cTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCTTCTGTGTCTCTGTATCTGGTTa  >  1:1857810/1‑70 (MQ=255)
                                                                    |                                                              
GCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCTTCTGTGTCTCTGTATCTGGTTA  >  Exported/2572727‑2572857

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: