Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,572,726 A→G 10.2% intergenic (+1494/‑285) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,572,7260AG10.2% 417.3 / 28.6 157intergenic (+1494/‑285)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (79/62);  new base G (10/6);  total (89/68)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 163 total aligned reads displayed.

GATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTC  >  Exported/2572659‑2572794
                                                                   |                                                                    
gATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGcc                                                                   >  1:1086747/1‑71 (MQ=255)
 aTCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCt                                                                  <  1:423499/71‑1 (MQ=255)
        ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                            <  1:1882051/70‑1 (MQ=255)
         cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                          >  1:28989/1‑71 (MQ=255)
           gagaTCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGatat                                                         <  1:3519381/67‑1 (MQ=255)
           gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtg                                                             <  1:3999812/66‑1 (MQ=255)
            aaaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                                 >  1:2153154/1‑61 (MQ=255)
              aTCTGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgt                                                              <  1:877913/62‑1 (MQ=255)
              aTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgt                                                              <  1:4897259/62‑1 (MQ=255)
              aTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCg                                                     <  1:1734648/71‑1 (MQ=255)
               tCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGatc                                                                    <  1:2424029/55‑3 (MQ=255)
               ccGGATTCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtga                                                            <  1:4494340/62‑2 (MQ=255)
                   attCAGCGGATACGCATACAGACTCCATCCCCCCTCCAGCCCTATCGGGTCCTGAGTGATATACCGCCCct                                                <  1:2205840/70‑1 (MQ=255)
                   attCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGat                                                                     >  1:320179/2‑49 (MQ=255)
                    ttCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCTCCCctg                                               >  1:4336450/1‑71 (MQ=255)
                    ttCAACGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATAcc                                                      >  1:1843433/1‑64 (MQ=255)
                    ttCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAgt                                                              <  1:5182467/56‑1 (MQ=255)
                    ttCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                                               >  1:1136220/1‑69 (MQ=255)
                    ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                          <  1:4411233/60‑1 (MQ=255)
                      cAACGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGCCCctgc                                              <  1:5249456/70‑1 (MQ=255)
                      cAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                                               >  1:3173499/1‑67 (MQ=255)
                      cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                          >  1:1947126/1‑58 (MQ=255)
                      cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                                             >  1:3994771/1‑71 (MQ=255)
                       aaCGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTg                                                                 >  1:3002682/1‑50 (MQ=255)
                       aaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                                            >  1:3513883/1‑71 (MQ=255)
                        aCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagcag                                            >  1:3815513/1‑68 (MQ=255)
                        aCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                       >  1:610003/1‑59 (MQ=255)
                         cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                            <  1:1588314/53‑1 (MQ=255)
                          ggATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtga                                                            <  1:5222836/52‑2 (MQ=255)
                          ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                            <  1:1911403/52‑1 (MQ=255)
                          ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                            <  1:241324/52‑1 (MQ=255)
                          ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                       >  1:5633958/1‑57 (MQ=255)
                          ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                       >  1:1156018/1‑57 (MQ=255)
                          ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                       >  1:666565/1‑57 (MQ=255)
                           gATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtgat                                                           >  1:1654182/1‑50 (MQ=255)
                           gATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGAta                                                          >  1:2706125/1‑53 (MQ=255)
                            aTACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGa                                                                      >  1:3232991/1‑40 (MQ=255)
                             tACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGATATATCGCCCCTGCAGCgggtc                                      <  1:2850444/71‑4 (MQ=255)
                             tACTGATAAAAATTCCATCCCCACTTCAGCCCAATCGGATCCTGAg                                                               >  1:4576623/1‑46 (MQ=255)
                             tACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGa                                                                      <  1:1103184/39‑1 (MQ=255)
                              aCTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGa                                                                >  1:240013/1‑44 (MQ=255)
                               cTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGat                                                                     <  1:1049524/38‑2 (MQ=16)
                               cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                          >  1:76943/1‑49 (MQ=255)
                                tGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                                  <  1:3630668/40‑1 (MQ=255)
                                tGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                                  <  1:2519766/40‑1 (MQ=255)
                                tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                                       <  1:4968894/67‑1 (MQ=255)
                                 gATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                                  <  1:5332975/39‑1 (MQ=16)
                                 gATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                                  <  1:1263718/39‑1 (MQ=16)
                                  aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtatat                                                       >  1:1130632/1‑48 (MQ=255)
                                  aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtatat                                                       >  1:5594015/1‑48 (MQ=255)
                                  aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtata                                                        >  1:3729263/1‑48 (MQ=255)
                                  aTAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGa                                                                >  1:4320881/1‑40 (MQ=255)
                                  aTAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGAta                                                          >  1:2750914/1‑46 (MQ=255)
                                  aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                       >  1:4383217/1‑49 (MQ=255)
                                  aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGtt                                 <  1:1445588/71‑1 (MQ=255)
                                   tAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCagcagc                                           <  1:1807372/60‑1 (MQ=255)
                                    aaaaaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGatat                                                         >  1:4746828/1‑45 (MQ=255)
                                    aaaaaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTAc                               >  1:1167399/1‑71 (MQ=255)
                                     aaaaTTCCATCACCCCTTCAGCCCAATCGGATCCTGAGTGAtata                                                        >  1:5468939/1‑45 (MQ=255)
                                       aaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGAta                                                          >  1:1699164/1‑41 (MQ=255)
                                       aaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAcctc                                                    <  1:4132582/47‑3 (MQ=255)
                                            cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                        <  1:967637/70‑1 (MQ=255)
                                            cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                        <  1:2640986/70‑1 (MQ=255)
                                             acctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                                           <  1:514564/50‑1 (MQ=255)
                                             acctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                                           <  1:309814/50‑1 (MQ=255)
                                              cctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt                     <  1:2116581/71‑1 (MQ=255)
                                               ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                                           <  1:399991/48‑1 (MQ=255)
                                               ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCgg                                         >  1:2143794/1‑50 (MQ=255)
                                               ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt                     >  1:2631853/1‑70 (MQ=255)
                                                tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                              <  1:5145287/44‑1 (MQ=255)
                                                tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                              <  1:158107/44‑1 (MQ=255)
                                                tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                              <  1:3548074/44‑1 (MQ=255)
                                                 acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                   >  1:2439551/1‑70 (MQ=255)
                                                  cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGa                            <  1:3796863/60‑1 (MQ=255)
                                                   ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                >  1:1943673/1‑71 (MQ=255)
                                                   ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                <  1:5513053/71‑1 (MQ=255)
                                                    ttCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTa                                >  1:815782/1‑54 (MQ=255)
                                                     tCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                >  1:1695969/1‑69 (MQ=255)
                                                      tcAGCCCAATCGGATCCTGAGTGATATACCGcccc                                                 <  1:4242852/34‑1 (MQ=37)
                                                      tcAGCCCAATCGGATCCTGAGTGATATACCGcccc                                                 <  1:3332834/34‑1 (MQ=37)
                                                      ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                              >  1:274979/1‑38 (MQ=255)
                                                      ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTTCTGCGGTTAATGTTACCGATAATAGTTATAGt               <  1:5113808/69‑1 (MQ=255)
                                                      ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa              <  1:5045414/70‑1 (MQ=255)
                                                       cAGCCCAATCGGATCCTGAGTGATATACCGCCCCTGCAGCgg                                         >  1:701797/1‑42 (MQ=255)
                                                       cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa              >  1:4094487/1‑69 (MQ=255)
                                                       cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa              >  1:2144679/1‑69 (MQ=255)
                                                          ccccATCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc         >  1:2273757/1‑71 (MQ=255)
                                                          ccccATCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc         <  1:1764035/71‑1 (MQ=255)
                                                           cccaTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGccc        >  1:1181916/4‑71 (MQ=255)
                                                           ccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGccc        >  1:2761371/1‑71 (MQ=255)
                                                            cAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata                          >  1:3425088/1‑52 (MQ=255)
                                                            cAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt       <  1:28268/71‑1 (MQ=255)
                                                             caTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTg      <  1:2543541/70‑1 (MQ=255)
                                                             aaTCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                                     <  1:4299080/40‑1 (MQ=255)
                                                              aTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGa     >  1:3363095/1‑71 (MQ=255)
                                                              aTCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGtt                                 >  1:2463427/1‑43 (MQ=255)
                                                               tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAt    >  1:4835131/1‑71 (MQ=255)
                                                               tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAt    >  1:2551937/1‑71 (MQ=255)
                                                                 ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataata                       >  1:4172927/1‑50 (MQ=255)
                                                                 ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt   >  1:5074099/1‑70 (MQ=255)
                                                                   |                                                                    
GATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTC  >  Exported/2572659‑2572794

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: