Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 832,692 A→T 11.3% M384L (ATG→TTG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported832,6920AT11.3% 165.6 / 12.0 62M384L (ATG→TTG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base T (3/4);  total (31/31)
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.97e-01

GGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGA  >  Exported/832625‑832762
                                                                   |                                                                      
ggTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATgg                                                                     >  1:1609186/1‑71 (MQ=255)
 gTTGATGCGCTGGGCCAGTAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAAt                                                                       <  1:2994825/68‑1 (MQ=255)
 gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGaaa                                                                        <  1:3994189/67‑1 (MQ=255)
 gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGaaa                                                                        <  1:4205509/67‑1 (MQ=255)
 gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAAt                                                                       <  1:3972065/68‑1 (MQ=255)
 gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAAt                                                                       <  1:2212994/68‑1 (MQ=255)
 gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATgg                                                                     <  1:4221993/70‑1 (MQ=255)
 gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATg                                                                      <  1:2192330/69‑1 (MQ=255)
     aTGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTgcg                                                                >  1:1871332/1‑71 (MQ=255)
         gcTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCg                                                             >  1:2032484/1‑70 (MQ=255)
             ggCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACgg                                                        >  1:2288422/1‑71 (MQ=255)
                 aGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCg                                                             <  1:1852836/62‑1 (MQ=255)
                 aGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCt                                                    <  1:486419/71‑1 (MQ=255)
                 aGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACAAATGGCGAAATGGTTGCGCCGAACGGTTc                                                     >  1:1381824/1‑70 (MQ=255)
                  ggAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGt                                                                    >  1:797062/1‑54 (MQ=255)
                  ggAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTgc                                                                 <  1:612997/57‑1 (MQ=255)
                   gAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGt                                                                    >  1:287084/1‑53 (MQ=255)
                    aaaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTgcgc                                                               >  1:781686/1‑57 (MQ=255)
                    aaaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTgcgc                                                               >  1:2697215/1‑57 (MQ=255)
                    aaaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTAtt                                                 >  1:497840/1‑71 (MQ=255)
                     aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTAt                                                  >  1:5661294/1‑69 (MQ=255)
                     aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTAt                                                  >  1:5679711/1‑69 (MQ=255)
                     aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTg                                                >  1:4676397/1‑71 (MQ=255)
                     aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTg                                                >  1:3755832/1‑71 (MQ=255)
                      aaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTgg                                               >  1:309927/1‑71 (MQ=255)
                      aaCAATCCGCAGTGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGAtcg                                                >  1:1829555/1‑68 (MQ=255)
                      aaCAATCCGCAGTGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGAtc                                                 >  1:1287352/1‑68 (MQ=255)
                           tCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGt                                                                    <  1:690011/45‑1 (MQ=255)
                              ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACgg                                                        <  1:222466/54‑1 (MQ=255)
                              ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACgg                                                        <  1:1353414/54‑1 (MQ=255)
                              ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACgg                                                        <  1:2235080/54‑1 (MQ=255)
                              ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCt                                       <  1:3559558/71‑1 (MQ=255)
                                agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                                     <  1:5462661/69‑1 (MQ=255)
                                agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                                     <  1:4568169/69‑1 (MQ=255)
                                agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                                     <  1:2090708/69‑1 (MQ=255)
                                aaTGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTgg                                     <  1:1443698/71‑1 (MQ=255)
                                 aTGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCt                                                    <  1:2674713/55‑1 (MQ=255)
                                  tGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAAc                                                          <  1:1467036/48‑1 (MQ=255)
                                  tGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTAtt                                                 <  1:3854500/57‑1 (MQ=255)
                                   ggAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTc                                                     <  1:5639540/52‑1 (MQ=255)
                                    gAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTgg                                               >  1:4543940/1‑57 (MQ=255)
                                    gAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGc                                  >  1:3160069/1‑70 (MQ=255)
                                      aaaCTGTCGCCTTTAATACCAATGGCGAAATGGtt                                                                   >  1:2026243/1‑35 (MQ=255)
                                       aaCTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTggg                                    <  1:4571973/65‑1 (MQ=255)
                                        aCTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                                        >  1:2670162/1‑60 (MQ=255)
                                              gCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCaa                       <  1:3167222/71‑1 (MQ=255)
                                              gCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCa                        >  1:3800650/1‑70 (MQ=255)
                                              gCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCa                        >  1:2257604/1‑70 (MQ=255)
                                               ccTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTc                                                     <  1:888479/40‑1 (MQ=255)
                                                   tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                                        >  1:5345206/1‑49 (MQ=255)
                                                   tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                                        >  1:1573237/1‑49 (MQ=255)
                                                   tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                                        >  1:4182093/1‑49 (MQ=255)
                                                   tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATgg                   >  1:4799910/1‑70 (MQ=255)
                                                      tACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCg                                 <  1:1178768/53‑1 (MQ=255)
                                                      tACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATc               <  1:4637853/71‑1 (MQ=255)
                                                       aCCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCg                                 <  1:2057519/52‑1 (MQ=255)
                                                         cAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCtt             <  1:1573428/70‑1 (MQ=255)
                                                         cAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTg            <  1:4891240/71‑1 (MQ=255)
                                                         cAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTg            <  1:1672421/71‑1 (MQ=255)
                                                              gtgAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAATCt              <  1:3070023/62‑1 (MQ=255)
                                                              gCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATg                    >  1:2187086/1‑58 (MQ=255)
                                                              gCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATg                    >  1:148097/1‑58 (MQ=255)
                                                              gCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATg                    >  1:4714306/1‑58 (MQ=255)
                                                                gAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTGGGGAGaaaaa                            >  1:5324752/1‑47 (MQ=16)
                                                                 aaaTGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTggg                                    >  1:1987587/1‑39 (MQ=255)
                                                                 aaaTGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTggg                                    >  1:4203655/1‑39 (MQ=255)
                                                                   aTGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGaa                             >  1:76384/1‑44 (MQ=255)
                                                                   aTGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGa  <  1:1441733/71‑1 (MQ=255)
                                                                   |                                                                      
GGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGA  >  Exported/832625‑832762

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: