Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 837,738 C→T 11.8% L71L (CTG→TTG) ‡ narI → nitrate reductase 1, gamma (cytochrome b(NR)) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported837,7380CT11.8% 152.6 / 9.6 59L71L (CTG→TTG) ‡narInitrate reductase 1, gamma (cytochrome b(NR)) subunit
Reads supporting (aligned to +/- strand):  ref base C (25/27);  new base T (2/5);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTG  >  Exported/837671‑837807
                                                                   |                                                                     
gAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa                                                                    >  1:4804580/1‑71 (MQ=255)
    cTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgcc                                                                <  1:1498757/71‑1 (MQ=255)
     tGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg                                                               >  1:1827984/1‑71 (MQ=255)
          tCGAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTg                                                          <  1:5667574/71‑1 (MQ=255)
            gAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGa                                                        >  1:1772263/1‑71 (MQ=255)
             aatcTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGAt                                                       >  1:4149972/4‑71 (MQ=255)
              atcTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACt                                                           <  1:4370513/64‑1 (MQ=255)
              atcTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATg                                                      <  1:4836129/69‑1 (MQ=255)
              atcTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATg                                                      <  1:3530681/69‑1 (MQ=255)
                 tGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTAt                                                   >  1:3813579/1‑71 (MQ=255)
                 tGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTAt                                                   >  1:3413289/1‑71 (MQ=255)
                  gTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGt                                                     <  1:1693949/68‑1 (MQ=255)
                       atatCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACCCCGCACTGGATg                                                      <  1:1901880/62‑1 (MQ=255)
                        tatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                                      <  1:1612111/45‑1 (MQ=255)
                        tatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                                      <  1:625380/45‑1 (MQ=255)
                        tatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGcc                                             <  1:2106917/70‑1 (MQ=255)
                           cGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGc                                         >  1:4708655/1‑71 (MQ=255)
                               aTTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGt                                                     <  1:3749098/55‑1 (MQ=255)
                               aTTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGcc                                     <  1:3701736/71‑1 (MQ=255)
                                ttCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                                      <  1:5320668/37‑1 (MQ=255)
                                ttCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                                      <  1:208765/37‑1 (MQ=255)
                                ttCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                                      <  1:3504586/37‑1 (MQ=255)
                                ttCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGcc                                     <  1:4949025/70‑1 (MQ=255)
                                  cTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                          >  1:2619851/1‑47 (MQ=255)
                                  cTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAt                                  >  1:974300/1‑71 (MQ=255)
                                  cTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAt                                  >  1:365107/1‑71 (MQ=255)
                                  cTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGACGCACTGGATGTATGAAGCCTGGCTGCCGAt                                  >  1:5258217/1‑71 (MQ=255)
                                     gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTg                                <  1:878421/70‑1 (MQ=255)
                                     gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGa                               <  1:3201461/71‑1 (MQ=255)
                                      ggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTAt                                                   <  1:15432/50‑1 (MQ=255)
                                       gATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAg                             >  1:3198810/1‑71 (MQ=255)
                                        aTTTTTGTCGGTCACTTCTTCGGTATGCTGAcgcc                                                                >  1:3579565/1‑35 (MQ=255)
                                        aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                          >  1:4697662/1‑41 (MQ=255)
                                        aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                          >  1:3820107/1‑41 (MQ=255)
                                        aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                            <  1:488473/71‑1 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                          <  1:2713587/39‑1 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                          <  1:4893634/39‑1 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                          <  1:4550046/39‑1 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAg                             >  1:174615/1‑68 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                            >  1:3932177/1‑69 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTg                           >  1:24118/1‑70 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTg                           <  1:3362383/70‑1 (MQ=255)
                                                cGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGa                                   >  1:3958079/1‑56 (MQ=255)
                                                cGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGa                                   >  1:1508084/1‑56 (MQ=255)
                                                cGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGa                                   >  1:5669370/1‑56 (MQ=255)
                                                 ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTg                                           >  1:3190265/1‑47 (MQ=255)
                                                 ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGa                    >  1:5684658/1‑70 (MQ=255)
                                                 ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGa                    >  1:2060130/1‑70 (MQ=255)
                                                 ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGa                    >  1:3720324/1‑70 (MQ=255)
                                                   tCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACATaaaa                 <  1:2201898/71‑1 (MQ=255)
                                                     aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAtt                                 <  1:5530285/53‑1 (MQ=255)
                                                       ttcttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaaa                 <  1:3202934/67‑1 (MQ=255)
                                                         cttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCa            <  1:868563/70‑1 (MQ=255)
                                                            cGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAg                             >  1:2293735/1‑50 (MQ=255)
                                                            cGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCaa           <  1:5480955/68‑1 (MQ=255)
                                                            cGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGt        >  1:1658256/1‑71 (MQ=255)
                                                                 tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa                        <  1:5179832/50‑1 (MQ=255)
                                                                 tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa                        <  1:3684330/50‑1 (MQ=255)
                                                                 tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGc    <  1:4568202/70‑1 (MQ=255)
                                                                   cTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTg     >  1:5624736/1‑67 (MQ=255)
                                                                   cTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTg     >  1:1058779/1‑67 (MQ=255)
                                                                   cTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtg  <  1:4619851/70‑1 (MQ=255)
                                                                   cTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtg  <  1:319788/70‑1 (MQ=255)
                                                                   |                                                                     
GAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTG  >  Exported/837671‑837807

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: