Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,028 C→G 19.2% A829G (GCC→GGC) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,0280CG19.2% 109.0 / 20.4 52A829G (GCC→GGC) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (20/22);  new base G (4/6);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

TTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCT  >  minE/833959‑834094
                                                                     |                                                                  
ttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc                                                                    >  1:583890/1‑70 (MQ=255)
   ttCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATc                                                                >  1:1281837/1‑71 (MQ=255)
       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                                            <  1:921612/70‑3 (MQ=255)
       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                                            <  1:781429/70‑3 (MQ=255)
       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                                            <  1:2048209/70‑3 (MQ=255)
       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                                            <  1:292052/70‑3 (MQ=255)
         cGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                                            >  1:3470978/1‑67 (MQ=255)
         cGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                                             >  1:2037592/1‑68 (MQ=255)
         cGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                                             >  1:1838794/1‑68 (MQ=255)
          gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaag                                                          >  1:1536389/1‑70 (MQ=255)
              tCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                                             <  1:2063502/63‑1 (MQ=255)
              tCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                                             <  1:2575807/63‑1 (MQ=255)
                     cggTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGaa                                              <  1:1211621/71‑1 (MQ=255)
                       gTCGATACGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAg                                             >  1:210192/1‑70 (MQ=255)
                        tCGATCCGGCCTGGGAAGCGAAACGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAAt                                                      <  1:1941501/60‑1 (MQ=255)
                        tCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCg                                                <  1:1616871/66‑1 (MQ=255)
                        tCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgtg                                           >  1:2782389/1‑71 (MQ=255)
                         cGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                                             <  1:2342961/52‑1 (MQ=255)
                           aTCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCcaaa                                                          <  1:3020136/53‑2 (MQ=255)
                            tCCGGCCTTGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCg                                       <  1:945227/71‑1 (MQ=255)
                             ccGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt                                      >  1:1433009/1‑71 (MQ=255)
                             ccGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt                                      >  1:3054163/1‑71 (MQ=255)
                                  cTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATtct                                                   >  1:3820388/1‑53 (MQ=255)
                                   ctggAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGc                                        >  1:2525992/3‑63 (MQ=255)
                                     ggAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATc                              >  1:597351/1‑71 (MQ=255)
                                     ggAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATc                              >  1:1586483/1‑71 (MQ=255)
                                      gAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCt                             <  1:2591893/71‑1 (MQ=255)
                                       aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTg                            <  1:1974148/71‑1 (MQ=255)
                                        aGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGaa                                              >  1:1634938/1‑52 (MQ=255)
                                        aGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTgg                           >  1:2398472/1‑71 (MQ=255)
                                         gCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTggg                          <  1:1754002/71‑1 (MQ=255)
                                           gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTgg                           <  1:2087617/68‑1 (MQ=255)
                                           gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa                        <  1:1517705/71‑1 (MQ=255)
                                           gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa                        <  1:3754192/71‑1 (MQ=255)
                                             aaaGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTaaa                      <  1:3922874/71‑1 (MQ=255)
                                               aGCGACTGGGAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGTGTGCGTGGGTCATCTTGGCaaa                      >  1:2475190/1‑69 (MQ=255)
                                                gCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGa                    <  1:3540844/70‑1 (MQ=255)
                                                   aCTGGGAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGTGTGCGTGGGTCATCTTGGCAAAGAAAcc                <  1:880791/71‑1 (MQ=255)
                                                    cTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa                        <  1:2596003/62‑1 (MQ=255)
                                                     tGGGAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGTGTGCGTGGGTCATCTTGGCAAAGaaa                  >  1:3043177/1‑67 (MQ=255)
                                                       ggAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGAc             >  1:1845737/1‑70 (MQ=255)
                                                       ggAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACa            >  1:378323/1‑71 (MQ=255)
                                                       ggAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACa            >  1:2880556/1‑71 (MQ=255)
                                                       ggAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACa            <  1:716772/71‑1 (MQ=255)
                                                        gAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGTGTGCGTGGGTCATCTTGGCAAAGAAACCGAcg            >  1:3090338/1‑69 (MQ=255)
                                                         aaaTCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATc          <  1:3558644/71‑1 (MQ=255)
                                                          aaTCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCg         <  1:2328367/71‑1 (MQ=255)
                                                             cTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATc                              >  1:335780/1‑47 (MQ=255)
                                                              tACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGaaa                  >  1:2993752/1‑58 (MQ=255)
                                                               aCAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa                        <  1:40990/51‑1 (MQ=255)
                                                               aCAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGaaa                  >  1:1087941/1‑57 (MQ=255)
                                                                  taGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCt                             <  1:354189/42‑1 (MQ=255)
                                                                  aaGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCt  <  1:1274669/70‑1 (MQ=255)
                                                                     |                                                                  
TTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCT  >  minE/833959‑834094

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: