Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 234,705 T→A 10.7% intergenic (+14/‑169) yaiI → / → aroL conserved hypothetical protein/shikimate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported234,7050TA10.7% 152.6 / 10.0 56intergenic (+14/‑169)yaiI/aroLconserved hypothetical protein/shikimate kinase II
Reads supporting (aligned to +/- strand):  ref base T (22/28);  new base A (2/4);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CCTTTGCCGCGGAGCTGGAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTGAGATTTTCACTTTAAGTGGAA  >  Exported/234636‑234773
                                                                     |                                                                    
ccTTTGCCGCGGAGCTGGAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTattatt                                                                     <  1:2410543/71‑1 (MQ=255)
    tGCCGCGGAGCTGGAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTac                                                                  <  1:3115352/70‑1 (MQ=255)
    tGCCGCGGAGCTGGAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTa                                                                   >  1:1422804/1‑69 (MQ=255)
     gCCGCGGAGCTGGAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTacac                                                                >  1:1006271/1‑71 (MQ=255)
      ccGCGGAGCTGGAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAAATAATAAATTacatt                                                               <  1:364988/71‑3 (MQ=38)
            aGCTGGAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTc                                                         <  1:1128180/71‑1 (MQ=255)
            aGCTGGAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTc                                                         <  1:2909705/71‑1 (MQ=255)
            aGCTGGAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTc                                                         <  1:1583035/71‑1 (MQ=255)
            aGCTGCAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTc                                                         <  1:1493151/71‑1 (MQ=255)
                ggAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTtattat                                                                      <  1:3644695/54‑1 (MQ=255)
                ggAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTg                                                      >  1:1196223/1‑70 (MQ=255)
                 gagaAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTg                                                      <  1:424835/69‑1 (MQ=255)
                  agaAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTCACACTTCATTCTTGaa                                                    >  1:103587/1‑70 (MQ=255)
                   gaAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAata                                                  >  1:3196971/1‑71 (MQ=255)
                   gaAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAata                                                  >  1:3151747/1‑71 (MQ=255)
                   gaAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAata                                                  >  1:2108233/1‑71 (MQ=255)
                     aGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAatat                                                 <  1:2826984/70‑1 (MQ=255)
                     aGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAAtatt                                                <  1:397570/71‑1 (MQ=255)
                         gtggCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAAtatttat                                             <  1:408897/70‑1 (MQ=255)
                          tggCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTattatt                                                                     <  1:145996/45‑1 (MQ=255)
                          tggCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTg                                           >  1:2872154/1‑71 (MQ=255)
                           ggCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTgg                                          >  1:3666513/1‑71 (MQ=255)
                           ggCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTgg                                          >  1:3060414/1‑71 (MQ=255)
                               ggAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGta                                        >  1:2583782/1‑69 (MQ=255)
                                gAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTg                                           >  1:2754735/1‑65 (MQ=255)
                                 aaGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAg                                     >  1:2190826/1‑70 (MQ=255)
                                  aGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTa                                   <  1:1911789/71‑1 (MQ=255)
                                  aGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTa                                   <  1:1210369/71‑1 (MQ=255)
                                  aGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTa                                   <  1:258144/71‑1 (MQ=255)
                                   gTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTaa                                  <  1:2514404/71‑1 (MQ=255)
                                    tGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACtt                                                              <  1:1672050/42‑1 (MQ=255)
                                    tGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACtt                                                              <  1:2891856/42‑1 (MQ=255)
                                    tGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAAtatttatt                                            >  1:387138/1‑60 (MQ=255)
                                      cAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAggg                               <  1:3095357/71‑1 (MQ=255)
                                      cAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAggg                               <  1:827348/71‑1 (MQ=255)
                                      cAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAggg                               >  1:3512008/1‑71 (MQ=255)
                                      cAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAggg                               >  1:1832722/1‑71 (MQ=255)
                                        aCGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGgt                             >  1:360836/1‑71 (MQ=255)
                                         cGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGgtg                            >  1:2420606/1‑71 (MQ=255)
                                             gTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAAtatttatt                                            >  1:2491591/1‑51 (MQ=255)
                                              tCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAg                                     <  1:2558154/57‑1 (MQ=255)
                                              tCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTg                       <  1:3174232/71‑1 (MQ=255)
                                              tCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTCAATATTTATTGGTATAg                                     <  1:2562433/57‑1 (MQ=255)
                                                 tGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTgaga                    <  1:3101280/71‑1 (MQ=255)
                                                 tGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTgaga                    <  1:2997575/71‑1 (MQ=255)
                                                 tGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTgaga                    <  1:1533548/71‑1 (MQ=255)
                                                   gCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTGAGAt                   <  1:1827138/70‑1 (MQ=255)
                                                      aaaTGTAATTTATTATTTACACTTCATTCTTGAATATTTATTgg                                          <  1:2763342/44‑1 (MQ=255)
                                                       agtGTAAATAATAAATTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTGAGAtttt                <  1:1982645/67‑1 (MQ=38)
                                                       agtGTAAATAATAAATTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTGAGAtttt                <  1:1713497/67‑1 (MQ=38)
                                                       agtGTAAATAATAAATTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTGAGAtttt                <  1:1149981/67‑1 (MQ=38)
                                                       aaTGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTGAGATTTTCa              >  1:821139/1‑71 (MQ=255)
                                                        gtGTAAATAATAAATTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTGAGATTTTCAc             >  1:1909310/2‑71 (MQ=255)
                                                        gtGTAAATAATAAATTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTGAGATTTTCAc             >  1:1768834/2‑71 (MQ=255)
                                                          gTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGt                                    <  1:2519966/46‑1 (MQ=255)
                                                          gTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGt                                    <  1:1454980/46‑1 (MQ=255)
                                                          gTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGt                                    <  1:1375425/46‑1 (MQ=255)
                                                            aaTTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTGAGATTTTCACttt          >  1:1109020/1‑70 (MQ=255)
                                                             aTTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTGAGATTTTCACttt          >  1:3453155/1‑69 (MQ=255)
                                                                   tattTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTGAGATTTTCACTTTAAGTGGaa  >  1:1188195/1‑71 (MQ=255)
                                                                     |                                                                    
CCTTTGCCGCGGAGCTGGAGAAGTGGTGGCTGGAAGTGCAACGTAGTCGTGGCTAAATGTAATTTATTATTTACACTTCATTCTTGAATATTTATTGGTATAGTAAGGGGTGTATTGAGATTTTCACTTTAAGTGGAA  >  Exported/234636‑234773

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: