Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 834,842 C→A 34.9% S1100S (TCC→TCA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported834,8420CA34.9% 87.3 / 48.2 63S1100S (TCC→TCAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (18/23);  new base A (8/14);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

ATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCC  >  Exported/834776‑834911
                                                                  |                                                                     
aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCAcc                                                                    >  1:4653409/1‑70 (MQ=255)
   cAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTa                                                                  >  1:5375499/1‑69 (MQ=255)
      cGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAg                                                               >  1:2368751/1‑69 (MQ=255)
      cGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGc                                                              >  1:4637303/1‑70 (MQ=255)
       ggCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGa                                                            >  1:2051331/1‑71 (MQ=255)
         cAACCAGGAAAAAGCGCTCAACTTTCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa                                                          <  1:3852092/71‑1 (MQ=255)
          aaCCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa                                                          <  1:2695136/70‑1 (MQ=255)
             ccggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAcc                                                       <  1:4368715/68‑1 (MQ=255)
             cAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACaa                                                         >  1:1402299/1‑68 (MQ=255)
              cggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAc                                                        <  1:3157486/67‑1 (MQ=255)
              cggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAc                                                        <  1:567457/67‑1 (MQ=255)
              cggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAc                                                        <  1:189225/67‑1 (MQ=255)
              aGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctg                                                     >  1:3144394/1‑71 (MQ=255)
                gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct                                                   <  1:1125150/71‑1 (MQ=255)
                gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct                                                   <  1:916451/71‑1 (MQ=255)
                gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct                                                   <  1:3603640/71‑1 (MQ=255)
                gAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTAGAGTGAAAACct                                                      >  1:2600314/1‑68 (MQ=255)
                gAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAgtg                                                             <  1:2542890/61‑3 (MQ=255)
                  aaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctga                                                 >  1:4887880/1‑71 (MQ=255)
                   aaaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaat                                                >  1:4634689/1‑71 (MQ=255)
                    aaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatg                                               <  1:3896675/71‑1 (MQ=255)
                    aaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatg                                               <  1:5477393/71‑1 (MQ=255)
                     aGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatgc                                              >  1:305566/1‑71 (MQ=255)
                      gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactccac                                                                     >  1:1098225/1‑47 (MQ=255)
                      gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactccac                                                                     >  1:839167/1‑47 (MQ=255)
                      gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactccac                                                                     >  1:5211438/1‑47 (MQ=255)
                      gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactccac                                                                     >  1:5217717/1‑47 (MQ=255)
                      gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactcca                                                                      >  1:1030869/1‑46 (MQ=255)
                      gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                                             <  1:778149/71‑1 (MQ=255)
                      gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                                             <  1:3647857/71‑1 (MQ=255)
                      gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                                             <  1:4230928/71‑1 (MQ=255)
                       cgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                                             >  1:3873271/1‑70 (MQ=255)
                         cttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAc                                          >  1:132462/4‑71 (MQ=255)
                         cttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAc                                          >  1:3654334/4‑71 (MQ=255)
                          ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACcttct                                                   <  1:2322660/59‑1 (MQ=255)
                          ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcg                                         >  1:2673952/3‑70 (MQ=255)
                          ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcg                                         >  1:2037702/3‑70 (MQ=255)
                            aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGCCtt                                        >  1:456030/1‑70 (MQ=255)
                            aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt                                       <  1:5208427/71‑1 (MQ=255)
                            aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt                                       <  1:3100376/71‑1 (MQ=255)
                            aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt                                       <  1:2071783/71‑1 (MQ=255)
                            aaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct                                       >  1:270191/1‑68 (MQ=255)
                             aCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACtt                                        >  1:1076048/1‑69 (MQ=255)
                             aCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACtt                                        >  1:1695595/1‑69 (MQ=255)
                              cTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAg                                     >  1:5416809/1‑71 (MQ=255)
                                tCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGGCa                                                          <  1:2015197/48‑1 (MQ=255)
                                tCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGt                                   >  1:3679281/1‑71 (MQ=255)
                                tCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGaaa                                                          <  1:4457032/48‑3 (MQ=19)
                                  cTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTcg                                 <  1:1511231/71‑1 (MQ=255)
                                    gACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgctaa                                                 <  1:3006596/53‑3 (MQ=255)
                                     aCGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcgg                              <  1:4046506/71‑1 (MQ=255)
                                     aCGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct                                       <  1:1805862/62‑4 (MQ=255)
                                      cgccgcACCAGAAATGGGGCATTCAGTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct                                       <  1:1920594/61‑4 (MQ=255)
                                       gccgcACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgct                                                   <  1:1119767/48‑1 (MQ=19)
                                       gccgcACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgct                                                   <  1:1606080/48‑1 (MQ=19)
                                        ccgcACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtgg                           <  1:4735454/71‑1 (MQ=255)
                                         cgcACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggt                          >  1:2559416/1‑71 (MQ=255)
                                          gcACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtc                         >  1:4439869/1‑71 (MQ=255)
                                            aCCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtccg                       <  1:368621/71‑1 (MQ=255)
                                              cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg                                               <  1:3101428/45‑1 (MQ=255)
                                              cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg                                               <  1:434018/45‑1 (MQ=255)
                                              cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg                                               <  1:5671764/45‑1 (MQ=255)
                                               aGAAGTAGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcgg                              <  1:3939067/61‑1 (MQ=255)
                                                gAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAgg                                    >  1:2987148/1‑54 (MQ=255)
                                                   gTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCt                >  1:1825984/1‑71 (MQ=255)
                                                    tGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                                             <  1:5218110/41‑1 (MQ=255)
                                                         tATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtc                         <  1:4076261/56‑1 (MQ=255)
                                                          aTCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGtt            <  1:825208/68‑1 (MQ=255)
                                                            ccactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtccg                       <  1:3825328/55‑1 (MQ=255)
                                                             cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                      >  1:4400425/1‑39 (MQ=255)
                                                             cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                      >  1:3318768/1‑39 (MQ=255)
                                                             cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                      >  1:4990428/1‑39 (MQ=255)
                                                             cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                      >  1:4188080/1‑39 (MQ=255)
                                                             cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                      >  1:5580994/1‑39 (MQ=255)
                                                             cactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTcg                                 >  1:5123155/1‑44 (MQ=255)
                                                             cactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTg       <  1:4499858/70‑1 (MQ=255)
                                                             cactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTg       >  1:2545110/1‑70 (MQ=255)
                                                             cactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTg       >  1:4805053/1‑70 (MQ=255)
                                                                tccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGt                                   >  1:1228432/1‑39 (MQ=255)
                                                                  cacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGcc  <  1:5009583/70‑1 (MQ=255)
                                                                  cacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGcc  <  1:4441411/70‑1 (MQ=255)
                                                                  cacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGcc  <  1:3818974/70‑1 (MQ=255)
                                                                  |                                                                     
ATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCC  >  Exported/834776‑834911

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: