Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,184,926 A→G 11.1% E212E (GAA→GAG gltD → glutamate synthase, 4Fe‑4S protein, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,184,9260AG11.1% 119.8 / 7.1 45E212E (GAA→GAGgltDglutamate synthase, 4Fe‑4S protein, small subunit
Reads supporting (aligned to +/- strand):  ref base A (22/18);  new base G (2/3);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

CCGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACGATCTGCTGAGTGATTACGATG  >  Exported/2184858‑2184990
                                                                    |                                                                
ccGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAAtt                                                                >  1:594712/1‑71 (MQ=255)
ccGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAAt                                                                 <  1:356897/70‑1 (MQ=255)
     cTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAAc                                                           <  1:2523366/71‑1 (MQ=255)
     cTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAAc                                                           >  1:2475317/1‑71 (MQ=255)
     cTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAAc                                                           >  1:713448/1‑71 (MQ=255)
       tCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACt                                                          >  1:1891348/1‑70 (MQ=255)
        cAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTc                                                         >  1:2710120/1‑70 (MQ=255)
        cAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTc                                                         <  1:1614486/70‑1 (MQ=255)
        cAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTc                                                         >  1:899802/1‑70 (MQ=255)
         aaGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCaa                                                       <  1:2869729/71‑1 (MQ=255)
         aaGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCaa                                                       <  1:655852/71‑1 (MQ=255)
              ggAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCAATAcc                                                   >  1:152157/1‑70 (MQ=255)
                aaaaaGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGaa                                                <  1:3436043/71‑1 (MQ=255)
                 aaaaGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGCTATTGAATTCAAACTCAATACCGAAg                                               >  1:91998/1‑71 (MQ=255)
                       ggTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAGTTTGAATTCAAt                                                      <  1:642513/58‑1 (MQ=37)
                       ggTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAGTTTGAATTCAAt                                                      <  1:1393422/58‑1 (MQ=37)
                       ggTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAGTTTGAATTCAAt                                                      <  1:2971377/58‑1 (MQ=37)
                         tAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAGTTTGAATTCAATAc                                                    >  1:1667051/1‑58 (MQ=25)
                         tAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCcgc                                       <  1:3543364/71‑1 (MQ=255)
                           aTGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGa                                     <  1:1393146/71‑1 (MQ=255)
                            tGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGAc                                    <  1:2744980/71‑1 (MQ=255)
                            tGACGCGTCGCCGTGAAATCTTCACCGGCATGGTATTGAATTCAAACTCTATACCGAAGTGGGCCGCGACg                                   >  1:3632903/1‑71 (MQ=255)
                              aCGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGt                                  <  1:3713283/71‑1 (MQ=255)
                                gcgTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTAc                                <  1:2208508/71‑1 (MQ=255)
                                 cgTCGCCGTTAAATCTTCACCGGCATGGGTATTGAATTCaa                                                             <  1:367648/41‑1 (MQ=255)
                                    cGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCt                            >  1:1743736/1‑71 (MQ=255)
                                        gTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGAc                        >  1:3074059/1‑71 (MQ=255)
                                         tGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCaa                                                       <  1:3134649/39‑1 (MQ=255)
                                         tGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACg                       <  1:2016355/71‑1 (MQ=255)
                                           aaaTCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAg                              >  1:2574364/1‑62 (MQ=255)
                                           aaaTCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGc                             >  1:3428586/1‑63 (MQ=255)
                                               cTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACGATctgc                 >  1:2816672/1‑71 (MQ=255)
                                               cTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACGATctgc                 >  1:2217924/1‑71 (MQ=255)
                                                  cACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACGATctgc                 <  1:577511/68‑1 (MQ=255)
                                                   aCCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACGATCTGCTGAg             <  1:3436563/71‑1 (MQ=255)
                                                     cGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTgg                          <  1:1047854/56‑1 (MQ=255)
                                                             gTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTa                                 >  1:1239885/1‑41 (MQ=255)
                                                             gTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTa                                 >  1:2737139/1‑41 (MQ=255)
                                                             gTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTa                                 >  1:1870738/1‑41 (MQ=255)
                                                              tATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACGATCTGCTGAGTGATTACGATg  >  1:2580327/1‑71 (MQ=255)
                                                              tATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACGATCTGCTGAGTGATTACGATg  >  1:2435378/1‑71 (MQ=255)
                                                               tttGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACGa                      <  1:1691901/49‑1 (MQ=255)
                                                               aTTGAGTTTGAATTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACGATCTGCTGAg             >  1:919780/1‑59 (MQ=37)
                                                                  gAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACGATctgc                 >  1:1953958/1‑52 (MQ=255)
                                                                  gAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACGATCTGCTGAGTGATTACGa    >  1:3449711/1‑65 (MQ=255)
                                                                    aTTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTg                           <  1:3495753/40‑1 (MQ=255)
                                                                    aTTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTg                           <  1:3680535/40‑1 (MQ=255)
                                                                    aTTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTg                           <  1:2645989/40‑1 (MQ=255)
                                                                    aTTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTg                           <  1:971533/40‑1 (MQ=255)
                                                                    |                                                                
CCGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAATTCAAACTCAATACCGAAGTGGGCCGCGACGTACAGCTGGACGATCTGCTGAGTGATTACGATG  >  Exported/2184858‑2184990

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: