Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 832,220 T→C 23.5% P226P (CCT→CCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported832,2200TC23.5% 117.0 / 35.2 68P226P (CCT→CCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (26/26);  new base C (6/10);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

GTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTA  >  Exported/832150‑832288
                                                                      |                                                                    
gTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                                      <  1:747071/71‑1 (MQ=255)
     cGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtc                                                                 <  1:168697/71‑1 (MQ=255)
       cTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtctc                                                               <  1:1356289/71‑1 (MQ=255)
          tctgtcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCg                                                             <  1:3323946/70‑1 (MQ=255)
          tctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                              >  1:5193665/1‑69 (MQ=255)
           ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                                            >  1:5165428/1‑70 (MQ=255)
            tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                                      >  1:725684/1‑59 (MQ=255)
            tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                                      >  1:1252817/1‑59 (MQ=255)
            tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                                      >  1:687189/1‑59 (MQ=255)
             ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                              >  1:1235824/1‑66 (MQ=255)
             ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                              >  1:1526610/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)
                gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                                                       <  1:569964/70‑1 (MQ=255)
                gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                                                       <  1:2597733/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:4451249/1‑40 (MQ=255)
                               ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                                                      >  1:1414845/1‑40 (MQ=255)
                               ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                                          <  1:5615593/52‑1 (MQ=255)
                                tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                              >  1:1723349/1‑47 (MQ=255)
                                tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                              >  1:4719233/1‑47 (MQ=255)
                                tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                              >  1:1707352/1‑47 (MQ=255)
                                  cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACg                                    <  1:1722190/71‑1 (MQ=255)
                                    tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGt                                   <  1:3375919/70‑1 (MQ=255)
                                    tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGt                                   <  1:3382378/70‑1 (MQ=255)
                                    tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                                               >  1:2239762/1‑58 (MQ=255)
                                    tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                                               >  1:5256285/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)
                                                       ttGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGaa                            <  1:2171945/57‑1 (MQ=255)
                                                       ttGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGAc                   >  1:1859496/2‑67 (MQ=255)
                                                       cTGCGTCTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACt                  <  1:1194552/68‑1 (MQ=255)
                                                       cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                                               >  1:3899584/1‑39 (MQ=255)
                                                       cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                                            <  1:2842309/42‑1 (MQ=255)
                                                       cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                                            <  1:3609614/42‑1 (MQ=255)
                                                         gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa              >  1:3097322/1‑70 (MQ=255)
                                                         gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa              >  1:180735/1‑70 (MQ=255)
                                                         gCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGa                                      >  1:3131241/1‑46 (MQ=255)
                                                         gCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAAt                           >  1:3352922/1‑57 (MQ=255)
                                                         gCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAAt                           >  1:5092699/1‑57 (MQ=255)
                                                         gCGACATGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGACTGGTa              >  1:2420630/1‑70 (MQ=255)
                                                          cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTg                 >  1:2286052/1‑66 (MQ=255)
                                                          cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc             <  1:1238824/70‑1 (MQ=255)
                                                          cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc             <  1:3566441/70‑1 (MQ=255)
                                                          cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACa            <  1:4997379/71‑1 (MQ=255)
                                                          cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACa            <  1:945171/71‑1 (MQ=255)
                                                           gACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACaa           <  1:3450789/71‑1 (MQ=255)
                                                            aCTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAAc          >  1:3771725/1‑71 (MQ=255)
                                                              ttGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTc        >  1:418955/1‑71 (MQ=255)
                                                                 cctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAg     >  1:2857378/1‑71 (MQ=255)
                                                                  ctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTa                                  >  1:4903174/1‑41 (MQ=255)
                                                                  ctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTa                                  >  1:4872348/1‑41 (MQ=255)
                                                                     ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTa  >  1:5526511/1‑70 (MQ=255)
                                                                     ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTa  >  1:3282601/1‑70 (MQ=255)
                                                                     ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGATa  >  1:2386521/1‑70 (MQ=255)
                                                                      |                                                                    
GTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTA  >  Exported/832150‑832288

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: