Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 373,682 G→T 6.3% A178S (GCG→TCG)  fepE → regulator of length of O‑antigen component of lipopolysaccharide chains

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE373,6820GT 6.3% 317.5 / 4.9 94A178S (GCG→TCG) fepEregulator of length of O‑antigen component of lipopolysaccharide chains
Reads supporting (aligned to +/- strand):  ref base G (53/35);  new base T (3/3);  total (56/38)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

ACAATGCCAGTAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTATATCTCTACGTTGGTGGTGAA  >  minE/373620‑373752
                                                              |                                                                      
aCAATGCCAGTAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAAc                                                                 <  1:4379779/70‑1 (MQ=255)
  aaTGCCAGTAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAg                                                              >  1:938849/1‑71 (MQ=255)
  aaTGCCAGTAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAg                                                              >  1:3612599/1‑71 (MQ=255)
   aTGCCAGTAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGt                                                             >  1:1596492/1‑71 (MQ=255)
      ccAGTAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGaa                                                          <  1:2342775/71‑1 (MQ=255)
          tAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGcc                                                                    <  1:3404709/57‑1 (MQ=255)
          tAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGcc                                                                    <  1:2240228/57‑1 (MQ=255)
          tAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGa                                                           <  1:1760517/66‑1 (MQ=255)
          tAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGa                                                           <  1:1437887/66‑1 (MQ=255)
            aGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAag                                                         >  1:109609/1‑66 (MQ=255)
            aGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAag                                                         >  1:1772854/1‑66 (MQ=255)
            aGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGgcg                                                    >  1:4644084/1‑71 (MQ=255)
            aGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGgcg                                                    >  1:1576524/1‑71 (MQ=255)
            aGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGgcg                                                    >  1:3731025/1‑71 (MQ=255)
            aGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGgcg                                                    >  1:58897/1‑71 (MQ=255)
             gAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGgcg                                                    >  1:3322123/1‑70 (MQ=255)
               aaaaaGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCa                                                  <  1:3661420/70‑1 (MQ=255)
                    gATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGt                                            <  1:1449248/71‑1 (MQ=255)
                    gATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGt                                            <  1:1543279/71‑1 (MQ=255)
                    gATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGt                                            <  1:4378999/71‑1 (MQ=255)
                    gATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGt                                            <  1:640976/71‑1 (MQ=255)
                      tGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAAcc                                                                <  1:2569378/49‑1 (MQ=255)
                        aaCCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTg                                        <  1:2020203/71‑1 (MQ=255)
                         aCCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTg                                        <  1:3133432/70‑1 (MQ=255)
                         aCCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTg                                        <  1:3546437/70‑1 (MQ=255)
                         aCCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGa                                       <  1:1369877/71‑1 (MQ=255)
                          ccGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAg                                      <  1:3297352/71‑1 (MQ=255)
                          ccGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAg                                      <  1:707873/71‑1 (MQ=255)
                          ccGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAg                                      <  1:615491/71‑1 (MQ=255)
                             tCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGCCCGTTTTGAGCgg                                   >  1:565770/1‑71 (MQ=255)
                             tCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGa                                       >  1:2481968/1‑67 (MQ=255)
                             tCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCgg                                   >  1:4286484/1‑71 (MQ=255)
                             tCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCgg                                   >  1:20444/1‑71 (MQ=255)
                                cTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGgcg                                                    <  1:1230508/51‑1 (MQ=255)
                                cTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGgcg                                                    <  1:904526/51‑1 (MQ=255)
                                cTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGgcg                                                    <  1:1033828/51‑1 (MQ=255)
                                cTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGgcg                                                    <  1:2453952/51‑1 (MQ=255)
                                cTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGgcg                                                    <  1:3643662/51‑1 (MQ=255)
                                 tGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGaa                                                          >  1:3829657/1‑44 (MQ=255)
                                 tGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGaa                                                          >  1:988523/1‑44 (MQ=255)
                                 tGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGaa                                                          >  1:3569533/1‑44 (MQ=255)
                                 tGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGaa                                                          >  1:2578221/1‑44 (MQ=255)
                                 tGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGaa                                                          >  1:15024/1‑44 (MQ=255)
                                 tGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGaa                                                          >  1:2500364/1‑44 (MQ=255)
                                 tGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGaa                                                          >  1:1184084/1‑44 (MQ=255)
                                   tataCCTCCTGGACGCTAAGTTTTACCTCGCCAACCAGTGAAGAGGCGCAGACCGt                                            >  1:4000525/1‑56 (MQ=255)
                                   tataCCTCCTGGACGCTAAGTTTTACCTCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTatat                             >  1:1594792/1‑71 (MQ=255)
                                   tataCCTCCTGGACGCTAAGTTTTACCTCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTatat                             >  1:5044160/1‑71 (MQ=255)
                                   tataCCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTatat                             >  1:594630/1‑71 (MQ=255)
                                     taCCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAaga                                                        >  1:2313734/1‑42 (MQ=255)
                                     taCCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAaga                                                        >  1:4245293/1‑42 (MQ=255)
                                     taCCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAaga                                                        >  1:4374766/1‑42 (MQ=255)
                                     taCCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAaga                                                        >  1:4408883/1‑42 (MQ=255)
                                     taCCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAaga                                                        >  1:3576534/1‑42 (MQ=255)
                                     taCCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAaga                                                        >  1:2248681/1‑42 (MQ=255)
                                     taCCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAaga                                                        >  1:5043498/1‑42 (MQ=255)
                                     taCCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAaga                                                        >  1:1048385/1‑42 (MQ=255)
                                     taCCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAag                                                         >  1:3479299/1‑41 (MQ=255)
                                       cctcctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGAt                         >  1:5025839/1‑71 (MQ=255)
                                       cctcctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGAt                         >  1:1341546/1‑71 (MQ=255)
                                       cctcctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGAt                         >  1:1199031/1‑71 (MQ=255)
                                       cctcctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGAt                         >  1:1918316/1‑71 (MQ=255)
                                       cctcctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGAt                         >  1:3095979/1‑71 (MQ=255)
                                       cctcctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGAt                         >  1:3900618/1‑71 (MQ=255)
                                        ctcctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGAt                         >  1:107499/1‑70 (MQ=255)
                                         tcctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGAt                         <  1:4115847/69‑1 (MQ=255)
                                         tcctGGACGCTAAGTTTTACCGCGCCAACCAATGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGAt                         <  1:2500763/69‑1 (MQ=255)
                                           ctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATtata                     >  1:1025875/1‑71 (MQ=255)
                                           ctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATtata                     >  1:1133023/1‑71 (MQ=255)
                                           ctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATtata                     >  1:2849340/1‑71 (MQ=255)
                                           ctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATtat                      <  1:4786096/70‑1 (MQ=255)
                                           ctGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATtat                      <  1:3627810/70‑1 (MQ=255)
                                              gACGCTAAGTTTTACCTCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTatat                    <  1:4844735/69‑1 (MQ=255)
                                              gACGCTAAGTTTTACCTCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTatat                    <  1:2164970/69‑1 (MQ=255)
                                              gACGCTAAGTTTTACCTCGCCAACCAATGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTatat                    <  1:1564589/69‑1 (MQ=255)
                                              gACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTatat                    >  1:3036361/1‑69 (MQ=255)
                                              gACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTATAtct                  >  1:1599833/1‑71 (MQ=255)
                                              gACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTATAtc                   >  1:5191253/1‑70 (MQ=255)
                                                cGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGtt                                           <  1:4330211/44‑1 (MQ=255)
                                                cGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGtt                                           <  1:5170908/44‑1 (MQ=255)
                                                     aGTTTTACCGCGCCAACCGGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTatat                    >  1:4688263/1‑62 (MQ=255)
                                                     aGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCTTTTTGAGCGGGTATATCGATTATATCTCTACGtt           <  1:3314932/71‑1 (MQ=255)
                                                     aGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTatat                    >  1:3295924/1‑62 (MQ=255)
                                                     aGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTatat                    >  1:4291392/1‑62 (MQ=255)
                                                     aGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTatat                    >  1:133636/1‑62 (MQ=255)
                                                     aGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTATATCTCTACGtt           <  1:1799457/71‑1 (MQ=255)
                                                      gTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTATATCTCTACGtt           <  1:2363098/70‑1 (MQ=255)
                                                       ttttACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTATATCTCTACGTtgg         <  1:3406397/71‑1 (MQ=255)
                                                         ttAGCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGAt                         <  1:2026365/53‑1 (MQ=255)
                                                         ttACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGa                          >  1:3565523/1‑52 (MQ=255)
                                                         ttACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGa                          >  1:4944753/1‑52 (MQ=255)
                                                         ttACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGa                          >  1:1167888/1‑52 (MQ=255)
                                                         ttACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGa                          >  1:3630533/1‑52 (MQ=255)
                                                         ttACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGa                          >  1:3034865/1‑52 (MQ=255)
                                                             cgcgCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTATATCTCTACGTTGGTGGTGa   <  1:4947270/71‑1 (MQ=255)
                                                             cgcgCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTATATCTCTACGTTGGTGGTGa   <  1:3382250/71‑1 (MQ=255)
                                                              gcgcCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTATATCTCTACGTTGGTGGTGaa  <  1:2392641/71‑1 (MQ=255)
                                                              gcgcCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTATATCTCTACGTTGGTGGTGaa  <  1:4381882/71‑1 (MQ=255)
                                                              |                                                                      
ACAATGCCAGTAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCTAAGTTTTACCGCGCCAACCAGTGAAGAGGCGCAGACCGTTTTGAGCGGGTATATCGATTATATCTCTACGTTGGTGGTGAA  >  minE/373620‑373752

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: