Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 756,778 C→T 52.0% I68I (ATC→ATT lolE → outer membrane‑specific lipoprotein transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE756,7780CT52.0% ‑1.7 / 67.2 75I68I (ATC→ATTlolEouter membrane‑specific lipoprotein transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (25/11);  new base T (31/8);  total (56/19)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCGCGATGAACGGCTTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGG  >  minE/756711‑756843
                                                                   |                                                                 
gcgcGATGAACGGCTTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATc                                                                   <  1:2745262/68‑1 (MQ=255)
 cgcgATGAACGGCTTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGGAATCGa                                                                 >  1:2960978/1‑69 (MQ=255)
 cgcgATGAACGGCTTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGa                                                                 >  1:2979784/1‑69 (MQ=255)
 cgcgATGAACGGCTTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGa                                                                 >  1:1931116/1‑69 (MQ=255)
 cgcgATGAACGGCTTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGa                                                                 >  1:3107532/1‑69 (MQ=255)
 cgcgATGAACGGCTTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGa                                                                 >  1:1978506/1‑69 (MQ=255)
 cgcgATGAACGGCTTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGa                                                                 >  1:1419040/1‑69 (MQ=255)
 cgcgATGAACGGCTTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGa                                                                 >  1:42839/1‑69 (MQ=255)
 cgcgATGAACGGCTTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGa                                                                 >  1:523449/1‑69 (MQ=255)
              tttGAACGCTAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:1273454/1‑69 (MQ=255)
              tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:1277254/1‑69 (MQ=255)
              tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:2570872/1‑69 (MQ=255)
              tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:2574019/1‑69 (MQ=255)
              tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:2419856/1‑69 (MQ=255)
              tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:561060/1‑69 (MQ=255)
              tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:3339514/1‑69 (MQ=255)
              tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGCGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:1771270/1‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAg                                                                >  1:3057622/3‑41 (MQ=16)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGATCAACTGGc                                    >  1:1962097/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:2496292/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:1851772/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:2485396/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:2480785/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:1604241/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:1726757/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:3097310/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:3205318/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:3348333/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:3369553/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:3395298/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:1370479/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:431358/3‑69 (MQ=255)
                              ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:643823/3‑69 (MQ=255)
                               caacCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGc                                    >  1:92325/1‑68 (MQ=255)
                               caacCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGc                                    >  1:2692861/1‑68 (MQ=255)
                               caacCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGc                                    >  1:463153/1‑68 (MQ=255)
                               caacCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGc                                    >  1:535677/1‑68 (MQ=255)
                               caacCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGc                                    >  1:1641517/1‑68 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAg                              >  1:574101/2‑68 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1766529/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1549970/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1446554/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:573864/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:799365/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1412418/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1268327/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1859803/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:2000885/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1194334/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:2169032/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:2176237/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:48117/2‑69 (MQ=255)
                                     caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:118646/2‑69 (MQ=255)
                                         cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACt                          >  1:1086869/1‑68 (MQ=255)
                                         cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:1148013/69‑1 (MQ=255)
                                         cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:386464/69‑1 (MQ=255)
                                         cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         >  1:352955/1‑69 (MQ=255)
                                         cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:1599146/69‑1 (MQ=255)
                                         cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:3068324/69‑1 (MQ=255)
                                         cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:1758355/69‑1 (MQ=255)
                                         cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         >  1:1828168/1‑69 (MQ=255)
                                         cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:2940535/69‑1 (MQ=255)
                                         cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:254092/69‑1 (MQ=255)
                                         cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:2415337/69‑1 (MQ=255)
                                         cTGGCGGTAGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         >  1:351712/1‑69 (MQ=255)
                                          tGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:1576130/68‑1 (MQ=255)
                                                            gTGAAATTGAGGCGGTAGATCAGCCGTTGACCAACTGGCAGGAAGCACTGGATAACGTGCAGAAAGTGc      <  1:2656106/69‑1 (MQ=255)
                                                            gTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGCACTGGATAACGTGCAGTAAGTGc      <  1:439973/69‑1 (MQ=255)
                                                            gTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGCACTGGATAACGTGCAGAAAGTGc      <  1:1134456/69‑1 (MQ=255)
                                                            gTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGCACTGGATAACGTGCAGAAAGTGc      <  1:2208094/69‑1 (MQ=255)
                                                            gTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGCACTGGATAACGTGCAGAAAGTGc      <  1:353259/69‑1 (MQ=255)
                                                            gTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGCACTGGATAACGTGCAGAAAGTGc      <  1:1442706/69‑1 (MQ=255)
                                                            gTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGCACTGGATAACGTGCAGAAAGTGc      <  1:1571300/69‑1 (MQ=255)
                                                            gTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGCACTGGATAACGTGCAGAAAGTGc      <  1:3147865/69‑1 (MQ=255)
                                                                 aTCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCggg  <  1:3047302/68‑1 (MQ=255)
                                                                 aTCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCggg  <  1:1013386/68‑1 (MQ=255)
                                                                   |                                                                 
GCGCGATGAACGGCTTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGG  >  minE/756711‑756843

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: