Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 703,826 T→C 5.9% L273L (TTA→TTG ymcA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported703,8260TC 5.9% 194.8 / 2.2 67L273L (TTA→TTGymcAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (43/20);  new base C (2/2);  total (45/22)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGCGCAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGC  >  Exported/703758‑703896
                                                                    |                                                                      
ggCGCAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTaa                                                                      >  1:4913628/1‑71 (MQ=255)
    cAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTg                                                                  <  1:5627962/71‑1 (MQ=255)
    cAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTg                                                                  <  1:2158183/71‑1 (MQ=255)
     aGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTg                                                                  >  1:1840702/1‑70 (MQ=255)
           atatGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCa                                                            <  1:877951/70‑1 (MQ=255)
           atatGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCa                                                            <  1:4355681/70‑1 (MQ=255)
             atGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCaaa                                                          >  1:5579130/1‑70 (MQ=255)
             atGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCaaa                                                          >  1:2177605/1‑70 (MQ=255)
                ccGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTg                                                                  <  1:4099818/59‑1 (MQ=255)
                ccGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTg                                                                  <  1:394051/59‑1 (MQ=255)
                 cGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATc                                                      >  1:999812/1‑70 (MQ=255)
                 cGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATc                                                      >  1:1110835/1‑70 (MQ=255)
                 cGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCa                                                     <  1:1662944/71‑1 (MQ=255)
                  ggTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCaa                                                    >  1:109960/1‑71 (MQ=255)
                  ggTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCaa                                                    >  1:3381006/1‑71 (MQ=255)
                     tttGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGa                                                                 >  1:2501632/1‑55 (MQ=255)
                     tttGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGa                                                                 >  1:4840384/1‑55 (MQ=255)
                         gCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACg                                                  <  1:5654407/66‑1 (MQ=255)
                          cGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTAc                                            <  1:2031798/71‑1 (MQ=255)
                            tATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGc                                                              <  1:2895887/51‑1 (MQ=255)
                            tATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACgt                                                 <  1:2692703/64‑1 (MQ=255)
                            tATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACgt                                                 <  1:2771650/64‑1 (MQ=255)
                              tCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTaa                                                                      >  1:5331475/1‑41 (MQ=255)
                              tCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACcgcg                                        >  1:5338065/1‑71 (MQ=255)
                                gcgcAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGtt                                      <  1:4107375/71‑1 (MQ=255)
                                gcgcAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGtt                                      <  1:2274816/71‑1 (MQ=255)
                                   cagggcagggCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                                   <  1:4726849/71‑1 (MQ=255)
                                    agggcagggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                                   <  1:1171639/70‑1 (MQ=255)
                                     gggcagggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATaa                                 >  1:1433550/1‑71 (MQ=255)
                                     gggcagggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATaa                                 >  1:1789026/1‑71 (MQ=255)
                                       gacgggcGCAGATCGTTAAAGTTGGTGCGCAACGTAACGCCAAACATAAAGGTGTTGCCACGTTCATAGCt                               <  1:1149429/68‑1 (MQ=255)
                                        cagggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACt                               >  1:5634674/1‑70 (MQ=255)
                                        acgggcGCAGATCGTTAAAGTTGGTGCGCAACGTAACGCCAAACATAAAGgtg                                                >  1:3028354/3‑53 (MQ=255)
                                         agggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACt                               >  1:2925019/1‑69 (MQ=255)
                                          gggcGCAGATCGTTAAAGTTGGTGCGCAACGTAACGCCAAACATaaagg                                                  <  1:3648392/49‑3 (MQ=255)
                                            gcgcAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTAcc                                           >  1:3580591/1‑54 (MQ=255)
                                            gcgcAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTAcc                                           >  1:2101668/1‑54 (MQ=255)
                                            gcgcAGATCGTTAAAGTTGGTGCGCAACGTAACGCCAAACATAAAGGTGTTGCCACGTTCATAGCTAAGGt                          >  1:5142751/1‑71 (MQ=255)
                                              gcAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTg                        <  1:4706880/71‑1 (MQ=255)
                                               cAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACgtg                                                <  1:3158008/46‑1 (MQ=255)
                                                aGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTgag                      <  1:3768694/71‑1 (MQ=255)
                                                 gATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGtt                                              <  1:2110083/46‑1 (MQ=255)
                                                 gATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGtt                         >  1:2571267/1‑67 (MQ=255)
                                                 gATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGtt                         >  1:1829395/1‑67 (MQ=255)
                                                  aTCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGtt                         >  1:185409/1‑66 (MQ=255)
                                                  aTCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACAGCGTTCATAACTCAGGtt                         >  1:5658889/1‑66 (MQ=255)
                                                    cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCaa                                                    >  1:628646/1‑37 (MQ=255)
                                                    cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCaa                                                    >  1:1864379/1‑37 (MQ=255)
                                                    cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAAc                                                   >  1:921970/1‑38 (MQ=255)
                                                    cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAAc                                                   >  1:3437851/1‑38 (MQ=255)
                                                    cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgg                           >  1:80547/1‑62 (MQ=255)
                                                    cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgg                           >  1:1087217/1‑62 (MQ=255)
                                                    cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgg                           >  1:2300236/1‑62 (MQ=255)
                                                    cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgc                           >  1:76829/1‑61 (MQ=255)
                                                    cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCt                  >  1:4514070/1‑71 (MQ=255)
                                                     gTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                       >  1:3238848/1‑65 (MQ=255)
                                                     gTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                       >  1:4265829/1‑65 (MQ=255)
                                                       tGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGcc               >  1:3591559/1‑71 (MQ=255)
                                                       tGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGcc               >  1:3713938/1‑71 (MQ=255)
                                                        gAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                       >  1:5049547/1‑62 (MQ=255)
                                                        gAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                       >  1:2545410/1‑62 (MQ=255)
                                                        gAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                       >  1:644927/1‑62 (MQ=255)
                                                             tGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                                   >  1:2886105/1‑45 (MQ=255)
                                                             tGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                                   >  1:3567367/1‑45 (MQ=255)
                                                             tGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTg         >  1:3031926/1‑71 (MQ=255)
                                                                  cGTAACGTGAAGCCAAACATCAACTTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAgc    >  1:2481980/1‑71 (MQ=255)
                                                                  cGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAgc    >  1:1838828/1‑71 (MQ=255)
                                                                    tAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACt                               <  1:5390072/42‑1 (MQ=255)
                                                                    tAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAgcgc  <  1:3595754/71‑1 (MQ=255)
                                                                    tAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAgcg   >  1:4762669/1‑70 (MQ=255)
                                                                    |                                                                      
GGCGCAGGTTGATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGC  >  Exported/703758‑703896

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: