Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,572,487 A→G 13.3% intergenic (+1255/‑524) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,572,4870AG13.3% 200.3 / 20.4 75intergenic (+1255/‑524)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base G (4/6);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

TATGCCGTCAGGTTTTCTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTG  >  Exported/2572417‑2572555
                                                                      |                                                                    
tATGCCGTCAGGTTTTCTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTa                                                                      <  1:1706495/71‑1 (MQ=255)
  tGCCGTCAGGTTTTCTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAg                                                                     <  1:2969598/70‑1 (MQ=255)
  tGCCGTCAGGTTTTCTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGc                                                                    <  1:2118040/71‑1 (MQ=255)
            ttttCTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTta                                                          <  1:2094796/71‑1 (MQ=255)
            ttttCTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTta                                                          >  1:695270/1‑71 (MQ=255)
            ttttCTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTta                                                          <  1:1980591/71‑1 (MQ=255)
            ttttCTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTta                                                          <  1:1711207/71‑1 (MQ=255)
             tttCTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTta                                                          >  1:3178750/1‑70 (MQ=255)
               tCTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTtata                                                        <  1:351166/70‑1 (MQ=255)
               tCTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTtata                                                        <  1:3327674/70‑1 (MQ=255)
                cTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATaaa                                                      >  1:3037807/1‑71 (MQ=255)
                   ttgcttgcTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTtat                                                         <  1:1403811/65‑1 (MQ=255)
                   ttgcttgcTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTtat                                                         <  1:2599088/65‑1 (MQ=255)
                      cttgctCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATaa                                                       <  1:2766502/64‑1 (MQ=255)
                        tgctCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTAc                                               >  1:2076148/1‑70 (MQ=255)
                         gctCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTGGCCATCTATTTATAAATGTGTACaa                                             >  1:1072681/1‑71 (MQ=255)
                         gctCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACa                                              <  1:2887650/70‑1 (MQ=255)
                          ctCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTGGCCATCTATTTATAAATGTGTACaa                                             <  1:2698685/70‑1 (MQ=255)
                          ctCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTGGCCATCTATTTATAAATGTGTACaa                                             <  1:2141018/70‑1 (MQ=255)
                          ctCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTGGCCATCTATTTATAAATGTGTACAAc                                            >  1:1615545/1‑71 (MQ=255)
                          ctCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAAc                                            >  1:1835470/1‑71 (MQ=255)
                          ctCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAAc                                            >  1:101386/1‑71 (MQ=255)
                           tCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACg                                           >  1:2050285/1‑71 (MQ=255)
                             ccGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACg                                           <  1:1744106/69‑1 (MQ=255)
                              cGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACggg                                         >  1:903323/1‑70 (MQ=255)
                               gggTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTGGCCATCTATTTATAAATGTGTACAACGGGGa                                       <  1:163864/71‑1 (MQ=255)
                               gggTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAg                                                                     >  1:1425145/1‑41 (MQ=255)
                               gggTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGa                                       <  1:1542032/71‑1 (MQ=255)
                               gggTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGa                                       <  1:613559/71‑1 (MQ=255)
                                ggTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGat                                      <  1:2577989/71‑1 (MQ=255)
                                  tCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGatt                                     <  1:63402/70‑1 (MQ=255)
                                  tCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGatt                                     <  1:2633292/70‑1 (MQ=255)
                                    aTTAACATTCCAACTCCATAGGCCTTGAGGTATTGGCCATCTATTTATAAATGTGTACAACGGGGattat                                   <  1:869440/70‑1 (MQ=255)
                                    aTTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGattat                                   <  1:3670738/70‑1 (MQ=255)
                                    aTTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGattat                                   <  1:2667569/70‑1 (MQ=255)
                                     ttAACATTCCAACTCCATAGGCCTTGAGGTATTGGCCATCTAt                                                             <  1:2864669/43‑1 (MQ=255)
                                     ttAACATTCCAACTCCATAGGCCTTGAGGTATTGGCCATCTAt                                                             <  1:1820636/43‑1 (MQ=255)
                                     ttAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTAttt                                 >  1:3190274/1‑71 (MQ=255)
                                     ttAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTAttt                                 >  1:3377502/1‑71 (MQ=255)
                                       aaCATTCCAACTCCATAGGCCTTGAGGTCTTAGCCATCTATTTATAAATGTGTACAACGGGGATTAtttt                                <  1:1823261/70‑1 (MQ=255)
                                        aCATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTa                              >  1:760695/1‑71 (MQ=255)
                                         cATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTa                              >  1:79855/1‑70 (MQ=255)
                                            tCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTAc                                               >  1:612440/1‑50 (MQ=255)
                                            tCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTAc                                               >  1:128734/1‑50 (MQ=255)
                                            tCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTAc                                               >  1:2270410/1‑50 (MQ=255)
                                             ccAACTCCATAGGCCTTGAGGTATTAGCCATCTATTta                                                          >  1:2433082/1‑38 (MQ=255)
                                             ccAACTCCATAGGCCTTGAGGTATTAGCCATCTATTta                                                          >  1:2307531/1‑38 (MQ=255)
                                             ccAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATaa                                                       <  1:3692908/41‑1 (MQ=255)
                                              cAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTg                         <  1:320315/70‑1 (MQ=255)
                                              cAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGc                        <  1:3070558/71‑1 (MQ=255)
                                               aaCTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGattatt                                  <  1:251878/60‑1 (MQ=255)
                                                 cTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCata                     <  1:1076396/71‑1 (MQ=255)
                                                 cTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCata                     <  1:2590508/71‑1 (MQ=255)
                                                 cTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCat                      >  1:3129248/1‑70 (MQ=255)
                                                   ccATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACgg                                          >  1:1425500/1‑48 (MQ=255)
                                                   ccATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGatta                                    <  1:1702629/54‑1 (MQ=255)
                                                    cATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTa                              <  1:2447418/59‑1 (MQ=255)
                                                      tAGGCCTTGAGGTATTGGCCATCTATTTATAAATGTGTACAACGGGGattat                                   >  1:3372007/1‑52 (MQ=255)
                                                      tAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACgggg                                        >  1:3582783/1‑47 (MQ=255)
                                                      tAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCAtattat                 >  1:2027299/1‑70 (MQ=255)
                                                        ggCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCata                     >  1:304899/1‑64 (MQ=255)
                                                         gCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTAc              >  1:452515/1‑70 (MQ=255)
                                                          ccTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTAc              >  1:821776/1‑69 (MQ=255)
                                                           cTTGAGGTATTGGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTAc              >  1:2858601/1‑68 (MQ=255)
                                                           cTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAg           <  1:1201218/71‑1 (MQ=255)
                                                             tGAGGTATTAGCCATCTATTTATAAATGTGTACAACg                                           >  1:2885249/1‑37 (MQ=255)
                                                             tGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAg           >  1:3199139/1‑69 (MQ=255)
                                                             tGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAg           >  1:30467/1‑69 (MQ=255)
                                                             tGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAg           >  1:728570/1‑69 (MQ=255)
                                                              gAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACAt                           >  1:290535/1‑52 (MQ=255)
                                                              gAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGaa         <  1:2409516/70‑1 (MQ=255)
                                                               aGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGc       >  1:2386953/1‑71 (MQ=255)
                                                                 gTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGa     >  1:2896732/1‑71 (MQ=255)
                                                                   aTTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACt   >  1:799605/1‑71 (MQ=255)
                                                                    ttAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGa          <  1:745963/63‑1 (MQ=255)
                                                                     tAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTg  >  1:1452331/1‑70 (MQ=255)
                                                                      aGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCAtattat                 >  1:1399105/1‑54 (MQ=255)
                                                                      aGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCg      <  1:2629929/65‑1 (MQ=255)
                                                                      aGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTg  <  1:833357/69‑1 (MQ=255)
                                                                      |                                                                    
TATGCCGTCAGGTTTTCTATTGCTTGCTCCCGGGTCATTAACATTCCAACTCCATAGGCCTTGAGGTATTAGCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTG  >  Exported/2572417‑2572555

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: