Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 757,001 G→T 8.8% A143S (GCG→TCG)  lolE → outer membrane‑specific lipoprotein transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE757,0010GT 8.8% 228.3 / 5.4 68A143S (GCG→TCG) lolEouter membrane‑specific lipoprotein transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (19/43);  new base T (3/3);  total (22/46)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAGGGGGATTG  >  minE/756933‑757068
                                                                    |                                                                   
cGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCg                                                                   <  1:1760918/71‑1 (MQ=255)
 gCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCgg                                                                  >  1:1015018/1‑71 (MQ=255)
    aCAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCg                                                               <  1:3289084/71‑1 (MQ=255)
    aCAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCg                                                               <  1:1833081/71‑1 (MQ=255)
      aGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAATCGGGCGaa                                                             <  1:84537/71‑1 (MQ=255)
      aGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGaa                                                             <  1:1018430/71‑1 (MQ=255)
           cAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAATCGGGCGaa                                                             <  1:3410259/66‑1 (MQ=255)
             gCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCg                                                               <  1:180105/62‑1 (MQ=255)
              cGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAAt                                                     <  1:3684380/71‑1 (MQ=255)
              cGTCTGAGCGCATTACCCTAGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAAt                                                     <  1:933661/71‑1 (MQ=255)
               ggcTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCg                                                               <  1:1562504/58‑1 (MQ=255)
                 cTGAGCGCATTACCCTCGTTTGTACAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGc                                                         >  1:1899133/1‑64 (MQ=255)
                     gcgcATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTatcatc                                              <  1:1234089/71‑1 (MQ=255)
                      cgcATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCg                                             <  1:3595031/71‑1 (MQ=255)
                       gcATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCgg                                            >  1:1248329/1‑71 (MQ=255)
                         aTTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTTAAGCGGGCGAACAGCAAATTATCATCGGCa                                          <  1:1797253/71‑1 (MQ=255)
                         aTTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCg                                             <  1:2167293/68‑1 (MQ=255)
                         aTTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCa                                          <  1:3299888/71‑1 (MQ=255)
                         aTTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCa                                          <  1:1167855/71‑1 (MQ=255)
                          ttACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCgg                                                                  <  1:1352648/46‑1 (MQ=255)
                          ttACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCgg                                            <  1:1758023/68‑1 (MQ=255)
                           tACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCa                                          <  1:1070573/69‑1 (MQ=255)
                                tCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCaaa                                                      >  1:633405/1‑52 (MQ=255)
                                  gTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCa                                                        <  1:2107283/48‑1 (MQ=255)
                                  gTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCa                                                        <  1:3467720/48‑1 (MQ=255)
                                  gTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTg                                  <  1:3441740/70‑1 (MQ=255)
                                    ttGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGaa                                                             >  1:81232/1‑41 (MQ=255)
                                    ttGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGaa                                                             >  1:2857015/1‑41 (MQ=255)
                                    ttGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGaa                                                             >  1:133422/1‑41 (MQ=255)
                                    ttGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGaa                                                             >  1:1553046/1‑41 (MQ=255)
                                    ttGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTatcatc                                              <  1:3432838/56‑1 (MQ=255)
                                     tGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCg                               <  1:2378149/70‑1 (MQ=255)
                                        ccAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTat                                                  <  1:2526918/48‑1 (MQ=255)
                                        ccAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTat                                                  <  1:2587511/48‑1 (MQ=255)
                                            ggCGATGCGTGGCGCAATTTTAAATCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCt                       >  1:2393580/1‑71 (MQ=255)
                                            ggCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCt                       >  1:293226/1‑71 (MQ=255)
                                            ggCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCt                       >  1:1078709/1‑71 (MQ=255)
                                              cGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGGAAAGGCGTGGCGGATGCGCTGa                     <  1:878856/71‑1 (MQ=255)
                                              cGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCtca                     <  1:3004704/71‑3 (MQ=255)
                                              cGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGa                     <  1:528863/71‑1 (MQ=255)
                                               gATGCGTGGCGCAATTTTAAATCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGaa                    <  1:2865324/71‑1 (MQ=255)
                                               gATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGaa                    <  1:2333835/71‑1 (MQ=255)
                                                aTGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGc                                     >  1:3340388/1‑53 (MQ=255)
                                                aTGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGaaa                   <  1:922169/71‑1 (MQ=255)
                                                 tGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTatca                                                <  1:14307/41‑1 (MQ=255)
                                                  gCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCaa                                         <  1:951033/47‑1 (MQ=255)
                                                  gCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCaa                                         <  1:418431/47‑1 (MQ=255)
                                                  gCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTgg                                 <  1:2364850/55‑1 (MQ=255)
                                                    gTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGa               <  1:3057235/71‑1 (MQ=255)
                                                       gcgcAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGa                     <  1:2685402/62‑1 (MQ=255)
                                                       gcgcAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAg             >  1:2935008/1‑70 (MQ=255)
                                                          cAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAgg         <  1:473636/71‑1 (MQ=255)
                                                           aaTTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAggg        >  1:3262136/1‑71 (MQ=255)
                                                           aaTTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAggg        <  1:3491116/71‑1 (MQ=255)
                                                            aTTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGa                     <  1:246180/57‑1 (MQ=255)
                                                            aTTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAgggg       >  1:1420697/1‑71 (MQ=255)
                                                            aTTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAgggg       >  1:1391446/1‑71 (MQ=255)
                                                             ttttAAATCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAgggg       >  1:2825616/1‑70 (MQ=255)
                                                             ttttAAATCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAgggg       >  1:1926723/1‑70 (MQ=255)
                                                               ttAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTg                      >  1:2548789/1‑53 (MQ=255)
                                                               ttAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTg                      >  1:1970880/1‑53 (MQ=255)
                                                                tAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAGGGGGAt    <  1:3302270/70‑1 (MQ=255)
                                                                tAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGCGGATGCGCTGAAAGTGAAGCAGGGGGAtt   <  1:541073/70‑1 (MQ=255)
                                                                 aaaGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAGGGGGAtt   <  1:858024/70‑1 (MQ=255)
                                                                 aaaGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAGGGGGATTg  >  1:2918983/1‑71 (MQ=255)
                                                                 aaaGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAGGGGGATTg  >  1:3012563/1‑71 (MQ=255)
                                                                  aaGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAGGGGGATTg  <  1:466265/70‑1 (MQ=255)
                                                                  aaGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAGGGGGATTg  <  1:1473184/70‑1 (MQ=255)
                                                                    gCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAGGGGGATTg  <  1:1733238/68‑1 (MQ=255)
                                                                    |                                                                   
CGCAACAGGAACAGCGTCTGAGCGCATTACCCTCGTTTGTCCAGGGCGATGCGTGGCGCAATTTTAAAGCGGGCGAACAGCAAATTATCATCGGCAAAGGCGTGGCGGATGCGCTGAAAGTGAAGCAGGGGGATTG  >  minE/756933‑757068

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