Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,705 C→T 26.5% L60L (CTG→TTG)  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,7050CT26.5% 113.1 / 40.1 68L60L (CTG→TTG) narInitrate reductase 1, gamma (cytochrome b(NR)) subunit
Reads supporting (aligned to +/- strand):  ref base C (20/30);  new base T (3/15);  total (23/45)
Fisher's exact test for biased strand distribution p-value = 8.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAT  >  minE/837638‑837775
                                                                   |                                                                      
ggcgTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTg                                                                      <  1:5327366/70‑1 (MQ=255)
     ccagccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTc                                                                        >  1:4871605/1‑63 (MQ=255)
      cagccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttt                                                               >  1:3464297/1‑71 (MQ=255)
      cagccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttt                                                               >  1:3620679/1‑71 (MQ=255)
      cagccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttt                                                               <  1:5195918/71‑1 (MQ=255)
      cagccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttt                                                               >  1:2337270/1‑71 (MQ=255)
      cagccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAtt                                                                <  1:3364014/70‑1 (MQ=255)
             aTGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCgg                                                        <  1:20724/71‑1 (MQ=255)
              tGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttttt                                                             >  1:2796961/1‑65 (MQ=255)
               gCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGt                                                       >  1:3104546/1‑70 (MQ=255)
                cTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTggg                                                                    <  1:1202644/56‑1 (MQ=255)
                  ggATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCt                                                                       >  1:5453522/1‑51 (MQ=255)
                  ggATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCt                                                                       >  1:4506375/1‑51 (MQ=255)
                  ggATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCt                                                                       >  1:1940993/1‑51 (MQ=255)
                   gATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCActt                                                  <  1:537561/71‑1 (MQ=255)
                   gATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCAc                                                    <  1:753902/69‑1 (MQ=255)
                    aTCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttttt                                                             <  1:3764596/59‑1 (MQ=255)
                    aTCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttttt                                                             <  1:3595563/59‑1 (MQ=255)
                    aTCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttttt                                                             <  1:3333333/59‑1 (MQ=255)
                     tCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCActtct                                                <  1:2740018/71‑1 (MQ=255)
                     tCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCActtct                                                >  1:5613210/1‑71 (MQ=255)
                      cGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttt                                                               <  1:5295160/55‑1 (MQ=255)
                       gCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCa                                                     <  1:3957735/64‑1 (MQ=255)
                       gCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCActtctt                                               <  1:4692638/70‑1 (MQ=255)
                         aaaGGGATGAACCTGGCGTCGAACCTGTTCCATGTCGGGATTCTGGGGAt                                                                 >  1:4840042/1‑50 (MQ=255)
                         aaaGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCgg                                            <  1:1799919/71‑1 (MQ=255)
                           aGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTa                                          >  1:633968/1‑71 (MQ=255)
                            gggATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCg                                             >  1:4080490/1‑67 (MQ=255)
                            gggATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCg                                             >  1:1547980/1‑67 (MQ=255)
                              gATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTa                                          >  1:1127883/1‑68 (MQ=255)
                               aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTc                                                                        <  1:1394125/37‑1 (MQ=255)
                                 gAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa                                    >  1:4804580/1‑71 (MQ=255)
                                    ccTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACttcttc                                              <  1:5253222/58‑1 (MQ=255)
                                     cTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgcc                                <  1:1498757/71‑1 (MQ=255)
                                      tGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg                               >  1:1827984/1‑71 (MQ=255)
                                           tCGAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCAc                                                    <  1:1052377/45‑1 (MQ=16)
                                           tCGAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCAc                                                    <  1:5684795/45‑1 (MQ=16)
                                           tCGAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCAc                                                    <  1:2456482/45‑1 (MQ=16)
                                           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)
                                               aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTa                                          <  1:3078725/51‑1 (MQ=255)
                                               aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTa                                          <  1:4547925/51‑1 (MQ=255)
                                                  tGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTAt                   >  1:3813579/1‑71 (MQ=255)
                                                  tGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTAt                   >  1:3413289/1‑71 (MQ=255)
                                                   gTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGt                     <  1:1693949/68‑1 (MQ=255)
                                                    ttCCATATCGGGATTCTGGGGATTTTTGTCGGTCActtctt                                               >  1:2709699/1‑41 (MQ=255)
                                                    ttCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCgg                                            <  1:3601887/44‑1 (MQ=255)
                                                    ttCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCgg                                            <  1:2088660/44‑1 (MQ=255)
                                                    ttCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCAt                                         <  1:1660383/47‑1 (MQ=16)
                                                    ttCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCAt                                         <  1:4474413/47‑1 (MQ=16)
                                                      ccATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATgt                                       >  1:453948/1‑46 (MQ=255)
                                                        atatCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACCCCGCACTGGATg                      <  1:1901880/62‑1 (MQ=255)
                                                         tatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                      <  1:625380/45‑1 (MQ=255)
                                                         tatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                      <  1:1612111/45‑1 (MQ=255)
                                                         tatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGcc             <  1:2106917/70‑1 (MQ=255)
                                                         tatCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCgg                                            <  1:2237207/39‑1 (MQ=16)
                                                         tatCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATgtt                                      <  1:4838486/45‑3 (MQ=255)
                                                         tatCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATgtt                                      <  1:2480542/45‑3 (MQ=255)
                                                         tatCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATgtt                                      <  1:1282650/45‑3 (MQ=255)
                                                            cGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGc         >  1:4708655/1‑71 (MQ=255)
                                                                aTTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGt                     <  1:3749098/55‑1 (MQ=255)
                                                                aTTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGcc     <  1:3701736/71‑1 (MQ=255)
                                                                 ttCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                      <  1:5320668/37‑1 (MQ=255)
                                                                 ttCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                      <  1:208765/37‑1 (MQ=255)
                                                                 ttCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                      <  1:3504586/37‑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)
                                                                   |                                                                      
GGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAT  >  minE/837638‑837775

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: