Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 861,111 G→A 10.4% L640L (CTG→TTG)  ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048861,1110GA10.4% 234.0 / 13.6 77L640L (CTG→TTG) ybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base G (24/45);  new base A (4/4);  total (28/49)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTC  >  AP009048/861044‑861181
                                                                   |                                                                      
gaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtc                                                                     >  1:3768545/1‑71 (MQ=255)
gaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtc                                                                     >  1:3195879/1‑71 (MQ=255)
     aCCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                                                >  1:390344/1‑71 (MQ=255)
       ccGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                              <  1:3529165/71‑1 (MQ=255)
       ccGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                              <  1:2558341/71‑1 (MQ=255)
        cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGatat                                                             <  1:3367507/71‑1 (MQ=255)
          cgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAg                                                           <  1:3002994/71‑1 (MQ=255)
          cgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAg                                                           >  1:1833014/1‑71 (MQ=255)
            tAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGt                                                         >  1:3256275/1‑71 (MQ=255)
             aaTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTc                                                        >  1:504897/1‑71 (MQ=255)
              aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                                                       <  1:906763/71‑1 (MQ=255)
               tAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                                                       >  1:984924/1‑70 (MQ=255)
               tAGTTGACGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTg                                                      <  1:2061960/71‑1 (MQ=255)
                 gTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGAtata                                                            >  1:2456791/1‑63 (MQ=255)
                 gTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGAt                                                    <  1:3879246/71‑1 (MQ=255)
                  ttGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTTGATATTGGTCTGATa                                                   <  1:1626353/71‑1 (MQ=255)
                  ttGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                                                       <  1:1837963/67‑1 (MQ=255)
                  ttGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATa                                                   <  1:1916990/71‑1 (MQ=255)
                  ttGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATa                                                   <  1:1275690/71‑1 (MQ=255)
                  ttGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGa                                                                <  1:1821212/58‑1 (MQ=255)
                   ttccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATaa                                                  <  1:3625976/69‑1 (MQ=255)
                    gccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATaa                                                  <  1:3046347/70‑1 (MQ=255)
                     ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATaa                                                  <  1:2431432/69‑1 (MQ=255)
                        ccAACCGGACCACGACTGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGAt                                                    <  1:663312/64‑1 (MQ=255)
                         cAACGGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGt                                                         <  1:1156506/58‑1 (MQ=255)
                         cAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGcgcg                                             <  1:131758/70‑1 (MQ=255)
                         cAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGa                                            <  1:3583355/71‑1 (MQ=255)
                          aaCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAg                                           >  1:1378833/1‑71 (MQ=255)
                          aaCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAg                                           >  1:1193506/1‑71 (MQ=255)
                           aCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGc                                          >  1:2242999/1‑71 (MQ=255)
                           aCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGc                                          <  1:816599/71‑1 (MQ=255)
                           aCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTAGTCGATATAGGTCTGATAAGCGCGAGc                                          >  1:3211280/1‑71 (MQ=255)
                             cGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                                                <  1:1916978/47‑1 (MQ=255)
                             cGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGc                                          >  1:1416510/1‑69 (MQ=255)
                             cGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCca                                        <  1:2601823/71‑1 (MQ=255)
                              ggACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcag                                       <  1:2666780/71‑1 (MQ=255)
                               gACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGt                                                                       <  1:95220/38‑1 (MQ=255)
                               gACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcag                                       >  1:1137198/1‑70 (MQ=255)
                                aCCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                                                <  1:2882320/44‑1 (MQ=255)
                                 ccACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                              >  1:4384140/1‑45 (MQ=255)
                                 ccACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                              >  1:1857583/1‑45 (MQ=255)
                                 ccACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCcagcag                                    >  1:1723073/1‑71 (MQ=255)
                                   aCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg                                  <  1:4167001/71‑1 (MQ=255)
                                    cGATCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                                                >  1:2764122/1‑40 (MQ=255)
                                    cGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCATCAGCGt                                 <  1:4187210/71‑1 (MQ=255)
                                    cGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt                                 <  1:3844294/71‑1 (MQ=255)
                                    cGACAGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt                                 <  1:315362/71‑1 (MQ=255)
                                     gACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa                                >  1:3409938/1‑71 (MQ=255)
                                      aCCGTAGCGCGGATTATGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgct                                             <  1:1708916/56‑2 (MQ=255)
                                       ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                              >  1:4215939/1‑39 (MQ=255)
                                       ccGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATc                              <  1:1400007/71‑1 (MQ=255)
                                         gTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                            <  1:3347149/71‑1 (MQ=255)
                                            gcgcGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACa                          <  1:2042956/70‑1 (MQ=255)
                                            gcgcGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAg                         >  1:4481664/1‑71 (MQ=255)
                                            gcgcGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACa                          <  1:2242564/70‑1 (MQ=255)
                                                ggATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTat                      >  1:381740/1‑70 (MQ=255)
                                                 gATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTatca                    <  1:771253/71‑1 (MQ=255)
                                                  aTTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTatcat                   <  1:3040909/71‑1 (MQ=255)
                                                  aTTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTatcat                   <  1:276061/71‑1 (MQ=255)
                                                   ttATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTatcatc                  >  1:3982577/1‑71 (MQ=255)
                                                   ttATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTatcatc                  <  1:3620446/71‑1 (MQ=255)
                                                      tGGTACTGTTTCAGTTCGTCGATATAGGTCTGCTAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgtc               <  1:1960569/71‑1 (MQ=255)
                                                      tGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgt                <  1:447600/70‑1 (MQ=255)
                                                        gTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgtcg              <  1:116310/70‑1 (MQ=255)
                                                         tACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGa                                            <  1:268216/39‑1 (MQ=255)
                                                         tACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTa                       <  1:161664/60‑1 (MQ=255)
                                                         tACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGtt            <  1:2774723/71‑1 (MQ=255)
                                                         tACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCCTATCGACAGTa                       <  1:4532778/60‑1 (MQ=255)
                                                         tACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGAGAGCcagca                                     <  1:3503057/46‑1 (MQ=255)
                                                           cTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGc          <  1:4238413/71‑1 (MQ=255)
                                                           cTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGc          >  1:3597478/1‑71 (MQ=255)
                                                             gTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCg        >  1:4454455/1‑71 (MQ=255)
                                                              tttCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGcc         <  1:1339209/69‑1 (MQ=255)
                                                                tCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACa                          >  1:2412600/1‑50 (MQ=255)
                                                                tCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTAc     >  1:4392299/1‑71 (MQ=255)
                                                                tCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTAc     >  1:2381007/1‑71 (MQ=255)
                                                                  aGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACtt   <  1:1373143/71‑1 (MQ=255)
                                                                   gTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTc  <  1:2921677/71‑1 (MQ=255)
                                                                   |                                                                      
GACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTC  >  AP009048/861044‑861181

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