Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 551,606 G→A 15.0% L640L (CTG→TTG)  ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE551,6060GA15.0% 180.6 / 20.8 73L640L (CTG→TTG) ybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base A (5/6);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

ATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTT  >  minE/551537‑551675
                                                                     |                                                                     
atgaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGt                                                                      >  1:568524/1‑71 (MQ=255)
 tgaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGt                                                                      <  1:3287217/70‑1 (MQ=255)
   aCGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcg                                                                   >  1:3199100/1‑71 (MQ=255)
     gTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgtc                                                                 >  1:2841901/1‑71 (MQ=255)
     gTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgtc                                                                 >  1:690985/1‑71 (MQ=255)
       aCCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                                               >  1:2461637/1‑71 (MQ=255)
        cccGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                                               <  1:1406341/70‑1 (MQ=255)
             gTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgtcg                                                                <  1:1929855/64‑1 (MQ=255)
               aaTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTc                                                       >  1:299013/1‑71 (MQ=255)
                aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                                                      >  1:820950/1‑71 (MQ=255)
                aTAGTTGCCGCCAACCGGACAACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTc                                                       >  1:1056025/1‑70 (MQ=255)
                 tAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTg                                                     <  1:3754005/71‑1 (MQ=255)
                   gTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGAt                                                   <  1:2336962/71‑1 (MQ=255)
                    ttGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtcgtc                                                                 <  1:1872860/56‑1 (MQ=255)
                     tGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGAt                                                   >  1:546058/1‑69 (MQ=255)
                     tGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATa                                                  <  1:3023056/70‑1 (MQ=255)
                      gccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATaa                                                 <  1:2269083/70‑1 (MQ=255)
                      gacgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGa                                                    >  1:797005/3‑67 (MQ=255)
                       ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAg                                                <  1:2839973/70‑1 (MQ=255)
                       ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtcgtc                                                                 >  1:2531550/1‑53 (MQ=255)
                       ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtcgtc                                                                 >  1:1226436/1‑53 (MQ=255)
                         gccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgc                                             <  1:630789/71‑1 (MQ=255)
                         gccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgc                                             >  1:1563408/1‑71 (MQ=255)
                         gccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgc                                             <  1:266129/71‑1 (MQ=255)
                          ccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgc                                             >  1:896987/1‑70 (MQ=255)
                          ccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgc                                               >  1:2072533/1‑68 (MQ=255)
                           cAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgc                                             >  1:3929534/1‑69 (MQ=255)
                            aaCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAg                                          >  1:1799281/1‑71 (MQ=255)
                            aaCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAg                                          >  1:965548/1‑71 (MQ=255)
                             aCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGc                                         >  1:1636360/1‑71 (MQ=255)
                              ccGGACCACGACCGTAGCGCGTATTATGGTACTGTTTCAGTTCGTCGATATAg                                                          >  1:220235/1‑53 (MQ=255)
                              ccGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGc                                         >  1:768676/1‑70 (MQ=255)
                                ggACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATaa                                                 >  1:37049/1‑60 (MQ=255)
                                 gACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgtc                                                                 <  1:1792764/43‑1 (MQ=255)
                                  aCCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagca                                    <  1:1239914/71‑1 (MQ=255)
                                  aCCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagca                                    <  1:1120119/71‑1 (MQ=255)
                                     aCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGa                                           >  1:3559369/1‑61 (MQ=255)
                                      cGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGatat                                                            >  1:451259/1‑43 (MQ=255)
                                      cGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATaa                                                 >  1:2553388/1‑54 (MQ=255)
                                      cGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt                                <  1:480832/71‑1 (MQ=255)
                                       gACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa                               <  1:1679345/71‑1 (MQ=255)
                                        aCCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATa                                                  <  1:3407214/51‑1 (MQ=255)
                                         ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATc                             >  1:231580/1‑71 (MQ=255)
                                           gTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                           <  1:212669/71‑1 (MQ=255)
                                           gTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                           <  1:3637067/71‑1 (MQ=255)
                                            tAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                           >  1:988280/1‑70 (MQ=255)
                                            tAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                           >  1:2956162/1‑70 (MQ=255)
                                            tAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                           >  1:1599210/1‑70 (MQ=255)
                                              gcgcGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACa                         <  1:553460/70‑1 (MQ=255)
                                              gcgcGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGAAGCGTATCGACAg                        <  1:3852110/71‑1 (MQ=255)
                                               cgcgGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAg                        <  1:1513443/70‑1 (MQ=255)
                                               cgcgGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGt                       <  1:2121937/71‑1 (MQ=255)
                                               cgcgGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGa                                           <  1:307639/51‑1 (MQ=255)
                                                  ggATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTAt                              >  1:3364302/1‑61 (MQ=255)
                                                   gATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTatca                   >  1:149624/1‑71 (MQ=255)
                                                           tACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGtt           <  1:1118933/71‑1 (MQ=255)
                                                           tACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGtt           <  1:3392097/71‑1 (MQ=255)
                                                           tACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGtt           <  1:765023/71‑1 (MQ=255)
                                                           tACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGtt           >  1:638977/1‑71 (MQ=255)
                                                            aCTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                           >  1:36055/1‑54 (MQ=255)
                                                              tGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTatca                   >  1:2315831/1‑60 (MQ=255)
                                                              tGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGcc        >  1:894711/1‑71 (MQ=255)
                                                               gTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACa                         >  1:719086/1‑53 (MQ=255)
                                                               gTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                           <  1:674927/51‑1 (MQ=255)
                                                               gTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                           <  1:1593032/51‑1 (MQ=255)
                                                                tttCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACa                         <  1:394473/52‑1 (MQ=255)
                                                                tttCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCg       <  1:2749467/70‑1 (MQ=255)
                                                                tttCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCg       <  1:2754120/70‑1 (MQ=255)
                                                                tttCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCg       <  1:3131170/70‑1 (MQ=255)
                                                                tttCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGt      <  1:3256480/71‑1 (MQ=255)
                                                                 ttCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGt      >  1:1889762/1‑70 (MQ=255)
                                                                 ttCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTa     >  1:2448938/1‑71 (MQ=255)
                                                                    aGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACtt  <  1:1307973/71‑1 (MQ=255)
                                                                     |                                                                     
ATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTT  >  minE/551537‑551675

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: