Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,223 G→C 21.3% A227A (GCG→GCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,2230GC21.3% 148.3 / 34.7 75A227A (GCG→GCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (33/26);  new base C (6/10);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACAT  >  minE/832155‑832291
                                                                    |                                                                    
cGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtc                                                                    <  1:168697/71‑1 (MQ=255)
  cTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtctc                                                                  <  1:1356289/71‑1 (MQ=255)
     tctgtcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCg                                                                <  1:3323946/70‑1 (MQ=255)
     tctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                 >  1:5193665/1‑69 (MQ=255)
      ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                                               >  1:5165428/1‑70 (MQ=255)
        ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                 >  1:1235824/1‑66 (MQ=255)
        ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                 >  1:1526610/1‑66 (MQ=255)
        ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaaa                                                            >  1:5052580/1‑71 (MQ=255)
         tcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCgatg                                                             <  1:761100/69‑4 (MQ=255)
         tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAc                                                           >  1:3511473/1‑71 (MQ=255)
           gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                                                          <  1:569964/70‑1 (MQ=255)
           gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                                                          <  1:2597733/70‑1 (MQ=255)
              gATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                                                          <  1:3422109/67‑1 (MQ=255)
                ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTggg                                                      >  1:5620596/1‑69 (MQ=255)
                ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGc                                                    <  1:1200017/71‑1 (MQ=255)
                  ggCGGTACATGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCg                                                   <  1:1009608/70‑1 (MQ=255)
                      gTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCaa                                              <  1:3585810/71‑1 (MQ=255)
                        aCTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                                               <  1:663507/68‑1 (MQ=255)
                         cTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAAc                                            <  1:4644879/70‑1 (MQ=255)
                         cTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACt                                           <  1:163062/71‑1 (MQ=255)
                          ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                                             <  1:5615593/52‑1 (MQ=255)
                           tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                 >  1:1723349/1‑47 (MQ=255)
                           tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                 >  1:1707352/1‑47 (MQ=255)
                           tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                 >  1:4719233/1‑47 (MQ=255)
                             cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACg                                       <  1:1722190/71‑1 (MQ=255)
                               tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGt                                      <  1:3382378/70‑1 (MQ=255)
                               tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGt                                      <  1:3375919/70‑1 (MQ=255)
                               tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                                                  >  1:5256285/1‑58 (MQ=255)
                               tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                                                  >  1:2239762/1‑58 (MQ=255)
                                aaGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTa                                     <  1:1006427/70‑1 (MQ=255)
                                aaGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                                               >  1:3515463/1‑44 (MQ=255)
                                aaGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGt                                      <  1:4166571/69‑1 (MQ=255)
                                   tttCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGc                                                    <  1:3806603/51‑1 (MQ=255)
                                     tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGaa                               >  1:1655704/1‑71 (MQ=255)
                                     tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGa                                <  1:3303868/70‑1 (MQ=255)
                                      cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                <  1:999977/69‑1 (MQ=255)
                                      cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                <  1:2123804/69‑1 (MQ=255)
                                      cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                              >  1:3042141/1‑71 (MQ=255)
                                       tACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                              >  1:1807426/1‑70 (MQ=255)
                                       tACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATc                             >  1:1283628/1‑71 (MQ=255)
                                       tACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATc                             >  1:1396551/1‑71 (MQ=255)
                                        aCGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                              >  1:4469270/1‑69 (MQ=255)
                                         cGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATc                             <  1:5684819/69‑1 (MQ=255)
                                         cGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                              <  1:90272/68‑1 (MQ=255)
                                         cGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATct                            >  1:4196601/1‑70 (MQ=255)
                                            cTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATctgctg                        <  1:1065204/71‑1 (MQ=255)
                                             tGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCgg                                 <  1:1009180/61‑1 (MQ=255)
                                              ggTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGAc                      >  1:5140846/1‑71 (MQ=255)
                                                  ttGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGaa                               <  1:2171945/57‑1 (MQ=255)
                                                  ttGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGAc                      >  1:1859496/2‑67 (MQ=255)
                                                  cTGCGTCTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACt                     <  1:1194552/68‑1 (MQ=255)
                                                  cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                                                  >  1:3899584/1‑39 (MQ=255)
                                                  cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                                               <  1:3609614/42‑1 (MQ=255)
                                                  cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                                               <  1:2842309/42‑1 (MQ=255)
                                                    gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa                 >  1:180735/1‑70 (MQ=255)
                                                    gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa                 >  1:3097322/1‑70 (MQ=255)
                                                    gCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGa                                         >  1:3131241/1‑46 (MQ=255)
                                                    gCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAAt                              >  1:3352922/1‑57 (MQ=255)
                                                    gCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAAt                              >  1:5092699/1‑57 (MQ=255)
                                                    gCGACATGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGACTGGTa                 >  1:2420630/1‑70 (MQ=255)
                                                     cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTg                    >  1:2286052/1‑66 (MQ=255)
                                                     cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                <  1:3566441/70‑1 (MQ=255)
                                                     cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                <  1:1238824/70‑1 (MQ=255)
                                                     cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACa               <  1:4997379/71‑1 (MQ=255)
                                                     cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACa               <  1:945171/71‑1 (MQ=255)
                                                      gACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACaa              <  1:3450789/71‑1 (MQ=255)
                                                       aCTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAAc             >  1:3771725/1‑71 (MQ=255)
                                                         ttGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTc           >  1:418955/1‑71 (MQ=255)
                                                            cctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAg        >  1:2857378/1‑71 (MQ=255)
                                                             ctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTa                                     >  1:4872348/1‑41 (MQ=255)
                                                             ctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTa                                     >  1:4903174/1‑41 (MQ=255)
                                                                ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTa     >  1:5526511/1‑70 (MQ=255)
                                                                ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTa     >  1:3282601/1‑70 (MQ=255)
                                                                ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGATa     >  1:2386521/1‑70 (MQ=255)
                                                                   cGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcat  >  1:5357351/1‑70 (MQ=255)
                                                                    |                                                                    
CGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACAT  >  minE/832155‑832291

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: