Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,146 T→G 6.1% S202A (TCG→GCG) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,1460TG 6.1% 171.2 / 6.2 66S202A (TCG→GCG) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base G (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCC  >  minE/832076‑832216
                                                                      |                                                                      
cAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCAt                                                                        >  1:1310227/1‑71 (MQ=255)
 aaaaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATc                                                                       >  1:694688/1‑71 (MQ=255)
 aaaaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATc                                                                       >  1:1724244/1‑71 (MQ=255)
  aaaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCg                                                                      <  1:2911752/71‑1 (MQ=255)
      cTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTg                                                                  <  1:1575124/71‑1 (MQ=255)
           gCCCGTACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGc                                                             >  1:2774026/1‑71 (MQ=255)
            cccGGACCGTGTTGCTGGTTTCCCCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGt                                                                   >  1:1361467/1‑63 (MQ=255)
              cGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAt                                                          >  1:221494/1‑71 (MQ=255)
                 aCCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtctc                                                       <  1:537180/71‑1 (MQ=255)
                 aCCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtct                                                        >  1:1812175/1‑70 (MQ=255)
                  ccGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtctct                                                      >  1:2979674/1‑71 (MQ=255)
                     tgtTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGc                                                   >  1:3946193/1‑71 (MQ=255)
                       ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                  <  1:1586102/70‑1 (MQ=255)
                       ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                  >  1:459835/1‑70 (MQ=255)
                       ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                  >  1:3406514/1‑70 (MQ=255)
                       ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                  >  1:3746947/1‑70 (MQ=255)
                           tGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                  <  1:2733639/66‑1 (MQ=255)
                             gTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCt                                                            <  1:1784274/54‑1 (MQ=255)
                              tttCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACTCTATCTCTCGCTGATTGGCg                                          <  1:4012040/71‑1 (MQ=255)
                                tctcGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCgg                                         >  1:1282069/1‑70 (MQ=255)
                                    gCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTgct                                               >  1:1212902/1‑58 (MQ=255)
                                    gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                >  1:570453/1‑59 (MQ=255)
                                    gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                >  1:590717/1‑59 (MQ=255)
                                    gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTa                                       <  1:3483918/68‑1 (MQ=255)
                                    gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACt                                     <  1:471673/70‑1 (MQ=255)
                                        aTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt                                <  1:3176962/71‑1 (MQ=255)
                                         ttCCGGCAATGTCGATGGTTTCTTACGCATCGGGTg                                                                  <  1:4108291/36‑1 (MQ=255)
                                         ttCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                  <  1:1805796/52‑1 (MQ=255)
                                          tCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCgg                                         >  1:3544688/1‑60 (MQ=255)
                                           ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa                              <  1:4095222/70‑1 (MQ=255)
                                           ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa                              <  1:407164/70‑1 (MQ=255)
                                           ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa                              <  1:1000774/70‑1 (MQ=255)
                                           ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa                              <  1:1498431/70‑1 (MQ=255)
                                             ggCAATGTCGATGGTTTCTTACGCATCGGGTTCACGCt                                                            <  1:3525044/38‑1 (MQ=255)
                                               cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTa                                       >  1:3461843/1‑57 (MQ=255)
                                               cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTc                         <  1:3395428/71‑1 (MQ=255)
                                                aaTGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCt                        <  1:1805607/71‑1 (MQ=255)
                                                 aTGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                  <  1:3450453/44‑1 (MQ=255)
                                                  tGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTAc                      <  1:2054079/71‑1 (MQ=255)
                                                  tGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTAc                      <  1:1468191/71‑1 (MQ=255)
                                                  tGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTAc                      <  1:1677780/71‑1 (MQ=255)
                                                   gTCGATGGTTTCTTACGCATCGGGTGCACGCTAtct                                                        <  1:1254508/36‑1 (MQ=255)
                                                   gTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCgg                                         >  1:269240/1‑51 (MQ=255)
                                                    tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa                              <  1:1209788/61‑1 (MQ=255)
                                                    tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGa                    >  1:1810687/1‑71 (MQ=255)
                                                    tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGa                    >  1:2714946/1‑71 (MQ=255)
                                                    tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACATGCTTAAGCTTCTACGa                    >  1:2844402/1‑71 (MQ=255)
                                                     cGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTTCTTGGCGGCACCTGTTTAAGTTTc                         <  1:1857182/65‑1 (MQ=255)
                                                     cGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGAc                   >  1:1667866/1‑71 (MQ=255)
                                                        tGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTgg                >  1:1195525/1‑71 (MQ=255)
                                                         ggTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTc                         <  1:1603052/61‑1 (MQ=255)
                                                         ggTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCgg                                         >  1:540171/1‑45 (MQ=255)
                                                          gTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTg                                                 <  1:3454078/36‑1 (MQ=255)
                                                          gTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTg                                                 <  1:2334187/36‑1 (MQ=255)
                                                          gTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTa              >  1:3624862/1‑71 (MQ=255)
                                                           tttCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTAc             >  1:3111587/1‑71 (MQ=255)
                                                            ttCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                >  1:108775/1‑35 (MQ=255)
                                                            ttCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCgg                                         >  1:1441310/1‑42 (MQ=255)
                                                             tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACt            >  1:2482458/1‑70 (MQ=255)
                                                             tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTg           >  1:2111103/1‑71 (MQ=255)
                                                              cTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTg                                   >  1:1470241/1‑46 (MQ=255)
                                                              cTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTg           <  1:93326/70‑1 (MQ=255)
                                                                   gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTAc                      <  1:2627211/54‑1 (MQ=255)
                                                                   gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt     <  1:3078385/71‑1 (MQ=255)
                                                                   gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt     <  1:185304/71‑1 (MQ=255)
                                                                     aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGc   <  1:526869/71‑1 (MQ=255)
                                                                     aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGc   <  1:12984/71‑1 (MQ=255)
                                                                      tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTgg                <  1:3933207/57‑1 (MQ=255)
                                                                      tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc  <  1:1493655/71‑1 (MQ=255)
                                                                      |                                                                      
CAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCC  >  minE/832076‑832216

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: