Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 832,205 C→T 16.2% Y221Y (TAC→TAT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported832,2050CT16.2% 187.4 / 24.7 80Y221Y (TAC→TATnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base T (4/9);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACT  >  Exported/832135‑832272
                                                                      |                                                                   
tttCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTAc                                                                     >  1:586300/1‑71 (MQ=255)
  tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACt                                                                    <  1:4581840/70‑1 (MQ=255)
  tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACt                                                                    <  1:428408/70‑1 (MQ=255)
  tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTg                                                                   >  1:3910588/1‑71 (MQ=255)
  tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTg                                                                   >  1:3482011/1‑71 (MQ=255)
      aCGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGc                                                                  <  1:4645435/68‑1 (MQ=255)
       cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACt                                                              <  1:4281585/71‑1 (MQ=255)
        gcttCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGc                                                                  <  1:798671/63‑1 (MQ=255)
        gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                             <  1:3784066/71‑1 (MQ=255)
          aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGAc                                                               >  1:3373802/1‑67 (MQ=255)
          aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGAc                                                               >  1:2123407/1‑67 (MQ=255)
          aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGc                                                           >  1:3964838/1‑71 (MQ=255)
           tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc                                                          <  1:121244/71‑1 (MQ=255)
           tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc                                                          <  1:4910086/71‑1 (MQ=255)
           tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc                                                          <  1:2835482/71‑1 (MQ=255)
            cGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACt                                                                    <  1:2174563/60‑1 (MQ=255)
             gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc                                                        >  1:4599508/1‑71 (MQ=255)
             gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcct                                                         <  1:4498985/70‑1 (MQ=255)
               gTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                      <  1:747071/71‑1 (MQ=255)
                 gCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcc                                                       <  1:2590620/68‑1 (MQ=255)
                    cGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtc                                                 <  1:168697/71‑1 (MQ=255)
                      cTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtctc                                               <  1:1356289/71‑1 (MQ=255)
                       tATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                             >  1:1171821/1‑56 (MQ=255)
                        aTCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGAc                                                               >  1:3205980/1‑53 (MQ=19)
                        aTCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGAc                                                               >  1:3175116/1‑53 (MQ=19)
                         tctgtcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCg                                             <  1:3323946/70‑1 (MQ=255)
                         tctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                              >  1:5193665/1‑69 (MQ=255)
                          ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                             <  1:3387117/53‑1 (MQ=255)
                          ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                             <  1:5300695/53‑1 (MQ=255)
                          ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                             <  1:1298291/53‑1 (MQ=255)
                          ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                            >  1:5165428/1‑70 (MQ=255)
                           tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                      >  1:687189/1‑59 (MQ=255)
                           tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                      >  1:725684/1‑59 (MQ=255)
                           tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                      >  1:1252817/1‑59 (MQ=255)
                            ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                              >  1:1526610/1‑66 (MQ=255)
                            ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                              >  1:1235824/1‑66 (MQ=255)
                            ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaaa                                         >  1:5052580/1‑71 (MQ=255)
                             tcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCgatg                                          <  1:761100/69‑4 (MQ=255)
                             tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAc                                        >  1:3511473/1‑71 (MQ=255)
                               gctgctTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTg                                                            >  1:3583570/1‑49 (MQ=16)
                               gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                                       <  1:2597733/70‑1 (MQ=255)
                               gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                                       <  1:569964/70‑1 (MQ=255)
                                  gATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                                       <  1:3422109/67‑1 (MQ=255)
                                    ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTggg                                   >  1:5620596/1‑69 (MQ=255)
                                    ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGc                                 <  1:1200017/71‑1 (MQ=255)
                                      ggCGGTACATGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCg                                <  1:1009608/70‑1 (MQ=255)
                                          gTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCaa                           <  1:3585810/71‑1 (MQ=255)
                                            aCTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                            <  1:663507/68‑1 (MQ=255)
                                             cTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAAc                         <  1:4644879/70‑1 (MQ=255)
                                             cTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACt                        <  1:163062/71‑1 (MQ=255)
                                              ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                      >  1:1414845/1‑40 (MQ=255)
                                              ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                      >  1:4451249/1‑40 (MQ=255)
                                              ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                          <  1:5615593/52‑1 (MQ=255)
                                               tGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc                                                        >  1:631224/1‑37 (MQ=255)
                                               tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                              >  1:4719233/1‑47 (MQ=255)
                                               tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                              >  1:1723349/1‑47 (MQ=255)
                                               tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                              >  1:1707352/1‑47 (MQ=255)
                                                 cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACg                    <  1:1722190/71‑1 (MQ=255)
                                                   tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGt                   <  1:3382378/70‑1 (MQ=255)
                                                   tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGt                   <  1:3375919/70‑1 (MQ=255)
                                                   tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                               >  1:5256285/1‑58 (MQ=255)
                                                   tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                               >  1:2239762/1‑58 (MQ=255)
                                                    aaGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTa                  <  1:1006427/70‑1 (MQ=255)
                                                    aaGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                            >  1:3515463/1‑44 (MQ=255)
                                                    aaGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGt                   <  1:4166571/69‑1 (MQ=255)
                                                       tttCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGc                                 <  1:3806603/51‑1 (MQ=255)
                                                         tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGaa            >  1:1655704/1‑71 (MQ=255)
                                                         tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGa             <  1:3303868/70‑1 (MQ=255)
                                                          cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa             <  1:999977/69‑1 (MQ=255)
                                                          cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa             <  1:2123804/69‑1 (MQ=255)
                                                          cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt           >  1:3042141/1‑71 (MQ=255)
                                                           tACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt           >  1:1807426/1‑70 (MQ=255)
                                                           tACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATc          >  1:1396551/1‑71 (MQ=255)
                                                           tACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATc          >  1:1283628/1‑71 (MQ=255)
                                                            aCGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt           >  1:4469270/1‑69 (MQ=255)
                                                             cGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATc          <  1:5684819/69‑1 (MQ=255)
                                                             cGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt           <  1:90272/68‑1 (MQ=255)
                                                             cGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATct         >  1:4196601/1‑70 (MQ=255)
                                                                cTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATctgctg     <  1:1065204/71‑1 (MQ=255)
                                                                 tGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCgg              <  1:1009180/61‑1 (MQ=255)
                                                                  ggTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGAc   >  1:5140846/1‑71 (MQ=255)
                                                                      cTGCGTCTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACt  <  1:1194552/68‑1 (MQ=255)
                                                                      cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                               >  1:3899584/1‑39 (MQ=255)
                                                                      cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                            <  1:3609614/42‑1 (MQ=255)
                                                                      cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                            <  1:2842309/42‑1 (MQ=255)
                                                                      |                                                                   
TTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACT  >  Exported/832135‑832272

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: