Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 836,443 A→G 32.1% V387V (GTA→GTG 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,4430AG32.1% 68.8 / 42.7 56V387V (GTA→GTGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base A (21/17);  new base G (6/12);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCT  >  Exported/836380‑836511
                                                               |                                                                    
agACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATc                                                               >  1:1267636/1‑71 (MQ=255)
agACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATc                                                               >  1:57277/1‑71 (MQ=255)
agACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATc                                                               >  1:2603199/1‑71 (MQ=255)
 gACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCt                                                              <  1:5278793/71‑1 (MQ=255)
 gACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCt                                                              <  1:56326/71‑1 (MQ=255)
         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)
                 gTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTa                                                                      >  1:4066757/1‑47 (MQ=255)
                 gTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTa                                                                      >  1:3994180/1‑47 (MQ=255)
                 gTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAcc                                              >  1:1523098/1‑71 (MQ=255)
                    gCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccgcc                                           >  1:2224747/1‑71 (MQ=255)
                      aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGa                                                         <  1:2640002/55‑1 (MQ=255)
                      aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGa                                                         <  1:136734/55‑1 (MQ=255)
                      aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGa                                                         <  1:5342355/55‑1 (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)
                                         tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                        <  1:2988772/37‑1 (MQ=21)
                                         tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                        <  1:5196247/37‑1 (MQ=21)
                                         tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                        <  1:1288919/37‑1 (MQ=21)
                                         tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                        <  1:5113463/37‑1 (MQ=21)
                                         tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                        <  1:210936/37‑1 (MQ=21)
                                         tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATatgtt                                                  <  1:5129562/43‑3 (MQ=16)
                                         tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGag                                               >  1:5576633/1‑45 (MQ=16)
                                         tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCggc                                        <  1:4324141/53‑2 (MQ=255)
                                         tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCggc                                        <  1:4330503/53‑2 (MQ=255)
                                          cGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgc                     >  1:3698793/1‑71 (MQ=255)
                                          cGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgc                     >  1:3473738/1‑71 (MQ=255)
                                           gAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                        >  1:3710952/1‑35 (MQ=16)
                                           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)
                                                               |                                                                    
AGACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCT  >  Exported/836380‑836511

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: