Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,145 A→C 6.2% A201A (GCA→GCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,1450AC 6.2% 163.3 / 6.5 65A201A (GCA→GCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base C (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGC  >  minE/832076‑832215
                                                                     |                                                                      
cAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCAt                                                                       >  1:1310227/1‑71 (MQ=255)
 aaaaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATc                                                                      >  1:1724244/1‑71 (MQ=255)
 aaaaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATc                                                                      >  1:694688/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:3746947/1‑70 (MQ=255)
                       ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                 <  1:1586102/70‑1 (MQ=255)
                       ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                 >  1:3406514/1‑70 (MQ=255)
                       ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                 >  1:459835/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:590717/1‑59 (MQ=255)
                                    gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                               >  1:570453/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:407164/70‑1 (MQ=255)
                                           ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa                             <  1:4095222/70‑1 (MQ=255)
                                           ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa                             <  1:1498431/70‑1 (MQ=255)
                                           ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa                             <  1:1000774/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)
                                                                     |                                                                      
CAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGC  >  minE/832076‑832215

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: