Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 862,250 G→T 13.8% T260N (ACT→AAT)  ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048862,2500GT13.8% 108.5 / 4.2 36T260N (ACT→AAT) ybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base G (15/16);  new base T (3/2);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGCCAGCAGTTCATCGCGACGGCTTT  >  AP009048/862181‑862319
                                                                     |                                                                     
gAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGt                                                                      >  1:895961/1‑71 (MQ=255)
 aGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTc                                                                     >  1:4142484/1‑71 (MQ=255)
  gTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCt                                                                    <  1:1869316/71‑1 (MQ=255)
  gTGACCGATAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCt                                                                    <  1:4570931/71‑1 (MQ=255)
   tGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTg                                                                   <  1:280741/71‑1 (MQ=255)
       cGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTgcggc                                                               >  1:3771604/1‑71 (MQ=255)
       cGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTgcg                                                                 >  1:1856071/1‑69 (MQ=255)
         ttAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGcggcgg                                                             >  1:2445919/1‑71 (MQ=255)
          tAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCggcggc                                                            >  1:2348638/1‑71 (MQ=255)
               tCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTc                                                                     >  1:4249415/1‑57 (MQ=255)
                 gATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTggg                                                     >  1:4090782/1‑71 (MQ=255)
                   tCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTggg                                                     <  1:884071/69‑1 (MQ=255)
                   tCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGc                                                    >  1:349244/1‑70 (MQ=255)
                    cTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGa                                                  >  1:4537404/1‑71 (MQ=255)
                       cAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAATTCTGCGGCGGCTGGTGGGCGATaa                                               <  1:1434073/71‑1 (MQ=255)
                       cAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGAt                                                 >  1:2750299/1‑69 (MQ=255)
                            tCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCg                                          <  1:1821618/71‑1 (MQ=255)
                                   gATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGtttt                                   <  1:3223090/71‑1 (MQ=255)
                                   gATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGtt                                     <  1:871611/69‑1 (MQ=255)
                                        aGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGcc                              <  1:928550/71‑1 (MQ=255)
                                        aGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGcc                              <  1:923228/71‑1 (MQ=255)
                                         gTAACACAGTTGCAGCGCCTGCCAGAAATTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCa                             <  1:2233826/71‑1 (MQ=255)
                                          tAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCg                                          >  1:1895048/1‑57 (MQ=255)
                                            acacAGTTGCAGCGCCTGCCAGAAATTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATc                           >  1:836392/1‑70 (MQ=255)
                                            acacAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTTGGCGATAAGATCGCAGTTTTCTGCCATc                           <  1:4284692/70‑1 (MQ=255)
                                                aGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGCcag                      <  1:1748054/71‑1 (MQ=255)
                                                           cTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGcc                        <  1:3130481/58‑1 (MQ=255)
                                                            tGCCAGAAATTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGCCAGCAGTTCATcgc           >  1:3218975/1‑70 (MQ=255)
                                                            tGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGcc                              <  1:2368725/51‑1 (MQ=255)
                                                             gacAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGc                         <  1:3192316/53‑1 (MQ=255)
                                                             gCCAGAAATTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGc                         >  1:1494127/1‑55 (MQ=255)
                                                               cAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGCCAGCAGtt                 >  1:1931230/1‑61 (MQ=255)
                                                               cAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGCCAGCAGTTCATCGCGACg       <  1:2182326/71‑1 (MQ=255)
                                                                  aaaGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGCCAGCAGTTCATCGCGACCGCt    >  1:1672775/1‑71 (MQ=255)
                                                                   aaGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGCCAGCAGTTCATCGCGACGGCtt   <  1:2678633/71‑1 (MQ=255)
                                                                    aGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGCCAGCAGTTCATCGCGACGGCttt  >  1:646183/1‑71 (MQ=255)
                                                                     |                                                                     
GAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGCCAGCAGTTCATCGCGACGGCTTT  >  AP009048/862181‑862319

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: