Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 834,249 A→C 26.7% I903L (ATC→CTC)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported834,2490AC26.7% 68.4 / 44.3 60I903L (ATC→CTC) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (16/28);  new base C (5/11);  total (21/39)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

AGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGA  >  Exported/834179‑834319
                                                                      |                                                                      
aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa                                                                        >  1:4131261/1‑71 (MQ=255)
    aaaaCCTCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                                    >  1:3373877/1‑71 (MQ=255)
    aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                                    <  1:537943/71‑1 (MQ=255)
    aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                                    >  1:2028724/1‑71 (MQ=255)
    aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                                    >  1:4059875/1‑71 (MQ=255)
    aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                                    >  1:3264816/1‑71 (MQ=255)
          gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCt                                                              <  1:151129/71‑1 (MQ=255)
                 acaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg                                                          >  1:2989700/1‑68 (MQ=255)
                  caTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa                                                                        >  1:3775983/1‑53 (MQ=255)
                    ttATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGaaa                                                    <  1:276779/71‑1 (MQ=255)
                       tGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTAt                                                                       >  1:1417611/1‑49 (MQ=255)
                         gTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                                               <  1:512139/71‑1 (MQ=255)
                         gTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                                               <  1:2863443/71‑1 (MQ=255)
                          tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGGGAAAATCGGt                                              <  1:918393/71‑1 (MQ=255)
                          tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                              <  1:2213556/71‑1 (MQ=255)
                          tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                              <  1:4618116/71‑1 (MQ=255)
                            gTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                            >  1:2423121/1‑71 (MQ=255)
                              agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg                                          >  1:2292926/1‑71 (MQ=255)
                              agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg                                          >  1:1718139/1‑71 (MQ=255)
                               gagCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATgg                                         <  1:2331647/71‑1 (MQ=255)
                               gagCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATgg                                         <  1:2770130/71‑1 (MQ=255)
                                agCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGc                                        <  1:2732854/71‑1 (MQ=255)
                                  cgcgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAAt                                                  <  1:2180099/59‑1 (MQ=255)
                                    cgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt                                     >  1:2191912/1‑70 (MQ=255)
                                      aTTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                                               <  1:3378964/58‑1 (MQ=255)
                                       ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGc                                                                   >  1:3829941/1‑37 (MQ=255)
                                       ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                            >  1:5053776/1‑60 (MQ=255)
                                       ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt                                     >  1:2175296/1‑67 (MQ=255)
                                       ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAg                                 >  1:4055910/1‑71 (MQ=255)
                                         aTCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCt                                                              <  1:4876335/40‑1 (MQ=255)
                                         aTCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAgg                                <  1:3738083/70‑1 (MQ=255)
                                            cGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAgg                                <  1:4685043/67‑1 (MQ=255)
                                             ggCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTg                           <  1:2770132/71‑1 (MQ=255)
                                                 aCTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGAAAGGGATTGCCTgg                      <  1:419840/71‑1 (MQ=255)
                                                   ttACAAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAggg                               <  1:3928344/61‑1 (MQ=255)
                                                     aCGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGaa                    <  1:291860/70‑1 (MQ=255)
                                                      tgAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAAt                                           <  1:271263/45‑1 (MQ=255)
                                                        aaCGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTgg                      <  1:267508/65‑1 (MQ=255)
                                                        aaCGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc                >  1:4518783/1‑71 (MQ=255)
                                                          cGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                              <  1:1646611/39‑1 (MQ=255)
                                                          cGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                              <  1:3049915/39‑1 (MQ=255)
                                                          cGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAggg                               <  1:584154/54‑1 (MQ=255)
                                                           gCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGa                              <  1:2556674/54‑1 (MQ=255)
                                                           gCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGa                              <  1:2753986/54‑1 (MQ=255)
                                                           gCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTgg                      <  1:3178416/62‑1 (MQ=255)
                                                            cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                            >  1:2996780/1‑39 (MQ=255)
                                                            cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                            >  1:3443090/1‑39 (MQ=255)
                                                            cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                            >  1:4420579/1‑39 (MQ=255)
                                                            cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGCTAATGGCGGTAAAGGGATTGCCTGGAACACCCaga            <  1:761213/71‑1 (MQ=255)
                                                             tttACCTCTATCGGCCCGCGGATGGAGAAAATCGGTa                                             >  1:2821064/1‑37 (MQ=255)
                                                              ttACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGc                          <  1:1780629/55‑1 (MQ=255)
                                                              ttACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGc                          <  1:2732616/55‑1 (MQ=255)
                                                               tACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTa                                    >  1:2957626/1‑44 (MQ=255)
                                                               tACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGc          <  1:3115838/70‑1 (MQ=255)
                                                               tACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATACCCaggat          <  1:4847850/70‑4 (MQ=255)
                                                                aCCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCg         <  1:5329022/70‑1 (MQ=255)
                                                                 ccTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga        <  1:2660199/70‑1 (MQ=255)
                                                                 ccTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga        <  1:3160474/70‑1 (MQ=255)
                                                                  ctctATCGGCCCGCTGATGGAGAAAATCGGTAATgg                                         <  1:2364127/36‑1 (MQ=255)
                                                                  ctctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga      <  1:4396547/71‑1 (MQ=255)
                                                                    ctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGc          <  1:210801/65‑1 (MQ=255)
                                                                     tATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATg    >  1:1766222/1‑70 (MQ=255)
                                                                      aTCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGa  >  1:894272/1‑71 (MQ=255)
                                                                      aTCGGCCCGCTGATGGAGAAAATAGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGa  >  1:4372801/1‑71 (MQ=255)
                                                                      |                                                                      
AGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGA  >  Exported/834179‑834319

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: