Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 832,212 T→C 20.2% L224L (TTG→CTG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported832,2120TC20.2% 160.4 / 34.0 82L224L (TTG→CTG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base C (6/10);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

CGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTC  >  Exported/832142‑832282
                                                                      |                                                                      
cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACt                                                                        <  1:4281585/71‑1 (MQ=255)
 gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                                       <  1:3784066/71‑1 (MQ=255)
   aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGc                                                                     >  1:3964838/1‑71 (MQ=255)
    tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc                                                                    <  1:4910086/71‑1 (MQ=255)
    tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc                                                                    <  1:2835482/71‑1 (MQ=255)
    tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc                                                                    <  1:121244/71‑1 (MQ=255)
      gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc                                                                  >  1:4599508/1‑71 (MQ=255)
      gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcct                                                                   <  1:4498985/70‑1 (MQ=255)
        gTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                                <  1:747071/71‑1 (MQ=255)
          gCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcc                                                                 <  1:2590620/68‑1 (MQ=255)
             cGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtc                                                           <  1:168697/71‑1 (MQ=255)
               cTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtctc                                                         <  1:1356289/71‑1 (MQ=255)
                tATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                                       >  1:1171821/1‑56 (MQ=255)
                  tctgtcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCg                                                       <  1:3323946/70‑1 (MQ=255)
                  tctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                        >  1:5193665/1‑69 (MQ=255)
                   ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                                       <  1:3387117/53‑1 (MQ=255)
                   ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                                       <  1:1298291/53‑1 (MQ=255)
                   ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                                       <  1:5300695/53‑1 (MQ=255)
                   ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                                      >  1:5165428/1‑70 (MQ=255)
                    tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                                >  1:725684/1‑59 (MQ=255)
                    tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                                >  1:1252817/1‑59 (MQ=255)
                    tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                                >  1:687189/1‑59 (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)
                        gctgctTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTg                                                                      >  1:3583570/1‑49 (MQ=16)
                        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)
                                       ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                                >  1:1414845/1‑40 (MQ=255)
                                       ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                                >  1:4451249/1‑40 (MQ=255)
                                       ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                                    <  1:5615593/52‑1 (MQ=255)
                                        tGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc                                                                  >  1:631224/1‑37 (MQ=255)
                                        tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                        >  1:4719233/1‑47 (MQ=255)
                                        tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                        >  1:1723349/1‑47 (MQ=255)
                                        tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                        >  1:1707352/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:2239762/1‑58 (MQ=255)
                                            tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                                         >  1:5256285/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:1396551/1‑71 (MQ=255)
                                                    tACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATc                    >  1:1283628/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:2842309/42‑1 (MQ=255)
                                                               cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                                      <  1:3609614/42‑1 (MQ=255)
                                                                 gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa        >  1:3097322/1‑70 (MQ=255)
                                                                 gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa        >  1:180735/1‑70 (MQ=255)
                                                                 gCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGa                                >  1:3131241/1‑46 (MQ=255)
                                                                 gCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAAt                     >  1:5092699/1‑57 (MQ=255)
                                                                 gCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAAt                     >  1:3352922/1‑57 (MQ=255)
                                                                 gCGACATGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGACTGGTa        >  1:2420630/1‑70 (MQ=255)
                                                                  cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTg           >  1:2286052/1‑66 (MQ=255)
                                                                  cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc       <  1:1238824/70‑1 (MQ=255)
                                                                  cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc       <  1:3566441/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)
                                                                      |                                                                      
CGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTC  >  Exported/832142‑832282

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: