Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,270 G→A 45.7% A576A (GCG→GCA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,2700GA45.7% 20.1 / 81.4 70A576A (GCG→GCAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (23/15);  new base A (18/14);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGG  >  minE/833203‑833339
                                                                   |                                                                     
aCTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCAc                                                                    >  1:3275974/1‑71 (MQ=255)
   aTCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTAt                                                                 >  1:1231917/1‑71 (MQ=255)
    tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTa                                                                  <  1:2326530/69‑1 (MQ=255)
    tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                                                                >  1:2480508/1‑71 (MQ=255)
    tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                                                                >  1:407809/1‑71 (MQ=255)
    tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGAGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                                                                <  1:472933/71‑1 (MQ=255)
         gTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                                                           <  1:4566378/71‑1 (MQ=255)
         gTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                                                           <  1:1916905/71‑1 (MQ=255)
          tCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                                                           >  1:806225/1‑70 (MQ=255)
              aTCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAg                                                        >  1:2821760/1‑69 (MQ=255)
              aTCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGa                                                      <  1:2944028/71‑1 (MQ=255)
               tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                                                  >  1:2691954/1‑58 (MQ=255)
               tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                                                  >  1:49515/1‑58 (MQ=255)
               tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                                                  >  1:3851377/1‑58 (MQ=255)
               tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                                                  >  1:5152176/1‑58 (MQ=255)
               tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                                                  >  1:3857956/1‑58 (MQ=255)
                 aataTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaa                                                   <  1:3016967/68‑1 (MQ=255)
                 aataTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaa                                                   <  1:4184778/68‑1 (MQ=255)
                 aataTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaa                                                    >  1:2981445/4‑70 (MQ=255)
                 aaCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa                                                    <  1:3493599/70‑1 (MQ=255)
                 aaCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa                                                    <  1:2167987/70‑1 (MQ=255)
                    aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                                                             >  1:4690972/1‑58 (MQ=255)
                    aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                                                             >  1:188112/1‑58 (MQ=255)
                    aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                                                             >  1:1262426/1‑58 (MQ=255)
                    aTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa                                                  >  1:4464413/1‑69 (MQ=255)
                    aTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa                                                  >  1:4304371/1‑69 (MQ=255)
                    aTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa                                                  >  1:3406557/1‑69 (MQ=255)
                     tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACt                                                                   <  1:5102118/51‑1 (MQ=255)
                     tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa                                                  >  1:4945892/1‑68 (MQ=255)
                     tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTg                                               >  1:3034003/1‑71 (MQ=255)
                          atcttCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa                                                  >  1:3438767/4‑63 (MQ=255)
                            tttCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCg                                         >  1:1400618/1‑70 (MQ=255)
                            cttCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGccca                                         >  1:949659/2‑69 (MQ=255)
                              tCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCa                                       <  1:5222416/70‑1 (MQ=255)
                                tGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCa                                                                     >  1:1644431/1‑38 (MQ=255)
                                tGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                                                 >  1:4963301/1‑58 (MQ=255)
                                tGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                                                 >  1:2936667/1‑58 (MQ=255)
                                 gCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa                                                    >  1:3734738/1‑54 (MQ=255)
                                   ggCTGTTTCGGTCAGAGCGGGGGCGGCTGGGCGCTCTATGTa                                                              >  1:2750740/1‑42 (MQ=38)
                                   ggCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                                                           >  1:3400123/1‑45 (MQ=255)
                                   ggCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                                                           >  1:3760021/1‑45 (MQ=255)
                                       gtgtCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTgg                             >  1:5529711/1‑71 (MQ=255)
                                        tgtCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTgg                             <  1:1003885/70‑1 (MQ=255)
                                        tgtCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTggt                            <  1:4121250/71‑2 (MQ=255)
                                        tgtCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTgg                             <  1:1119725/70‑1 (MQ=255)
                                          tCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                                                 >  1:1166919/1‑48 (MQ=255)
                                          tCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCaa                                      <  1:2999885/59‑1 (MQ=255)
                                          tCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCaa                                      <  1:1294849/59‑1 (MQ=255)
                                           cGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGc                         >  1:5440272/1‑71 (MQ=255)
                                            ggTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                                                 >  1:3649610/1‑46 (MQ=255)
                                            ggTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                                                 >  1:3089178/1‑46 (MQ=255)
                                              tCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGcccc                      <  1:3377181/71‑3 (MQ=255)
                                                aaagCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTg                    <  1:616480/69‑1 (MQ=255)
                                                  agCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGc                  <  1:4208992/71‑1 (MQ=255)
                                                   gCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGcc                 <  1:2457101/71‑2 (MQ=255)
                                                    cGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGcct                <  1:2909029/71‑3 (MQ=255)
                                                       tggcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTg             <  1:244647/70‑1 (MQ=255)
                                                         gcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGcg           >  1:4792371/1‑70 (MQ=255)
                                                         gcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGcg           <  1:2734165/71‑2 (MQ=255)
                                                         gcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGcg           <  1:590398/71‑2 (MQ=255)
                                                         gcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGcg           <  1:1126363/71‑2 (MQ=255)
                                                         gcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGc            >  1:96120/1‑70 (MQ=255)
                                                           ggcTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCt         >  1:5004806/1‑71 (MQ=255)
                                                            gcTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCtt        <  1:3890155/71‑1 (MQ=255)
                                                             cTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGCGCtcg       <  1:1484523/71‑3 (MQ=255)
                                                             cTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGCGCtc        >  1:385212/1‑69 (MQ=255)
                                                                ggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTgg                             >  1:907848/1‑46 (MQ=255)
                                                                ggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGAc     <  1:1988872/70‑1 (MQ=255)
                                                                ggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGAc     <  1:21541/70‑1 (MQ=255)
                                                                ggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACt    >  1:2654840/1‑71 (MQ=255)
                                                                   gcACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTgg  >  1:3041872/1‑70 (MQ=255)
                                                                   |                                                                     
ACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGG  >  minE/833203‑833339

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: