Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,740 G→A 11.8% L71L (CTG→CTA) ‡ narI → nitrate reductase 1, gamma (cytochrome b(NR)) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE837,7400GA11.8% 150.3 / 9.1 59L71L (CTG→CTA) ‡narInitrate reductase 1, gamma (cytochrome b(NR)) subunit
Reads supporting (aligned to +/- strand):  ref base G (28/24);  new base A (2/5);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGT  >  minE/837671‑837809
                                                                     |                                                                     
gAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa                                                                      >  1:4804580/1‑71 (MQ=255)
    cTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgcc                                                                  <  1:1498757/71‑1 (MQ=255)
     tGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg                                                                 >  1:1827984/1‑71 (MQ=255)
          tCGAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTg                                                            <  1:5667574/71‑1 (MQ=255)
            gAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGa                                                          >  1:1772263/1‑71 (MQ=255)
             aatcTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGAt                                                         >  1:4149972/4‑71 (MQ=255)
              atcTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACt                                                             <  1:4370513/64‑1 (MQ=255)
              atcTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATg                                                        <  1:3530681/69‑1 (MQ=255)
              atcTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATg                                                        <  1:4836129/69‑1 (MQ=255)
                 tGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTAt                                                     >  1:3413289/1‑71 (MQ=255)
                 tGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTAt                                                     >  1:3813579/1‑71 (MQ=255)
                  gTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGt                                                       <  1:1693949/68‑1 (MQ=255)
                       atatCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACCCCGCACTGGATg                                                        <  1:1901880/62‑1 (MQ=255)
                        tatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGcc                                               <  1:2106917/70‑1 (MQ=255)
                           cGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGc                                           >  1:4708655/1‑71 (MQ=255)
                               aTTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGt                                                       <  1:3749098/55‑1 (MQ=255)
                               aTTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGcc                                       <  1:3701736/71‑1 (MQ=255)
                                ttCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGcc                                       <  1:4949025/70‑1 (MQ=255)
                                  cTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                            >  1:2619851/1‑47 (MQ=255)
                                  cTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAt                                    >  1:365107/1‑71 (MQ=255)
                                  cTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAt                                    >  1:974300/1‑71 (MQ=255)
                                  cTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGACGCACTGGATGTATGAAGCCTGGCTGCCGAt                                    >  1:5258217/1‑71 (MQ=255)
                                     gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTg                                  <  1:878421/70‑1 (MQ=255)
                                     gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGa                                 <  1:3201461/71‑1 (MQ=255)
                                      ggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTAt                                                     <  1:15432/50‑1 (MQ=255)
                                       gATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAg                               >  1:3198810/1‑71 (MQ=255)
                                        aTTTTTGTCGGTCACTTCTTCGGTATGCTGAcgcc                                                                  >  1:3579565/1‑35 (MQ=255)
                                        aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                            >  1:4697662/1‑41 (MQ=255)
                                        aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                            >  1:3820107/1‑41 (MQ=255)
                                        aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                              <  1:488473/71‑1 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                            <  1:2713587/39‑1 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                            <  1:4893634/39‑1 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                            <  1:4550046/39‑1 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAg                               >  1:174615/1‑68 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                              >  1:3932177/1‑69 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTg                             >  1:24118/1‑70 (MQ=255)
                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTg                             <  1:3362383/70‑1 (MQ=255)
                                                cGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGa                                     >  1:5669370/1‑56 (MQ=255)
                                                cGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGa                                     >  1:1508084/1‑56 (MQ=255)
                                                cGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGa                                     >  1:3958079/1‑56 (MQ=255)
                                                 ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTg                                             >  1:3190265/1‑47 (MQ=255)
                                                 ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGa                      >  1:2060130/1‑70 (MQ=255)
                                                 ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGa                      >  1:3720324/1‑70 (MQ=255)
                                                 ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGa                      >  1:5684658/1‑70 (MQ=255)
                                                   tCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACATaaaa                   <  1:2201898/71‑1 (MQ=255)
                                                     aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAtt                                   <  1:5530285/53‑1 (MQ=255)
                                                       ttcttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaaa                   <  1:3202934/67‑1 (MQ=255)
                                                         cttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCa              <  1:868563/70‑1 (MQ=255)
                                                            cGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAg                               >  1:2293735/1‑50 (MQ=255)
                                                            cGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCaa             <  1:5480955/68‑1 (MQ=255)
                                                            cGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGt          >  1:1658256/1‑71 (MQ=255)
                                                                 tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa                          <  1:3684330/50‑1 (MQ=255)
                                                                 tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa                          <  1:5179832/50‑1 (MQ=255)
                                                                 tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGc      <  1:4568202/70‑1 (MQ=255)
                                                                   cTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTg       >  1:1058779/1‑67 (MQ=255)
                                                                   cTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTg       >  1:5624736/1‑67 (MQ=255)
                                                                   cTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtg    <  1:4619851/70‑1 (MQ=255)
                                                                   cTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtg    <  1:319788/70‑1 (MQ=255)
                                                                    tGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtggt  >  1:1184410/1‑71 (MQ=255)
                                                                     |                                                                     
GAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGT  >  minE/837671‑837809

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