Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,612 T→C 22.2% Y444H (TAC→CAC)  narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE836,6120TC22.2% 151.0 / 42.0 81Y444H (TAC→CAC) narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (35/28);  new base C (11/7);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

TGACGGTAAAGTTGATACCCGTGCGCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCA  >  minE/836545‑836680
                                                                   |                                                                    
tGACGGTAAAGTTGATACCCGTGCTCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTAcc                                                                   >  1:2997884/1‑71 (MQ=255)
  aCGGTAAAGTTGATACCCGTGCGCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGCACCGt                                                                 <  1:863890/71‑1 (MQ=255)
   cGGTAAAGTTGATACCCGTGCGCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGt                                                                 >  1:1878755/1‑70 (MQ=255)
       aaaGTTGATACCCGTGCGCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCt                                                            >  1:1780408/1‑71 (MQ=255)
        aaGTTGATACCCGTGCGCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCt                                                            <  1:3376845/70‑1 (MQ=255)
        aaGTTGATACCCGTGCGCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTg                                                           <  1:1157408/71‑1 (MQ=255)
             gATACCCGTGCGCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGa                                                       <  1:961623/70‑1 (MQ=255)
              aTACCCGTGCGCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGAtt                                                     <  1:1500196/71‑1 (MQ=255)
               tACCCGTGCGCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTg                                                    >  1:3216574/1‑71 (MQ=255)
                aCCCGTGCGCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGc                                                   >  1:1367429/1‑71 (MQ=255)
                    gTGCGCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGCACCGTTATCTGGCGATTGCTAAc                                               >  1:2233303/1‑71 (MQ=255)
                        gcTGGAAGAGGTCGGTCTGACCGAAGCCCAGGGACAGGAGATGTACCGTTATCt                                                            <  1:1776727/54‑1 (MQ=255)
                        gcTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCt                                                            <  1:418949/54‑1 (MQ=255)
                        gcTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGc                                                         <  1:2145397/57‑1 (MQ=255)
                         cTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACg                                           <  1:1261099/70‑1 (MQ=255)
                         cTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACg                                           <  1:334060/70‑1 (MQ=255)
                          tGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGtt                                                                <  1:930902/48‑1 (MQ=255)
                           ggAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACg                                           >  1:1206168/1‑68 (MQ=255)
                            gAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACg                                           >  1:162032/1‑67 (MQ=255)
                            gAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGCACCGTTATCTGGCGATTGCTAACTACGAAGa                                       <  1:3374650/71‑1 (MQ=255)
                             aaGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGCACCGTTATCTGGCGATTGCTAACTACg                                           >  1:2574324/1‑66 (MQ=255)
                              agagGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACg                                           >  1:497729/1‑65 (MQ=255)
                               gagGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGt                                                                      >  1:3497204/1‑37 (MQ=255)
                                agGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGCACCGTTATCTg                                                           >  1:1629306/1‑47 (MQ=255)
                                agGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGCACCGTTATCTg                                                           >  1:3188181/1‑47 (MQ=255)
                                agGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGCACCGTTATCTg                                                           >  1:2300576/1‑47 (MQ=255)
                                   tcggtcTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCttt                                <  1:1235938/71‑1 (MQ=255)
                                    cggtcTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGa                                                       >  1:1645453/1‑47 (MQ=255)
                                    cggtcTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGa                                                       >  1:2199123/1‑47 (MQ=255)
                                    cggtcTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCtt                                 >  1:47269/1‑69 (MQ=255)
                                    cggtcTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTtg                               >  1:3094574/1‑71 (MQ=255)
                                      gtcTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTtgtg                             <  1:824376/71‑1 (MQ=255)
                                      gtcTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTtgtg                             <  1:1467970/71‑1 (MQ=255)
                                      gtcTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTtgtg                             <  1:449511/71‑1 (MQ=255)
                                      gtcTGACCGAAGCCCAGGCACAGGAGATGCACCGTTATCTGGCGATTGCTAACTACGaa                                         <  1:3922392/59‑1 (MQ=255)
                                       tcTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTgtgg                            >  1:3314841/1‑71 (MQ=255)
                                       tcTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTgtgg                            >  1:2564295/1‑71 (MQ=255)
                                        cTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTgtggt                           >  1:152027/1‑71 (MQ=255)
                                         tGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTg                                                    >  1:2235632/1‑45 (MQ=255)
                                         tGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGaa                                         <  1:3489769/56‑1 (MQ=255)
                                         tGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTgtggtg                          >  1:558584/1‑71 (MQ=255)
                                           aCCGAAGCCCAGGCACAGGAGATGTACCGTTATCTg                                                           >  1:2699011/1‑36 (MQ=255)
                                           aCCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGAt                                                      <  1:1538584/41‑1 (MQ=255)
                                           aCCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGAt                                                      <  1:3538876/41‑1 (MQ=255)
                                           aCCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTaa                                                >  1:584162/1‑47 (MQ=255)
                                           aCCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGAt                                      >  1:1432204/1‑57 (MQ=255)
                                            ccGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGcc                        >  1:958585/1‑70 (MQ=255)
                                            ccGAAGCCCAGGCACAGGAGATGCACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGcc                        >  1:1199309/1‑70 (MQ=255)
                                             cGAAGCCCAGGCACAGGAGATGCACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTtg                               <  1:2520652/62‑1 (MQ=255)
                                             cGAAGCCCAGGCACAGGAGATGCACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTtg                               <  1:2380132/62‑1 (MQ=255)
                                              gAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTaa                                                >  1:2136182/1‑44 (MQ=255)
                                              gAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGag                     <  1:3142984/71‑1 (MQ=255)
                                               aaGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGag                     >  1:1893541/1‑70 (MQ=255)
                                               aaGCCCAGGCACAGGAGATGCACCGTTATCTGGCGATTg                                                    >  1:98782/1‑39 (MQ=255)
                                               aaGCCCAGGCACAGGAGATGCACCGTTATCTGGCGATTGCTAACTACGaa                                         <  1:515421/50‑1 (MQ=255)
                                               aaGCCCAGGCACAGGAGATGCACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGag                     >  1:1270479/1‑70 (MQ=255)
                                                aGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGag                     >  1:265921/1‑69 (MQ=255)
                                                 gCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATc                                     >  1:1218274/1‑52 (MQ=255)
                                                 gCCCAGGCACAGGAGATGCACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGagtag                  >  1:174001/1‑71 (MQ=255)
                                                    cAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTgtggtg                          <  1:2709056/60‑1 (MQ=255)
                                                    cAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGagt                    >  1:1714251/1‑66 (MQ=255)
                                                     aGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCAt              >  1:1651741/1‑71 (MQ=255)
                                                     aGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCAt              >  1:638637/1‑71 (MQ=255)
                                                     aGGCACAGGAGATGCACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCAt              >  1:182283/1‑71 (MQ=255)
                                                       gCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTgtgg                            >  1:505876/1‑55 (MQ=255)
                                                       gCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATc             <  1:588006/70‑1 (MQ=255)
                                                       gCACAGGAGATGCACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCg            <  1:156424/71‑1 (MQ=255)
                                                           aGGAGATTTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCttt                                >  1:902816/1‑47 (MQ=255)
                                                           aGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGaa        >  1:214990/1‑71 (MQ=255)
                                                           aGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGa         <  1:568056/70‑1 (MQ=255)
                                                             gagaTGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAAc       >  1:1729829/1‑70 (MQ=255)
                                                             gagaTGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACt      <  1:677847/71‑1 (MQ=255)
                                                             gagaTGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACt      >  1:1279638/1‑71 (MQ=255)
                                                             gagaTGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACt      >  1:1160337/1‑71 (MQ=255)
                                                             gagaTGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCGTTGTGGTGCCGAGTAGTCATCGTGAACt      <  1:568678/71‑1 (MQ=255)
                                                             gagaTGTACAGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAAc       >  1:2268586/1‑70 (MQ=255)
                                                              agaTGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTg     <  1:803686/71‑1 (MQ=255)
                                                              agaTGCACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTg     >  1:2790088/1‑71 (MQ=255)
                                                               gaTGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGcc                        <  1:270943/51‑1 (MQ=255)
                                                                 tGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCa               <  1:57929/58‑1 (MQ=255)
                                                                 tGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGACCTGGCa  <  1:511603/71‑1 (MQ=255)
                                                                 tGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCa  <  1:3826428/71‑1 (MQ=255)
                                                                   tACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGa         <  1:1774241/62‑1 (MQ=255)
                                                                   |                                                                    
TGACGGTAAAGTTGATACCCGTGCGCTGGAAGAGGTCGGTCTGACCGAAGCCCAGGCACAGGAGATGTACCGTTATCTGGCGATTGCTAACTACGAAGATCGCTTTGTGGTGCCGAGTAGTCATCGTGAACTGGCA  >  minE/836545‑836680

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: