Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,042,854 A→G 11.3% intergenic (‑202/‑67) glgS ← / → yqiJ predicted glycogen synthesis protein/predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,042,8540AG11.3% 147.4 / 8.8 53intergenic (‑202/‑67)glgS/yqiJpredicted glycogen synthesis protein/predicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (17/30);  new base G (3/3);  total (20/33)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATAT  >  Exported/2042786‑2042920
                                                                    |                                                                  
aaaTAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATa                                                                  <  1:3417513/71‑1 (MQ=255)
   tAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATaaa                                                                >  1:374216/1‑70 (MQ=255)
     aaaaaCCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATg                                                              <  1:1044298/70‑1 (MQ=255)
      aaaaCCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTa                                                            <  1:807150/71‑1 (MQ=255)
       aaaCCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTaa                                                           <  1:3104638/71‑1 (MQ=255)
       aaaCCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTaa                                                           <  1:1768872/71‑1 (MQ=255)
          ccATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCa                                                        <  1:877496/71‑1 (MQ=255)
          ccATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCa                                                        <  1:201624/71‑1 (MQ=255)
           cATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGAt                                                                   <  1:1640929/59‑1 (MQ=255)
             tatTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAgg                                                     <  1:3465889/71‑1 (MQ=255)
              atTAAAAGATTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATg                                                              >  1:1482031/1‑61 (MQ=255)
              atTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAgg                                                     >  1:1117647/1‑70 (MQ=255)
               ttAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGa                                                    <  1:3370330/70‑1 (MQ=255)
                tAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGaaa                                                  >  1:1551535/1‑71 (MQ=255)
                      aCTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTa                                            <  1:2843360/71‑1 (MQ=255)
                          tctcTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACa                                          >  1:3294860/1‑69 (MQ=255)
                              ttCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGAt                                                                   <  1:1643102/40‑1 (MQ=255)
                              ttCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGt                                      <  1:1618637/69‑1 (MQ=255)
                              ttCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGt                                      <  1:1413788/69‑1 (MQ=255)
                                ccAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAg                                                          >  1:540943/1‑47 (MQ=255)
                                  aGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAatat                                <  1:429296/71‑1 (MQ=255)
                                  aGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAatat                                <  1:280158/71‑1 (MQ=255)
                                     aCTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGc                             <  1:3688374/71‑1 (MQ=255)
                                          tctAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGc                             <  1:1184487/66‑1 (MQ=255)
                                            tAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGt                           >  1:2797279/1‑66 (MQ=255)
                                            tAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACt                       >  1:662451/1‑70 (MQ=255)
                                            tAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACt                       >  1:3545435/1‑70 (MQ=255)
                                             aaCCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTa                      >  1:1093924/1‑70 (MQ=255)
                                              aCCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTaaa                    >  1:829079/1‑71 (MQ=255)
                                              aCCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTaa                     <  1:2529947/70‑1 (MQ=255)
                                              aCCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTaa                     <  1:2417443/70‑1 (MQ=255)
                                                  gtgtTAATAATAACACAGGTAAATGTAAGCAAGGa                                                    <  1:749202/35‑1 (MQ=17)
                                                  gtgtTAATAATAACACAGGTAAATGTAAGCAAGGAAACTATTACAAAGTaaa                                   <  1:2270071/52‑1 (MQ=37)
                                                   tgtTAATAATAACACAGGTAAATGTAAGCAAGGAAACTATTACAAAGTaaa                                   <  1:2635767/51‑1 (MQ=37)
                                                    gtTAATAATAACACAGGTAAATGTAAGCAAGGaaa                                                  >  1:661771/1‑35 (MQ=18)
                                                    gtTAATAATAACACAGGTAAATGTAAGCAAGGaaa                                                  >  1:1419580/1‑35 (MQ=18)
                                                    gtTAATAATAACACAGGTAAATGTAAGCAAGGaaa                                                  >  1:3644511/1‑35 (MQ=18)
                                                     ttattattaACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACAt             <  1:2041158/71‑1 (MQ=255)
                                                     ttattattaACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACAt             <  1:1392834/71‑1 (MQ=255)
                                                     ttattattaACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACAt             >  1:3708411/1‑71 (MQ=255)
                                                          attaACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGa         <  1:2304077/70‑1 (MQ=255)
                                                            taACACAGATAAATGTAAGCAAGGAAACTATTAGAAAGTAAATATTGCGTCACTAAATGGACATTGGAGt       >  1:376208/1‑70 (MQ=255)
                                                             aaCACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGa     >  1:2756441/1‑71 (MQ=255)
                                                               cacaGATAAATGTAAGCAAGGAAACTATTACAAAgg                                      <  1:3010876/36‑2 (MQ=255)
                                                               cacaGATAAATGTAAGCAAGGAAACTATTACAAAgg                                      <  1:2965590/36‑2 (MQ=255)
                                                               cacaGATAAATGTAAGCAAGGAAACTATTACAAAgg                                      <  1:754802/36‑2 (MQ=255)
                                                               cacaGATAAATGTAAGCAAGGAAACTATTACAAAgg                                      <  1:2849037/36‑2 (MQ=255)
                                                                acaGATAAATGTAAGCAAGGAAACTATTACAAAGTAAata                                 >  1:1471127/1‑40 (MQ=255)
                                                                acaGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGata   >  1:1790404/1‑70 (MQ=255)
                                                                 caGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGatat  <  1:485252/70‑1 (MQ=255)
                                                                 caGATAAATGTAAGCAAGGAAACTATTACAAAGGaa                                    >  1:2275299/1‑36 (MQ=255)
                                                                  aGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGat    <  1:1821557/67‑1 (MQ=255)
                                                                  aGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGat    <  1:1344695/67‑1 (MQ=255)
                                                                    |                                                                  
AAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATAT  >  Exported/2042786‑2042920

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: