Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 498,954 C→T Q430Q (CAG→CAA modF ← fused subunits of molybdate transporter and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE498,9540CT93.4% 186.6 / 3.9 61Q430Q (CAG→CAAmodFfused subunits of molybdate transporter and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (24/33);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGC  >  minE/498886‑499016
                                                                    |                                                              
gTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTg                                                              <  1:591563/71‑1 (MQ=255)
 tttCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCATCGGATCAAGCCCTTGt                                                             <  1:2280382/71‑1 (MQ=255)
 tttCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTg                                                              <  1:2376055/70‑1 (MQ=255)
 tttCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGt                                                             <  1:459024/71‑1 (MQ=255)
   tCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAg                                                           >  1:1124280/1‑71 (MQ=255)
    cACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAg                                                           <  1:917083/70‑1 (MQ=255)
    cACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGt                                                          >  1:59564/1‑71 (MQ=255)
     aCCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGt                                                          <  1:2396725/70‑1 (MQ=255)
        ttCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGtt                                                      >  1:2394701/1‑71 (MQ=255)
           gCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCa                                                    <  1:1745639/70‑1 (MQ=255)
           gCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCAt                                                   <  1:222156/71‑1 (MQ=255)
              aaTCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATc                                                  >  1:305181/1‑69 (MQ=255)
               aTCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATc                                                  >  1:61081/1‑68 (MQ=255)
               aTCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCga                                                <  1:127925/70‑1 (MQ=255)
               aTCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCga                                                <  1:125066/70‑1 (MQ=255)
               aTCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCga                                                <  1:1985769/70‑1 (MQ=255)
               aTCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCg                                                 <  1:1949386/69‑1 (MQ=255)
                tCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCgaga                                              <  1:2105326/71‑1 (MQ=255)
                  aGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGaa                                             >  1:221492/1‑70 (MQ=255)
                  aGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAAt                                            <  1:1679124/71‑1 (MQ=255)
                       aTCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAg                                         >  1:1482065/1‑69 (MQ=255)
                        tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAAt                                            <  1:150089/65‑1 (MQ=255)
                           acaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGt                                   <  1:219365/71‑1 (MQ=255)
                              aaaCGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCaa                                      >  1:1636075/1‑65 (MQ=255)
                              aaaCGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCgg                                <  1:919783/71‑1 (MQ=255)
                                aCGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCAt                                                   <  1:170521/50‑1 (MQ=255)
                                aCGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCAt                                                   <  1:415472/50‑1 (MQ=255)
                                aCGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCAt                                                   <  1:2201913/50‑1 (MQ=255)
                                  ggcggATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCg                                                 >  1:1114652/1‑50 (MQ=255)
                                   gcggATAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGtt                           <  1:2159382/71‑1 (MQ=255)
                                       aTAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCAc                       >  1:2361019/1‑71 (MQ=255)
                                       aTAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCAc                       >  1:2291874/1‑71 (MQ=255)
                                       aTAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCAc                       >  1:1036353/1‑71 (MQ=255)
                                        tAAGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCAc                       >  1:2064602/1‑70 (MQ=255)
                                          aGCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATc                                                  <  1:1842507/41‑1 (MQ=255)
                                           gCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCgaga                                              <  1:228287/44‑1 (MQ=255)
                                           gCTGGCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCg                                 <  1:250471/57‑1 (MQ=255)
                                              ggCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGt                                   <  1:418092/52‑1 (MQ=255)
                                              ggCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGt                                   <  1:237503/52‑1 (MQ=255)
                                               gCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTgcg                <  1:2038997/70‑1 (MQ=255)
                                               gCGATTCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTgc                 <  1:1980505/69‑1 (MQ=255)
                                               gCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATaa                                          <  1:1999394/44‑1 (MQ=255)
                                                 gATTCAGCGGATCAACCCCCTGTAGTGGTTCATCGAGAATAAGc                                        >  1:6447/1‑44 (MQ=255)
                                                   ttCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGt                                   >  1:184186/1‑47 (MQ=255)
                                                   ttCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGt                                   >  1:2316306/1‑47 (MQ=255)
                                                   ttCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAg                    <  1:834618/62‑1 (MQ=255)
                                                   ttCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAg                    <  1:1151220/62‑1 (MQ=255)
                                                    tCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGttt                          >  1:410612/1‑55 (MQ=255)
                                                    tCAGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGttt                          >  1:554907/1‑55 (MQ=255)
                                                      aGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACaa         <  1:2193431/70‑1 (MQ=255)
                                                      aGCGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAAt        <  1:2208218/71‑1 (MQ=255)
                                                        cGGATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAAt        >  1:993457/1‑69 (MQ=255)
                                                         ggATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGAt                              >  1:1901162/1‑46 (MQ=255)
                                                         ggATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCa                        >  1:1429030/1‑52 (MQ=255)
                                                         ggATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCAc                       >  1:776653/1‑53 (MQ=255)
                                                         ggATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACa          >  1:485663/1‑66 (MQ=255)
                                                         ggATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAAt        >  1:750571/1‑68 (MQ=255)
                                                         ggATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAg     <  1:2065259/71‑1 (MQ=255)
                                                          gATCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAg     >  1:1052617/1‑70 (MQ=255)
                                                            tCAAGCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAgcgc  <  1:1458559/71‑1 (MQ=255)
                                                                gCCCTTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTgcgc               <  1:1053787/54‑1 (MQ=255)
                                                                    |                                                              
GTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGC  >  minE/498886‑499016

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: