Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,132 T→C 13.7% A530A (GCT→GCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,1320TC13.7% 176.8 / 17.9 73A530A (GCT→GCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base C (8/2);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GGCCGAGCAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGA  >  minE/833063‑833197
                                                                     |                                                                 
ggCCGAGCAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCt                                                                   >  1:496668/1‑70 (MQ=255)
    gAGCAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATaa                                                              <  1:998684/71‑1 (MQ=255)
    gAGCAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATaa                                                              <  1:812631/71‑1 (MQ=255)
     aGCAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATaaa                                                             <  1:392896/71‑1 (MQ=255)
       cAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAAc                                                           >  1:4503973/1‑71 (MQ=255)
              ccGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACgg                                                     >  1:2061527/1‑70 (MQ=255)
              ccGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGt                                                    >  1:989685/1‑71 (MQ=255)
                ggCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCg                                                  >  1:1446035/1‑71 (MQ=255)
                    tttCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:5458287/71‑1 (MQ=255)
                        tcGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCgatga                                          >  1:4160535/1‑71 (MQ=255)
                        tcGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCgatga                                          >  1:2032117/1‑71 (MQ=255)
                         cGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCgatgat                                         >  1:4495796/1‑71 (MQ=255)
                          gCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCgatgat                                         >  1:3294832/1‑70 (MQ=255)
                          gCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGAtt                                        <  1:1660017/71‑1 (MQ=255)
                          gCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGAtt                                        <  1:4631708/71‑1 (MQ=255)
                                aGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTAtcgtc                                  >  1:4424897/1‑71 (MQ=255)
                                  attattCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTAtcgtcg                                 >  1:3876212/1‑70 (MQ=255)
                                  attattCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCgg                                >  1:5486270/1‑71 (MQ=255)
                                    tattCGCATCGCCCGTGAATTTGCCGATAACGCTGa                                                                 >  1:129051/1‑36 (MQ=255)
                                    tattCGCATCGCCCGTGAATTTGCCGATAACGCTGa                                                                 >  1:3794782/1‑36 (MQ=255)
                                    tattCGCATCGCCCGTGAATTTGCCGATAACGCTGa                                                                 >  1:2120276/1‑36 (MQ=255)
                                    tattCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGt                                                    >  1:3776606/1‑49 (MQ=255)
                                     attCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGc                             >  1:2310739/1‑71 (MQ=255)
                                      ttCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGc                             <  1:4779345/70‑1 (MQ=255)
                                      ttCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCg                            >  1:448964/1‑71 (MQ=255)
                                      ttCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCg                            >  1:4875783/1‑71 (MQ=255)
                                      ttCGCATCGCACGTGAATTTGCCGATAACGCTGATaa                                                              <  1:388856/37‑1 (MQ=255)
                                         ccATCGCCCGTGAATTTGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTATCCTCGGCGCagg                          >  1:3470908/2‑70 (MQ=255)
                                            tCGCCCGTGAATTTGCCGATACCGCCCATAAAACgt                                                         <  1:5321883/36‑2 (MQ=21)
                                            tCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCt                       >  1:5553360/1‑70 (MQ=255)
                                              gCCCGTGAATTTGCCGATAGCGCTGATAAAACGCACGGTCGTTCGATGATTAtcgtc                                  >  1:4249936/1‑57 (MQ=255)
                                              gCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTAtcgtc                                  >  1:1539867/1‑57 (MQ=255)
                                              gCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGa                     <  1:4305386/70‑1 (MQ=255)
                                               cccGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAAc                   <  1:1926014/71‑1 (MQ=255)
                                                  gTGAATTTGCCGATACCGCCCATAAAACGCATGGGCGCTCgatg                                           >  1:4957753/1‑44 (MQ=16)
                                                  gTGAATTTGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTATCCTCGGCGCAGGTGTTAACCAc                <  1:3346124/71‑1 (MQ=255)
                                                   tGAATTTGCCGATAACGCTGATAAAACGCACGGTCg                                                  <  1:3733505/36‑1 (MQ=255)
                                                   tGAATTTGCCGATAACGCTGATAAAACGCACGGTCg                                                  <  1:776077/36‑1 (MQ=255)
                                                   tGAATTTGCCGATAACGCTGATAAAACGCACGGTCg                                                  <  1:2804959/36‑1 (MQ=255)
                                                   tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:4147449/40‑1 (MQ=255)
                                                   tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:5602389/40‑1 (MQ=255)
                                                   tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:4947925/40‑1 (MQ=255)
                                                   tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:1478583/40‑1 (MQ=255)
                                                   tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:5063282/40‑1 (MQ=255)
                                                   tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:4130145/40‑1 (MQ=255)
                                                   tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:4379060/40‑1 (MQ=255)
                                                   tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:2529688/40‑1 (MQ=255)
                                                    gAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTg                      <  1:2100598/63‑1 (MQ=255)
                                                    gAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTg                      <  1:725909/63‑1 (MQ=255)
                                                    gAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTg              >  1:3616516/1‑71 (MQ=255)
                                                     aaTTTGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTAtcctc                                  >  1:1112654/1‑50 (MQ=255)
                                                     aaTTTGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTAtcctc                                  >  1:2958871/1‑50 (MQ=255)
                                                     aaTTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTg                      <  1:3599317/62‑1 (MQ=255)
                                                     aaTTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTgg             >  1:3454848/1‑71 (MQ=255)
                                                      aTTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTAtcgtc                                  >  1:963441/1‑49 (MQ=255)
                                                      aTTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGt            <  1:5615322/71‑1 (MQ=255)
                                                        ttGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCgg                                >  1:2504301/1‑49 (MQ=255)
                                                         tggcGATACCGCCCATAAAACGCATGGGCGCTCGATGATTATCCTCGGCGCAGGTGTTAACCACTGGTAt          >  1:5308718/4‑70 (MQ=255)
                                                         tGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTAtc                                     >  1:1594714/1‑43 (MQ=16)
                                                         tGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTATCCTCGGCGCaggtg                        >  1:964578/1‑54 (MQ=255)
                                                         tGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTATCCTCGGCGCAGGTGTTAACCa                 >  1:3023558/1‑63 (MQ=255)
                                                         tGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTAt          >  1:2204781/1‑70 (MQ=255)
                                                              aTAACGCTGATAAAACGCACGGTCGTTCGATGATTAt                                      >  1:5328840/1‑37 (MQ=255)
                                                              aTAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    >  1:2843909/1‑71 (MQ=255)
                                                              aTAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    >  1:511904/1‑71 (MQ=255)
                                                              aTAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    >  1:5594759/1‑71 (MQ=255)
                                                              aTAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    >  1:3836548/1‑71 (MQ=255)
                                                                aaCGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    <  1:3811269/69‑1 (MQ=255)
                                                                aaCGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    <  1:4694658/69‑1 (MQ=255)
                                                                aaCGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    <  1:3817943/69‑1 (MQ=255)
                                                                aaCGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGa  >  1:3457660/1‑71 (MQ=255)
                                                                 aCGCTGATAAAACGCACGGTCGTTCGATGATTATCGTAGGTGGGGGGCTGAACCACTGGTATCACCTc    <  1:5588011/68‑1 (MQ=255)
                                                                  cGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCg   >  1:3274504/1‑68 (MQ=255)
                                                                   gCTGATAAAACGCACGGTCGTTCGATGATTAtcgt                                   <  1:2252131/35‑1 (MQ=255)
                                                                     |                                                                 
GGCCGAGCAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGA  >  minE/833063‑833197

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: