Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 836,459 C→A 25.4% L393M (CTG→ATG)  narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported836,4590CA25.4% 111.4 / 26.0 55L393M (CTG→ATG) narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (24/17);  new base A (8/6);  total (32/23)
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.94e-01

TGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATG  >  Exported/836389‑836518
                                                                      |                                                           
tGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATc                                                             <  1:4638752/71‑1 (MQ=255)
tGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATc                                                             <  1:3902555/71‑1 (MQ=255)
       ggTAGCAACGGCATTCTGCCAGACGTCGAATGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAc                                                      >  1:2362281/1‑71 (MQ=255)
       ggTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGa                                                       <  1:423820/70‑1 (MQ=255)
       ggTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGa                                                       <  1:2075581/70‑1 (MQ=255)
       ggTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAc                                                      >  1:5250598/1‑71 (MQ=255)
        gTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAcc                                                     >  1:1523098/1‑71 (MQ=255)
           gCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccgcc                                                  >  1:2224747/1‑71 (MQ=255)
             aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCgg                                                >  1:3747422/1‑71 (MQ=255)
              aCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATc                                                             <  1:1697224/57‑1 (MQ=255)
              aCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                               >  1:28786/1‑71 (MQ=255)
              aCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                               <  1:1931873/71‑1 (MQ=255)
                     tCTGCCAGACGTCAAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATAcc                                         >  1:3861722/1‑70 (MQ=255)
                       tGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCaaa                                      >  1:2621459/1‑71 (MQ=255)
                           aGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                                  <  1:2195677/71‑1 (MQ=255)
                            gACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGt                                 <  1:386904/71‑1 (MQ=255)
                               gTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACt                              >  1:2737099/1‑71 (MQ=255)
                               gTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACt                              >  1:2060746/1‑71 (MQ=255)
                                tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATatgtt                                                         <  1:5129562/43‑3 (MQ=16)
                                tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGag                                                      >  1:5576633/1‑45 (MQ=16)
                                tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCggc                                               <  1:4330503/53‑2 (MQ=255)
                                tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCggc                                               <  1:4324141/53‑2 (MQ=255)
                                 cGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgc                            >  1:3473738/1‑71 (MQ=255)
                                 cGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgc                            >  1:3698793/1‑71 (MQ=255)
                                  gAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTa                                >  1:1249682/1‑66 (MQ=255)
                                     aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcg                         >  1:1105263/1‑70 (MQ=255)
                                       ccTGAGTATTCCGGTACAGTATCTGGCGAATCTGCTGa                                                       >  1:3340767/1‑38 (MQ=255)
                                        cTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCt                     <  1:2254244/71‑1 (MQ=255)
                                        cTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACt                     <  1:2687358/71‑1 (MQ=255)
                                           cGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGa                   <  1:4432431/70‑1 (MQ=255)
                                           cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGaa                  <  1:4110428/71‑1 (MQ=255)
                                           cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                   >  1:3630906/1‑70 (MQ=255)
                                            gTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGaa                  >  1:1837881/1‑70 (MQ=255)
                                             tATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCt                     <  1:1188187/66‑1 (MQ=255)
                                             tATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACt                              <  1:1475079/57‑1 (MQ=255)
                                               ttCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACg               <  1:1596027/70‑1 (MQ=255)
                                               ttCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACg               >  1:401318/1‑70 (MQ=255)
                                                tCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTa             >  1:65595/1‑71 (MQ=255)
                                                 ccGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGAt            >  1:1989129/1‑71 (MQ=255)
                                                 ccGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTAt            >  1:4710526/1‑71 (MQ=255)
                                                 ccGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTAt            >  1:1161640/1‑71 (MQ=255)
                                                    gTACAGTATCTGGCGAATCTGCTGACCGCCGGAGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCt         <  1:537934/71‑1 (MQ=255)
                                                      aCAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCt         <  1:1362446/69‑1 (MQ=255)
                                                       caatATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGc      >  1:982427/4‑71 (MQ=255)
                                                        aatATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGc      >  1:5160913/3‑70 (MQ=255)
                                                         gTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                       >  1:2672883/1‑52 (MQ=255)
                                                         gTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                       >  1:4812093/1‑52 (MQ=255)
                                                         gTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                       >  1:603611/1‑52 (MQ=255)
                                                          tATCTGGCGAATCTGCTGACCGCCGGTGATACCAAAc                                     <  1:3942951/37‑1 (MQ=255)
                                                          tATCTGGCGAATCTGCTGACCGCCGGTGATACCAAAc                                     <  1:4700042/37‑1 (MQ=255)
                                                          tATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTa                                >  1:5315952/1‑42 (MQ=255)
                                                           aTCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                   >  1:236465/1‑54 (MQ=255)
                                                           aTCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGCGATg  >  1:295549/1‑71 (MQ=255)
                                                            tCTGGCGAATCTGCTGACCGCCGGTGATACCAAAc                                     <  1:1542139/35‑1 (MQ=255)
                                                              tGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                   <  1:1812489/51‑1 (MQ=255)
                                                              tGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                   <  1:727777/51‑1 (MQ=255)
                                                              tGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                   <  1:1838486/51‑1 (MQ=255)
                                                                   aaTCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGaaa                 >  1:5265714/1‑48 (MQ=255)
                                                                   aaTCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                   <  1:1351984/46‑1 (MQ=255)
                                                                      |                                                           
TGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATG  >  Exported/836389‑836518

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: