Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,644 C→T 8.4% I1034I (ATC→ATT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,6440CT 8.4% 218.9 / 8.1 71I1034I (ATC→ATTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (28/37);  new base T (4/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

CTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGA  >  minE/834578‑834712
                                                                  |                                                                    
ctggtctggtCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCAt                                                                  >  1:307503/1‑71 (MQ=255)
   gtctggtctgGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                                                                   >  1:2797452/1‑67 (MQ=255)
    tctggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCAt                                                                  >  1:3385649/1‑67 (MQ=255)
    tctggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                                                              >  1:2624022/1‑71 (MQ=255)
    tctggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                                                              <  1:2396592/71‑1 (MQ=255)
    tctggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                                                              >  1:2938135/1‑71 (MQ=255)
    tctggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                                                              >  1:4145513/1‑71 (MQ=255)
       ggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATcc                                                                     >  1:2496529/1‑61 (MQ=255)
       ggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATcc                                                                     <  1:2760870/61‑1 (MQ=255)
             gAAGATTAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctctc                                                     >  1:2752234/1‑71 (MQ=255)
             gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATcc                                                                     <  1:2804917/55‑1 (MQ=255)
             gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctct                                                      <  1:1371975/70‑1 (MQ=255)
             gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGct                                                        <  1:233608/68‑1 (MQ=255)
               aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATgg                                                                <  1:2570987/58‑1 (MQ=255)
               aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATgg                                                                <  1:735766/58‑1 (MQ=255)
               aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                                                   <  1:981661/71‑1 (MQ=255)
                gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATccc                                                                    <  1:3406873/53‑1 (MQ=255)
                gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTgg                                                  <  1:147719/71‑1 (MQ=255)
                gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTgg                                                  <  1:2564465/71‑1 (MQ=255)
                gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTgg                                                  <  1:2816892/71‑1 (MQ=255)
                gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTgg                                                  <  1:3062528/71‑1 (MQ=255)
                  tGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATg                                                                 <  1:1138669/54‑1 (MQ=255)
                  tGAACACGTTTCTTACAACGCCGGATACACCAACGTTCACGAGCTGATCCCATg                                                                 <  1:3252305/54‑1 (MQ=255)
                  tGAACACGTTTCTTACAACCCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtc                                                >  1:3910468/1‑71 (MQ=255)
                   gAACACGTTTCTTACAACGCCGGTTACACCCACGTTCACGAGCTGATCCCAt                                                                  >  1:3620640/1‑52 (MQ=255)
                    aaCACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtcgt                                              <  1:3338650/71‑1 (MQ=255)
                    aaCACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtcgt                                              >  1:1312727/1‑71 (MQ=255)
                            tCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACt                                      >  1:1033494/1‑71 (MQ=255)
                            tCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACt                                      >  1:1602219/1‑71 (MQ=255)
                               tACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGct                                                        <  1:1552652/50‑1 (MQ=255)
                               tACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTgg                                                  <  1:37391/56‑1 (MQ=255)
                                aCAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctctc                                                     >  1:3837979/1‑52 (MQ=255)
                                aCAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGt                                    <  1:3511710/69‑1 (MQ=255)
                                aCAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGt                                    <  1:2808065/69‑1 (MQ=255)
                                aCAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTa                                   <  1:2108158/70‑1 (MQ=255)
                                   aCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCAg                               >  1:430498/1‑71 (MQ=255)
                                   aCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCAg                               >  1:2179343/1‑71 (MQ=255)
                                   aCGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                        >  1:3689754/1‑62 (MQ=255)
                                   aCGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                        >  1:172358/1‑62 (MQ=255)
                                   aCGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAg                               <  1:3209890/71‑1 (MQ=255)
                                       cGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                        <  1:1037343/58‑1 (MQ=255)
                                       cGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATc                           <  1:226872/71‑1 (MQ=255)
                                         gTTACACCAACGTTCACGAGCTGATCCCATGGCGTACg                                                          <  1:833036/38‑1 (MQ=255)
                                         gTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAc                         <  1:454433/71‑1 (MQ=255)
                                         gTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAc                         <  1:441215/71‑1 (MQ=255)
                                          ttACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAcc                        <  1:2115458/71‑1 (MQ=255)
                                          ttACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATAAGGATCAcc                        <  1:2147240/71‑1 (MQ=255)
                                          atACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCAGGATCatc                        >  1:2216851/2‑69 (MQ=255)
                                           tACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCa                                            >  1:2353341/1‑50 (MQ=255)
                                           tACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCa                                            >  1:27537/1‑50 (MQ=255)
                                           tACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCa                          >  1:631924/1‑68 (MQ=255)
                                           tACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAc                         >  1:4120813/1‑69 (MQ=255)
                                           tACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAc                         >  1:786533/1‑69 (MQ=255)
                                            acacCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGa                             <  1:3736207/64‑1 (MQ=255)
                                            acacCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAg                      >  1:670615/1‑71 (MQ=255)
                                            acacCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAg                      >  1:272026/1‑71 (MQ=255)
                                             cacCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAgg                              <  1:504026/62‑1 (MQ=255)
                                              acCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGGAACTGTATCAGGATCACCAGTg                    <  1:1025442/71‑1 (MQ=255)
                                              acCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCa                                         <  1:3904627/50‑1 (MQ=255)
                                              acCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTg                    <  1:3576880/71‑1 (MQ=255)
                                              acCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTg                    >  1:3068521/1‑71 (MQ=255)
                                              acCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTg                    <  1:3505227/71‑1 (MQ=255)
                                               ccAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTg                    <  1:3910876/70‑1 (MQ=255)
                                                 aaCGTTGACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGa                  >  1:538444/1‑70 (MQ=255)
                                                 aaCGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCAGGATCATCCGTGGAt                 >  1:1990523/1‑71 (MQ=255)
                                                   cGTTCACGAGCTGATCCCATGGGGTACGCTCTCTGGTCGGCAGc                                          <  1:3624101/44‑1 (MQ=255)
                                                     ttCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCAGGATCATCCGTg                    <  1:1730817/64‑1 (MQ=255)
                                                     ttCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                        >  1:314352/1‑44 (MQ=255)
                                                     ttCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                        >  1:2833849/1‑44 (MQ=255)
                                                      tCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCAGGATCATCCGTg                    <  1:3558390/63‑1 (MQ=255)
                                                      tCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTa                                   >  1:2609143/1‑48 (MQ=255)
                                                          gAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCa                       <  1:1271156/56‑1 (MQ=255)
                                                           aGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTc       <  1:2607396/71‑1 (MQ=255)
                                                                aTCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGa  <  1:2740569/71‑1 (MQ=255)
                                                                 tCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGa  >  1:664217/1‑70 (MQ=255)
                                                                  cccATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGt             >  1:1653097/1‑58 (MQ=255)
                                                                  cccATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGt             >  1:2933068/1‑58 (MQ=255)
                                                                  |                                                                    
CTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGA  >  minE/834578‑834712

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: