Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 833,868 C→A 13.3% L776M (CTG→ATG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported833,8680CA13.3% 220.9 / 13.8 75L776M (CTG→ATG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (38/27);  new base A (5/5);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

AGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACA  >  Exported/833801‑833934
                                                                   |                                                                  
agaagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTg                                                                  >  1:5303677/1‑70 (MQ=255)
agaagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt                                                                 >  1:4526262/1‑71 (MQ=255)
agaagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt                                                                 >  1:4808825/1‑71 (MQ=255)
agaagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt                                                                 >  1:4832198/1‑71 (MQ=255)
agaagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt                                                                 >  1:2941702/1‑71 (MQ=255)
agaagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt                                                                 <  1:5625951/71‑1 (MQ=255)
     tGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt                                                                 <  1:1380689/66‑1 (MQ=255)
     tGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt                                                                 <  1:3412403/66‑1 (MQ=255)
                   aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGt                                                                 <  1:443498/52‑1 (MQ=255)
                   aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTc                                              >  1:1998311/1‑71 (MQ=255)
                    aTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTcc                                             <  1:5286047/71‑1 (MQ=255)
                      ggTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGc                                                            >  1:750484/1‑54 (MQ=255)
                       gTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAc                                          <  1:4024818/71‑1 (MQ=255)
                         attgAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGatat                                        <  1:2025996/68‑1 (MQ=255)
                          tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAg                                                             <  1:3298717/49‑1 (MQ=255)
                          tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAg                                                             <  1:663837/49‑1 (MQ=255)
                          tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg                                            <  1:2182199/66‑1 (MQ=255)
                           ggAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGc                                                            <  1:311481/49‑1 (MQ=255)
                           ggAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGc                                                            <  1:4466279/49‑1 (MQ=255)
                           ggAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca                                      >  1:5251164/1‑71 (MQ=255)
                            gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa                                           >  1:5042738/1‑65 (MQ=255)
                            gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatcat                                     >  1:5199391/1‑71 (MQ=255)
                            gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca                                      <  1:1612153/70‑1 (MQ=255)
                            gAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTcca                                                              >  1:557420/1‑44 (MQ=255)
                            gAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGAtatcg                                      <  1:4205290/70‑2 (MQ=255)
                             aaGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtt                                    >  1:2428842/1‑71 (MQ=255)
                                gCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTg                                 <  1:5113514/71‑1 (MQ=255)
                                 cAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCAc                                                          >  1:796553/1‑45 (MQ=255)
                                   aGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAg                                                             <  1:2521653/40‑1 (MQ=255)
                                   aGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACGTGTCTCTATTCCGACATCAttt                                   >  1:4358453/1‑66 (MQ=255)
                                   aGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGcc                               >  1:4173647/1‑70 (MQ=255)
                                    gCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca                                      >  1:110470/1‑62 (MQ=255)
                                    gCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGa                             >  1:5493622/1‑71 (MQ=255)
                                    gCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGa                             <  1:4925303/71‑1 (MQ=255)
                                     cTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGa                                                              >  1:4072121/1‑37 (MQ=255)
                                     cTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGa                                                              >  1:3182038/1‑37 (MQ=255)
                                       ggATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTACGAcatca                                      >  1:5185177/1‑59 (MQ=255)
                                        gATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGc                         <  1:5279210/71‑1 (MQ=255)
                                        gATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGc                         <  1:2478558/71‑1 (MQ=255)
                                         aTCTGGAGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGc                         >  1:728416/1‑70 (MQ=255)
                                           cTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa                                           >  1:1669180/1‑50 (MQ=255)
                                           cTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa                                           >  1:785298/1‑50 (MQ=255)
                                              gtggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGActt                    <  1:3872112/70‑3 (MQ=255)
                                                 gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtt                                    <  1:983705/51‑1 (MQ=255)
                                                 gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtt                                    <  1:3865070/51‑1 (MQ=255)
                                                  gtACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGcc                               <  1:3076230/54‑1 (MQ=255)
                                                    aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACg              <  1:1473821/70‑1 (MQ=255)
                                                    aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACg              >  1:2414707/1‑70 (MQ=255)
                                                    aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACg              <  1:2530617/70‑1 (MQ=255)
                                                    aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACg              <  1:4562829/70‑1 (MQ=255)
                                                     cGCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaa            <  1:774887/71‑1 (MQ=255)
                                                      gCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaa           <  1:5376024/71‑1 (MQ=255)
                                                      gCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa           <  1:3798898/71‑1 (MQ=255)
                                                        tcgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaaa          >  1:1865983/3‑70 (MQ=255)
                                                         cgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaa            >  1:2054885/2‑67 (MQ=255)
                                                          gACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaaa         <  1:2869982/69‑1 (MQ=255)
                                                          gACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaaa         <  1:227333/69‑1 (MQ=255)
                                                          gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaaaa         <  1:2516655/69‑1 (MQ=255)
                                                          gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAga       >  1:3497426/1‑71 (MQ=255)
                                                          gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAga       >  1:3490634/1‑71 (MQ=255)
                                                           aCTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGcc                               >  1:3902865/1‑46 (MQ=255)
                                                             ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAgacga    >  1:2829301/1‑71 (MQ=255)
                                                               ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc                      >  1:4541119/1‑51 (MQ=255)
                                                               ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa  >  1:3799207/1‑71 (MQ=255)
                                                               ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa  >  1:1729023/1‑71 (MQ=255)
                                                               ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa  >  1:252478/1‑71 (MQ=255)
                                                               ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa  >  1:5451873/1‑71 (MQ=255)
                                                               ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa  >  1:3290762/1‑71 (MQ=255)
                                                               ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa  >  1:3333857/1‑71 (MQ=255)
                                                               ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGCCGACCTGGTACGAAAAAGACGACa  >  1:3451766/1‑71 (MQ=255)
                                                               ccGTCTCTCGAGCACCTGTCTCTATTCCGAcatca                                      >  1:4583525/1‑35 (MQ=255)
                                                                cGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaa           >  1:4675448/1‑61 (MQ=255)
                                                                cGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa  >  1:590478/1‑70 (MQ=255)
                                                                 gTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaag           >  1:5404182/1‑59 (MQ=255)
                                                                 gTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaa          >  1:1978027/1‑61 (MQ=255)
                                                                   |                                                                  
AGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACA  >  Exported/833801‑833934

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: