Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,545,084 C→T 16.0% G231R (GGG→AGG)  narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,545,0840CT16.0% 126.3 / 15.2 50G231R (GGG→AGG) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (27/15);  new base T (4/4);  total (31/19)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCATTCGCCAGCGACAT  >  W3110S.gb/1545017‑1545153
                                                                   |                                                                     
tgaaGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCt                                                                      >  1:2185358/4‑69 (MQ=255)
tGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGa                                                                    >  1:2878616/1‑71 (MQ=255)
tGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTAGGCGCGTGAACTGGCGATATCATTCATCCCTGAc                                                                   >  1:255660/1‑71 (MQ=255)
          gTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATgcggc                                                           >  1:155390/1‑70 (MQ=255)
           tAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATgcggc                                                           >  1:3633215/1‑69 (MQ=255)
           tAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGa                                                         >  1:2210766/1‑71 (MQ=255)
                   ggAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCa                                                 <  1:77099/71‑1 (MQ=255)
                   ggAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCa                                                 <  1:3681899/71‑1 (MQ=255)
                   ggAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCa                                                 <  1:1465587/71‑1 (MQ=255)
                   ggAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGAGATCGTGGCa                                                 <  1:3090998/71‑1 (MQ=255)
                            cGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAg                                         >  1:1487071/1‑70 (MQ=255)
                             ggCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCTTGACCATGCGGCGATCGTGGCAATCGCCAg                                         >  1:3926261/1‑69 (MQ=255)
                             ggCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCTTGACCATGCGGCGATCGTGGCAATCGCCAg                                         >  1:1117059/1‑69 (MQ=255)
                             ggCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCTTGACCATGCGGCGATCGTGGCAATCGCCAg                                         >  1:2215110/1‑69 (MQ=255)
                                 aTTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGa                                                         <  1:3327757/49‑1 (MQ=255)
                                 aTTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGa                                                         <  1:2883444/49‑1 (MQ=255)
                                 aTTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCgg                                   <  1:2597682/71‑1 (MQ=255)
                                  ttGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGcc                                           <  1:630505/62‑1 (MQ=255)
                                  ttGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGcc                                           <  1:767901/62‑1 (MQ=255)
                                  ttGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGcc                                           <  1:2678121/62‑1 (MQ=255)
                                   tGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGt                                                     >  1:3614366/1‑51 (MQ=255)
                                    gAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCa                                                                 >  1:2158103/1‑38 (MQ=255)
                                    gAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCa                                                                 >  1:3140064/1‑38 (MQ=255)
                                    gAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCa                                                                 >  1:837029/1‑38 (MQ=255)
                                    gAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTAc                                >  1:3267169/1‑71 (MQ=255)
                                    gAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTAc                                >  1:1660580/1‑71 (MQ=255)
                                        cgcgTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCaa                            <  1:2442874/71‑1 (MQ=255)
                                        cgcgTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCaa                            <  1:740059/71‑1 (MQ=255)
                                           gTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAAt                          <  1:865845/70‑1 (MQ=255)
                                           gTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATc                         <  1:2038550/71‑1 (MQ=255)
                                                   gCGATATCATCCCTCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCgg                  <  1:1650403/70‑1 (MQ=255)
                                                       tatCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCAt                      >  1:3092203/1‑62 (MQ=255)
                                                        atCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACGCAAATCCa                       >  1:1638543/1‑60 (MQ=255)
                                                        atCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCa                       >  1:1169573/1‑60 (MQ=255)
                                                        atCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCa                       >  1:1235901/1‑60 (MQ=255)
                                                        atCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCa                       >  1:3719151/1‑60 (MQ=255)
                                                        atCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCATTCg            >  1:1056493/1‑71 (MQ=255)
                                                        atCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCATTCg            >  1:2206991/1‑71 (MQ=255)
                                                        atCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCATTCg            >  1:1232582/1‑71 (MQ=255)
                                                         tcattcatCCTTGACCATGCGGCGATCTTGGCAATCGCCAGTAGc                                     <  1:1957249/45‑1 (MQ=38)
                                                         tcattcatCCTTGACCATGCGGCGATCGTGGCAATCGCCAGTAGc                                     <  1:1254986/45‑1 (MQ=255)
                                                         tcattcatCCTTGACCATGCGGCGATCGTGGCAATCGCCAGTAGc                                     <  1:1491675/45‑1 (MQ=255)
                                                         tcattcatCCTTGACCATGCGGCGATCGTGGCAATCGCCAGTAGc                                     <  1:189136/45‑1 (MQ=255)
                                                         tcattcatCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTa                                 >  1:1054380/1‑49 (MQ=255)
                                                         tcattcatCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCa                >  1:3431927/1‑66 (MQ=255)
                                                         tcattcatCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCa                >  1:2867111/1‑66 (MQ=255)
                                                         tcattcatCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCa                >  1:2704696/1‑66 (MQ=255)
                                                         tcattcatCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCAtt              >  1:3495266/1‑68 (MQ=255)
                                                         tcattcatCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCAtt              >  1:2052922/1‑68 (MQ=255)
                                                           attcatCCTTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCgg                                   >  1:1961692/1‑45 (MQ=255)
                                                                  ccTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCATTCGCCAGCGACAt  >  1:2607149/1‑71 (MQ=255)
                                                                   |                                                                     
TGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGCGGCGATCGTGGCAATCGCCAGTAGCGGTACCCAAATCCATGCGGCATTCGCCAGCGACAT  >  W3110S.gb/1545017‑1545153

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: