Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 975,257:1 +A 68.4% coding (358/375 nt) ydeH ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported975,2571.A68.4% 145.1 / 90.9 76coding (358/375 nt)ydeHconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (10/14);  new base A (27/25);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCT  >  Exported/975188‑975319
                                                                      |                                                              
aGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTa                                                               <  1:835760/71‑1 (MQ=255)
      cgacgaCCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAt                                                              <  1:4661400/66‑1 (MQ=255)
      cgacgaCCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGtt                                                         <  1:4493219/71‑1 (MQ=255)
      cgacgaCCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGtt                                                         <  1:143029/71‑1 (MQ=255)
       gacgacCCGGCAATGCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGtt                                                         <  1:18858/71‑1 (MQ=255)
       gacgacCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGt                                                          >  1:2392563/1‑70 (MQ=255)
       gacgacCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGt                                                          >  1:3454128/1‑70 (MQ=255)
               ggCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCg                                                  >  1:4437141/1‑69 (MQ=255)
               ggCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCg                                                  >  1:744262/1‑69 (MQ=255)
                gCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGtt                                                         <  1:656069/62‑1 (MQ=255)
                   aTCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTa                                             <  1:2866552/71‑1 (MQ=255)
                   aTCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTaa                                            <  1:892241/71‑1 (MQ=255)
                    tCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGtt                                                         <  1:5506891/58‑1 (MQ=39)
                    tCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTa                                             >  1:268053/1‑70 (MQ=255)
                     cccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTaa                                            >  1:1746068/1‑70 (MQ=255)
                      ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTaa                                                       >  1:3034630/1‑58 (MQ=39)
                      ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTaa                                                       >  1:4182202/1‑58 (MQ=39)
                      ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTaa                                                       >  1:5108654/1‑58 (MQ=39)
                      ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTaaaa                                          <  1:4072864/71‑1 (MQ=255)
                      ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTaaaa                                          <  1:859248/71‑1 (MQ=255)
                      ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTaa                                            >  1:2911408/1‑69 (MQ=255)
                      ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTa                                             >  1:424484/1‑68 (MQ=255)
                       cGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaag                                         <  1:19517/71‑1 (MQ=255)
                       cGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaaga                                        <  1:1135870/71‑1 (MQ=255)
                        gTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaaga                                        <  1:5491840/71‑1 (MQ=255)
                        gTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaag                                         >  1:4250279/1‑70 (MQ=255)
                           aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTaaa                                           >  1:4204154/1‑65 (MQ=255)
                           aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTaaa                                           >  1:1498318/1‑65 (MQ=255)
                           aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTaaa                                           >  1:4288034/1‑65 (MQ=255)
                           aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaagaaag                                     >  1:3451618/1‑71 (MQ=255)
                           aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaagaaag                                     >  1:1421147/1‑71 (MQ=255)
                           aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTaa                                                       >  1:131692/1‑52 (MQ=255)
                           aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaagaaag                                     >  1:1674959/1‑70 (MQ=255)
                           aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaagaaag                                     >  1:4873859/1‑70 (MQ=255)
                           aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaagaaag                                     >  1:3682130/1‑70 (MQ=255)
                           aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGc                                    <  1:3931975/71‑1 (MQ=255)
                           aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGc                                    <  1:5041616/71‑1 (MQ=255)
                              aCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTaa                                                       <  1:4590777/50‑1 (MQ=255)
                              aCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTaa                                                       <  1:5222628/50‑1 (MQ=255)
                                aTCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTaaaa                                          <  1:2512629/61‑1 (MQ=255)
                                aTCCATATTGCTACGGATCGTCAACAAATAAATTTTGTAATCGGTTAATg                                                     >  1:855902/1‑50 (MQ=255)
                                     tatTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTaa                                            <  1:1691885/54‑1 (MQ=255)
                                     tatTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaaga                                        <  1:4013558/57‑1 (MQ=255)
                                      atTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGtt                                                         <  1:1273873/40‑1 (MQ=38)
                                      atTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAAc                                 >  1:4678301/1‑64 (MQ=255)
                                       ttGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTg                         >  1:2519605/1‑70 (MQ=255)
                                        tGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTaa                                                       >  1:681331/1‑40 (MQ=38)
                                        tGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAAccc                               >  1:2118662/1‑64 (MQ=255)
                                         gCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaaga                                        >  1:1238827/1‑54 (MQ=255)
                                         gCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGa                        >  1:4566857/1‑70 (MQ=255)
                                         gCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGaa                       <  1:1254714/70‑1 (MQ=255)
                                              ggATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaag                                         <  1:1478025/48‑1 (MQ=39)
                                               gATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGa                        >  1:1663202/1‑64 (MQ=255)
                                               gATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGa                        >  1:2509741/1‑64 (MQ=255)
                                                aTCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGc                                    <  1:2110897/51‑1 (MQ=255)
                                                 tCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa              <  1:896133/71‑1 (MQ=255)
                                                 tCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa              <  1:902803/71‑1 (MQ=255)
                                                 tCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa              <  1:4920244/71‑1 (MQ=255)
                                                 tCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa              <  1:32197/71‑1 (MQ=255)
                                                    tCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGaaa            >  1:2483408/1‑71 (MQ=255)
                                                     cAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaaga                                        <  1:400492/42‑1 (MQ=38)
                                                      aGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGaa             <  1:5474625/68‑1 (MQ=255)
                                                      aGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAAt           <  1:5386018/70‑1 (MQ=255)
                                                      aGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCg               >  1:5281952/1‑65 (MQ=255)
                                                        cAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGaaa                      >  1:3044935/1‑56 (MQ=255)
                                                        cAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATgcgc       >  1:2662177/1‑71 (MQ=255)
                                                         aaataaatTTTGTATCGGTTAATGCCGCAGTAaaaga                                        <  1:5092264/37‑1 (MQ=255)
                                                          aataaatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAAccc                               >  1:3222838/1‑46 (MQ=255)
                                                          aataaatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTc                <  1:155803/61‑1 (MQ=255)
                                                          aataaatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGcgcg      <  1:1227966/71‑1 (MQ=255)
                                                           ataaatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATgcgc       <  1:1336422/69‑1 (MQ=255)
                                                           ataaatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATgcgc       <  1:502446/69‑1 (MQ=255)
                                                              aatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGaa                       >  1:220899/1‑50 (MQ=255)
                                                              aatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCt  >  1:4872964/1‑71 (MQ=255)
                                                                ttttGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa              <  1:2690228/57‑1 (MQ=255)
                                                                ttttGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa              <  1:5189778/57‑1 (MQ=255)
                                                                ttttGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATgc         <  1:645381/62‑1 (MQ=255)
                                                                      |                                                              
AGAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCT  >  Exported/975188‑975319

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: