Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 836,455 G→C 31.8% A391A (GCG→GCC 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,4550GC31.8% 86.9 / 52.6 66A391A (GCG→GCCnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (23/22);  new base C (9/12);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-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:2075581/70‑1 (MQ=255)
       ggTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGa                                                       <  1:423820/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)
             aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGa                                                                <  1:136734/55‑1 (MQ=255)
             aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGa                                                                <  1:2640002/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:2060746/1‑71 (MQ=255)
                               gTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACt                              >  1:2737099/1‑71 (MQ=255)
                                tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                               <  1:210936/37‑1 (MQ=21)
                                tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                               <  1:5196247/37‑1 (MQ=21)
                                tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                               <  1:1288919/37‑1 (MQ=21)
                                tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                               <  1:2988772/37‑1 (MQ=21)
                                tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                               <  1:5113463/37‑1 (MQ=21)
                                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: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:1161640/1‑71 (MQ=255)
                                                 ccGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTAt            >  1:4710526/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:4812093/1‑52 (MQ=255)
                                                         gTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                       >  1:2672883/1‑52 (MQ=255)
                                                         gTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                       >  1:603611/1‑52 (MQ=255)
                                                          tATCTGGCGAATCTGCTGACCGCCGGTGATACCAAAc                                     <  1:4700042/37‑1 (MQ=255)
                                                          tATCTGGCGAATCTGCTGACCGCCGGTGATACCAAAc                                     <  1:3942951/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:1838486/51‑1 (MQ=255)
                                                              tGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                   <  1:1812489/51‑1 (MQ=255)
                                                              tGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                   <  1:727777/51‑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: