Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,620 T→A 6.5% G1026G (GGT→GGA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,6200TA 6.5% 166.6 / 2.7 61G1026G (GGT→GGAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (22/35);  new base A (2/2);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACC  >  minE/834550‑834690
                                                                      |                                                                      
aGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGt                                                                        <  1:2348408/71‑1 (MQ=255)
 gCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGt                                                                        >  1:2913607/1‑70 (MQ=255)
  ccGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATa                                                                      <  1:272517/71‑1 (MQ=255)
  ccGCGCAAAATTATCTCCAGCCCGACATGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTa                                                                      <  1:2163643/71‑1 (MQ=255)
           aTTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGt                                                             >  1:1232942/1‑71 (MQ=255)
           aTTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGt                                                             <  1:2594773/71‑1 (MQ=255)
            ttATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGtt                                                            >  1:731391/1‑71 (MQ=255)
                 tcCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTa                                                                      <  1:2579058/56‑1 (MQ=255)
                     gCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCa                                                                 <  1:593070/57‑1 (MQ=255)
                       ccGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTa                                                                      >  1:2796898/1‑50 (MQ=255)
                       ccGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGAt                                                 <  1:1347082/71‑1 (MQ=255)
                           ccTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCTGTTa                                                                      >  1:3585625/1‑46 (MQ=255)
                           ccTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTa                                                                      >  1:2234482/1‑46 (MQ=255)
                            ctggtctggtCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCAt                                            >  1:307503/1‑71 (MQ=255)
                             tggtctggtcTGGAAGATGAACACGTTTCTTACAACGCCGGTTa                                                                      <  1:1593722/44‑1 (MQ=255)
                              ggtctggtctGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTc                                                           <  1:1974179/54‑1 (MQ=255)
                               gtctggtctgGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                                             >  1:2797452/1‑67 (MQ=255)
                                tctggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCAt                                            >  1:3385649/1‑67 (MQ=255)
                                tctggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                                        >  1:4145513/1‑71 (MQ=255)
                                tctggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                                        >  1:2938135/1‑71 (MQ=255)
                                tctggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                                        >  1:2624022/1‑71 (MQ=255)
                                tctggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                                        <  1:2396592/71‑1 (MQ=255)
                                   ggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATcc                                               <  1:2760870/61‑1 (MQ=255)
                                   ggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATcc                                               >  1:2496529/1‑61 (MQ=255)
                                     tctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGtt                                                            <  1:1017930/46‑1 (MQ=255)
                                     tctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGtt                                                            <  1:186904/46‑1 (MQ=255)
                                       tggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCa                                                          >  1:189373/1‑46 (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:2564465/71‑1 (MQ=255)
                                            gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTgg                            <  1:147719/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:1312727/1‑71 (MQ=255)
                                                aaCACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtcgt                        <  1:3338650/71‑1 (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:441215/71‑1 (MQ=255)
                                                                     gTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAc   <  1:454433/71‑1 (MQ=255)
                                                                      ttACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAcc  <  1:2115458/71‑1 (MQ=255)
                                                                      ttACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATAAGGATCAcc  <  1:2147240/71‑1 (MQ=255)
                                                                      |                                                                      
AGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACC  >  minE/834550‑834690

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: