Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,281,601 Δ1 bp 7.7% coding (1738/3744 nt) narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,281,6010G. 7.7% 169.8 / 7.7 52coding (1738/3744 nt)narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (29/19);  new base . (2/2);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTG  >  NC_000913/1281461‑1281750
                                                                                                                                            |                                                                                                                                                     
cGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGAGTCAGGaaaaa                                                                                                                                            >  1:770641/1‑150 (MQ=255)
cGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGAGTCAGGaaaaa                                                                                                                                            >  1:239179/1‑150 (MQ=255)
cGCACGGGCGCTCGATGATTATCGTCGGGGCGGGGCGGAACCACTGGTCTCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa                                                                                                                                              <  1:228981/150‑1 (MQ=255)
  cACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaaaa                                                                                                                                            >  1:72741/1‑150 (MQ=255)
    cGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCTGGATAAACt                                                                                                                                          >  1:820677/1‑150 (MQ=255)
     ggTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTg                                                                                                                                         <  1:412591/150‑1 (MQ=255)
        cGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGt                                                                                                                                      >  1:712357/1‑150 (MQ=255)
         gTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTc                                                                                                                                     <  2:141428/150‑1 (MQ=255)
            cGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCTGGAAAAACTGCGTCCGc                                                                                                                                  >  2:805288/1‑150 (MQ=255)
             gatgatTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCa                                                                                                                                 >  1:23015/1‑150 (MQ=255)
                  ttATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCTGAGCGGGGGCGGCTGGGCGCACTATGAGTCAGGAAAAACTGCGTCCGCAAACCGGc                                                                                                                          <  2:239179/150‑1 (MQ=255)
                  ttATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGAGTCAGGAAAAACTGCGTCCGCAAACCGGc                                                                                                                          <  2:770641/150‑1 (MQ=255)
                   tATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCgg                                                                                                                           >  1:583888/1‑150 (MQ=255)
                     tcgtcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCt                                                                                                                         <  2:72741/150‑1 (MQ=255)
                      cgtcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTg                                                                                                                        >  1:1139916/1‑150 (MQ=255)
                       gtcgGGGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGAGCAACATGCTGATTTTCTGCGGCGGTGTCGGCCAGAGCGGGGGCGGCTGGGCGCACAATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTgg                                                                                                                       <  1:447126/150‑1 (MQ=255)
                          ggTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAg                                                                                                                    >  1:472714/1‑150 (MQ=255)
                              cGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGc                                                                                                                >  1:486368/1‑150 (MQ=255)
                              cGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGc                                                                                                                >  2:1100765/1‑150 (MQ=255)
                                   cTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCCGCCGCTGGCg                                                                                                           >  2:807337/1‑150 (MQ=255)
                                   cTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTggag                                                                                                           >  1:229894/1‑148 (MQ=255)
                                            tGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCtt                                                                                                  <  2:151661/150‑1 (MQ=255)
                                                  cACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTgg                                                                                            >  2:1042318/1‑150 (MQ=255)
                                                  cACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTgg                                                                                            >  2:691488/1‑150 (MQ=255)
                                                        gATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACGATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGt                                                                                      >  1:1056890/1‑150 (MQ=255)
                                                             gAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGgc                                                                                 >  2:116288/1‑150 (MQ=255)
                                                                cTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGcgcg                                                                              >  2:577786/1‑150 (MQ=255)
                                                                   tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGGCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTca                                                                           <  1:1286255/150‑1 (MQ=255)
                                                                      tGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGTAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTTACTGGCAGCGTCAGGCGCGTCACAt                                                                        >  1:285900/1‑150 (MQ=255)
                                                                      tGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACAt                                                                        <  2:740016/150‑1 (MQ=255)
                                                                       ggTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATg                                                                       <  2:774401/150‑1 (MQ=255)
                                                                            gATCATCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAATAACTGCGTCCGCAATCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGAGTGGCAGCGTCCGGCGCGTCCCGTGagctg                                                                  >  1:178219/1‑146 (MQ=255)
                                                                                         aTTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTAttta                                                     >  2:1196953/1‑149 (MQ=255)
                                                                                           tgactccGGCTGTGTCGGTCAGAGCGGGGGCGGCTGTGCGCACTATGGAGGTCAGGAAAAACTGCGGCCGCATACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGGCCGGCGCGTCACATGAACAGCACTTCTTATTTCta                                                   <  2:892604/144‑1 (MQ=255)
                                                                                               tGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAAc                                               >  2:37332/1‑150 (MQ=255)
                                                                                                gCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGGCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAAcc                                              <  1:151032/150‑1 (MQ=255)
                                                                                                  ggctggctCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCAc                                            <  1:675913/150‑1 (MQ=255)
                                                                                                   gCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACt                                           >  2:250401/1‑150 (MQ=255)
                                                                                                       tgtCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCCCTTCTTATTTCTATAACCACTccag                                       >  1:367365/1‑150 (MQ=255)
                                                                                                        gtCGGTCAGCGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACGTCTGATTTCTATACCCACTccagc                                      >  2:766871/1‑150 (MQ=255)
                                                                                                         tCGGTCAGAGCGGGGGCGGCTGGGCGCACGATGTAGGGCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagcc                                     <  2:198494/150‑1 (MQ=255)
                                                                                                         gcGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagcc                                     <  1:110103/149‑1 (MQ=255)
                                                                                                           ggTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagccag                                   >  2:467825/1‑150 (MQ=255)
                                                                                                             tCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTg                                 >  2:172997/1‑150 (MQ=255)
                                                                                                                 tgCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAATAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGt                             <  2:56626/149‑1 (MQ=255)
                                                                                                                     ggggCGGCTGGGCGCACTAAGTAGCTCAGGAAAAACTGCGCCCGCAACCCGGCGGGCAGCCGCGGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATg                         <  2:697043/150‑1 (MQ=255)
                                                                                                                                   cACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggcg           >  1:781096/1‑150 (MQ=255)
                                                                                                                                   cACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggcg           >  1:1295472/1‑150 (MQ=255)
                                                                                                                                     cTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGAGTAAGCCCATGACAGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATCACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggcggc         <  1:828667/150‑2 (MQ=255)
                                                                                                                                         gTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAagag     <  2:44650/150‑1 (MQ=255)
                                                                                                                                          tAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGATACCGTCACGGCGGAAGAGt    >  1:590891/1‑150 (MQ=255)
                                                                                                                                           aGGTCCGGACAAACTGCGTCCGCAAACCGGCAGGCAGCCGCTGGCGTTAGCCCTTGACTGGCTGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGtt   <  1:856429/150‑1 (MQ=255)
                                                                                                                                            ggTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCCGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTg  >  1:649695/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
CGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTG  >  NC_000913/1281461‑1281750

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