Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,069,493 G→A 17.9% M130I (ATG→ATA cfa → cyclopropane fatty acyl phospholipid synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,069,4930GA17.9% 186.1 / 31.2 78M130I (ATG→ATAcfacyclopropane fatty acyl phospholipid synthase
Reads supporting (aligned to +/- strand):  ref base G (41/23);  new base A (9/5);  total (50/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

GAGTAAAAAACGTGCCTGGATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATCTGGAATCTGCCCAGCA  >  minE/1069424‑1069561
                                                                     |                                                                    
gagTAAAAAACGTGCCTGGATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGc                                                                     <  1:1813654/71‑1 (MQ=255)
gagTAAAAAACGTGCCTGGATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGc                                                                     >  1:2460552/1‑71 (MQ=255)
   tAAAAAACGTGCCTGGATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCtt                                                                   >  1:1024703/1‑70 (MQ=255)
     aaaaaCGTGCCTGGATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATACTTGAt                                                                >  1:198769/1‑71 (MQ=255)
       aaaCGTGCCTGGATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCtt                                                                   >  1:2679405/1‑66 (MQ=255)
       aaaCGTGCCTGGATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATcc                                                              >  1:4105843/1‑71 (MQ=255)
        aaCGTGCCTGGATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATccc                                                             >  1:375520/1‑71 (MQ=255)
        aaCGTGCCTGGATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATACTTGATccc                                                             >  1:1704716/1‑71 (MQ=255)
         aCGTGCCTGGATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCt                                                            <  1:3418132/71‑1 (MQ=255)
               cTGGATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGc                                                      <  1:3799203/71‑1 (MQ=255)
                 ggATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCaa                                                    >  1:1717226/1‑71 (MQ=255)
                 ggATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCaa                                                    >  1:3928585/1‑71 (MQ=255)
                 ggATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCa                                                     >  1:2265775/1‑70 (MQ=255)
                  gATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGAt                                                                >  1:3750442/1‑58 (MQ=255)
                  gATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCaa                                                    >  1:1567167/1‑70 (MQ=255)
                  gATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCaa                                                    >  1:3438333/1‑70 (MQ=255)
                  gATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCaa                                                    >  1:870293/1‑70 (MQ=255)
                   aTAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCaa                                                    >  1:1019492/1‑69 (MQ=255)
                    tAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCtt                                                                   <  1:2835584/53‑1 (MQ=255)
                    tAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCtt                                                                   <  1:2830949/53‑1 (MQ=255)
                    tAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGAt                                                                >  1:1741468/1‑56 (MQ=255)
                      gTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTg                                                                  <  1:3197453/52‑1 (MQ=255)
                       tCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTcc                                              >  1:2892264/1‑71 (MQ=255)
                       tCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTcc                                              >  1:1020546/1‑71 (MQ=255)
                         ggCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCGGCATGCTTGATCCCTTCAt                                                        <  1:550297/59‑1 (MQ=255)
                         ggCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCtt                                                           <  1:2898610/56‑1 (MQ=255)
                         ggCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCa                                                         >  1:2817048/1‑58 (MQ=255)
                         ggCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTg                                            >  1:884685/1‑71 (MQ=255)
                          gCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGc                                                      >  1:2539036/1‑60 (MQ=255)
                            aaaGCGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATACTTGATCCCTTCATGCAATATTCCTgcgc                                         >  1:2043033/1‑71 (MQ=255)
                              agagCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCt                                        >  1:658093/1‑70 (MQ=255)
                              agagCATTACGATTTGGGTAATGACTTGTTCAGCCGCATACTTGATCCCTTCATGCAATAtt                                                <  1:48601/62‑1 (MQ=255)
                              agagCATTACGATTTGGGTAATGACTTGTTCAGCCGCATACTTGATCCCTTCATGCAATAtt                                                <  1:2056438/62‑1 (MQ=255)
                               gagCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTa                                      >  1:408425/1‑71 (MQ=255)
                                 gCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACt                                    >  1:1713244/1‑71 (MQ=255)
                                 gCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACt                                    <  1:3193848/71‑1 (MQ=255)
                                  cATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTg                                   >  1:880099/1‑71 (MQ=255)
                                   aTTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTg                                   >  1:94834/1‑70 (MQ=255)
                                    ttACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCt                                        <  1:2212272/64‑1 (MQ=255)
                                    ttACGATTTGGGTAATGACTTGTTCAGCCGCATACTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGa                                 <  1:4135431/71‑1 (MQ=255)
                                      aCGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGAt                                                                >  1:1208878/1‑38 (MQ=255)
                                      aCGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGAt                                                                >  1:2243402/1‑38 (MQ=255)
                                      aCGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATcc                                                              >  1:195264/1‑40 (MQ=255)
                                      aCGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCAt                                                        >  1:2434509/1‑46 (MQ=255)
                                      aCGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGaa                                <  1:854542/70‑1 (MQ=255)
                                      aCGATTTGGGTAATGACTTGTTCAGCCGCATACTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGa                                 >  1:2273422/1‑69 (MQ=255)
                                       cGATTTGGGTAATGACTTGTTCAGCCGCATACTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGaaa                               <  1:3052750/70‑1 (MQ=255)
                                        gATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCtt                                       <  1:2507826/61‑1 (MQ=255)
                                        gATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGaaa                               >  1:1050951/1‑69 (MQ=255)
                                        gATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGa                                 >  1:1809791/1‑67 (MQ=255)
                                        gATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGa                                 >  1:1408941/1‑67 (MQ=255)
                                        gATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGa                             <  1:119908/71‑1 (MQ=255)
                                           ttGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGc                          >  1:313005/1‑71 (MQ=255)
                                            tGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTa                                      >  1:2354589/1‑58 (MQ=255)
                                             gggTAATGACTTGTTCAGCCGCATACTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGcc                         <  1:3826405/70‑1 (MQ=255)
                                                tAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATa                     <  1:3231190/71‑1 (MQ=255)
                                                 aaTGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCt                                        <  1:462049/51‑1 (MQ=255)
                                                  aTGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAAt                   >  1:1342467/1‑71 (MQ=255)
                                                  aTGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAAt                   >  1:677143/1‑71 (MQ=255)
                                                  aTGACTTGTTCAGCCGCATACTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAAt                   >  1:3791019/1‑71 (MQ=255)
                                                   tGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATc                  >  1:2173444/1‑71 (MQ=255)
                                                   tGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATc                  <  1:803152/71‑1 (MQ=255)
                                                    gACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAat                                                   >  1:366113/1‑37 (MQ=255)
                                                     aCTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTc                                               >  1:3706786/1‑40 (MQ=255)
                                                     aCTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAAt                   >  1:3882050/1‑68 (MQ=255)
                                                     aCTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATCTg                >  1:2481163/1‑71 (MQ=255)
                                                     aCTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATCTg                <  1:846140/71‑1 (MQ=255)
                                                     aCTTGTTCAGCCGCATGCTTGATCCCTTAATGCAATATTc                                               >  1:2549616/1‑40 (MQ=255)
                                                       ttGTTCAGACGCATACTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATCTgg               >  1:2147141/1‑70 (MQ=255)
                                                         gTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATCTGGa              <  1:2193956/69‑1 (MQ=255)
                                                         gTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATCTGGa              <  1:3633145/69‑1 (MQ=255)
                                                         gTTCAGCCGCATACTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATc                  >  1:3530044/1‑65 (MQ=255)
                                                             aGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATCTGGAATCTg         <  1:3047708/70‑1 (MQ=255)
                                                             aGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATCTGGAATCTg         <  1:764922/70‑1 (MQ=255)
                                                               ccGCATACTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATCTGGAATCTGcc       >  1:822987/1‑70 (MQ=255)
                                                               ccGCATACTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATCTGGAATCTGc        >  1:1937018/1‑69 (MQ=255)
                                                                  cATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCTATAATCTGGAATCTGCCcagc   <  1:2723030/71‑1 (MQ=255)
                                                                   aTGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATCTGGAATCTGCCcagca  <  1:3872477/71‑1 (MQ=255)
                                                                     |                                                                    
GAGTAAAAAACGTGCCTGGATAGTCGGCAAAGAGCATTACGATTTGGGTAATGACTTGTTCAGCCGCATGCTTGATCCCTTCATGCAATATTCCTGCGCTTACTGGAAAGATGCCGATAATCTGGAATCTGCCCAGCA  >  minE/1069424‑1069561

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: