Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 756,787 G→A 57.6% V71V (GTG→GTA lolE → outer membrane‑specific lipoprotein transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE756,7870GA57.6% 17.9 / 48.8 66V71V (GTG→GTAlolEouter membrane‑specific lipoprotein transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (17/11);  new base A (30/8);  total (47/19)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGG  >  minE/756725‑756843
                                                              |                                                        
tttGAACGCTAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:1273454/1‑69 (MQ=255)
tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:561060/1‑69 (MQ=255)
tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:2419856/1‑69 (MQ=255)
tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:1277254/1‑69 (MQ=255)
tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:3339514/1‑69 (MQ=255)
tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:2570872/1‑69 (MQ=255)
tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:2574019/1‑69 (MQ=255)
tttGAACGCGAACTGAACAACCGTATTCTGGCGGTGGCGCCGCATGGTGAAATCGAGGCGGTGGATCAg                                                    >  1:1771270/1‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGATCAACTGGc                                    >  1:1962097/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:3205318/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:1851772/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:2480785/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:2485396/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:2496292/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:1726757/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:1604241/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:3097310/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:3348333/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:3369553/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:1370479/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:3395298/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:431358/3‑69 (MQ=255)
                ataacCGCATTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGc                                    >  1:643823/3‑69 (MQ=255)
                 caacCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGc                                    >  1:2692861/1‑68 (MQ=255)
                 caacCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGc                                    >  1:92325/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:2169032/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:799365/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:573864/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:48117/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:118646/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1194334/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1268327/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1412418/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1446554/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1549970/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:2000885/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1859803/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:2176237/2‑69 (MQ=255)
                       caTTCTGGCGGTAGTACCGCACGGTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGc                             >  1:1766529/2‑69 (MQ=255)
                           cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACt                          >  1:1086869/1‑68 (MQ=255)
                           cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         >  1:1828168/1‑69 (MQ=255)
                           cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:386464/69‑1 (MQ=255)
                           cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         >  1:352955/1‑69 (MQ=255)
                           cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:1758355/69‑1 (MQ=255)
                           cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:2415337/69‑1 (MQ=255)
                           cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:254092/69‑1 (MQ=255)
                           cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:1599146/69‑1 (MQ=255)
                           cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:1148013/69‑1 (MQ=255)
                           cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:3068324/69‑1 (MQ=255)
                           cTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTg                         <  1:2940535/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:353259/69‑1 (MQ=255)
                                              gTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGCACTGGATAACGTGCAGAAAGTGc      <  1:1442706/69‑1 (MQ=255)
                                              gTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGCACTGGATAACGTGCAGAAAGTGc      <  1:3147865/69‑1 (MQ=255)
                                              gTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGCACTGGATAACGTGCAGAAAGTGc      <  1:1134456/69‑1 (MQ=255)
                                              gTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGCACTGGATAACGTGCAGAAAGTGc      <  1:2208094/69‑1 (MQ=255)
                                              gTGAAATTGAGGCGGTAGATCAGCCGTGGACCAACTGGCAGGAAGCACTGGATAACGTGCAGAAAGTGc      <  1:1571300/69‑1 (MQ=255)
                                                   aTCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCggg  <  1:3047302/68‑1 (MQ=255)
                                                   aTCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCggg  <  1:1013386/68‑1 (MQ=255)
                                                              |                                                        
TTTGAACGCGAACTGAACAACCGTATTCTGGCGGTGGTGCCGCATGGTGAAATCGAGGCGGTGGATCAGCCGTGGACTAACTGGCAGGAAGCACTGGATCACGTGCAAAAAGTGCCGGG  >  minE/756725‑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: