Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,099,547 C→T 10.5% A143V (GCC→GTC)  ycdX → predicted zinc‑binding hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,099,5470CT10.5% 110.5 / 4.0 38A143V (GCC→GTC) ycdXpredicted zinc‑binding hydrolase
Reads supporting (aligned to +/- strand):  ref base C (11/23);  new base T (2/2);  total (13/25)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

CAGCGGCAATGTGCATATAATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTC  >  AP009048/1099485‑1099613
                                                              |                                                                  
cAGCGGCAATGTGCATATAATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTg                                                            <  1:1677972/71‑1 (MQ=255)
  gCGGCAATGTGCATATAATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGaa                                                          <  1:1713751/71‑1 (MQ=255)
         tgtgCATATAATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAagc                                                        <  1:4291409/66‑1 (MQ=255)
         tgtgCATATAATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAag                                                         <  1:3654468/65‑1 (MQ=255)
             cATATAATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGa                                                <  1:663984/70‑1 (MQ=255)
               tataATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGaaa                                              >  1:4100049/1‑70 (MQ=255)
               tataATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGaa                                               <  1:3825647/69‑1 (MQ=255)
                 taataaGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACAt                                           <  1:1830529/71‑1 (MQ=255)
                  aataaGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAagcag                                                      >  1:341859/1‑59 (MQ=255)
                  aataaGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGa                                                <  1:2243899/65‑1 (MQ=255)
                   ataaGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTg                                                            <  1:1655994/52‑1 (MQ=255)
                   ataaGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTg                                                            <  1:3370315/52‑1 (MQ=255)
                   ataaGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATc                                          <  1:954108/70‑1 (MQ=255)
                         cATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGaa                                                          <  1:41939/48‑1 (MQ=255)
                          aTCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTg                                                            <  1:3345207/45‑1 (MQ=255)
                           tCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGgcgc                                 >  1:1521729/1‑71 (MQ=255)
                           tCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGgcgc                                 >  1:3575261/1‑71 (MQ=255)
                           tCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGAGAAACATCAGGTGgcgc                                 >  1:3839133/1‑71 (MQ=255)
                             ccGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAag                                                         >  1:2112208/1‑45 (MQ=255)
                             ccGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGa                                                <  1:1730120/54‑1 (MQ=255)
                               ggAAATCCCAAATATGAAATAGATGTGAAAGTCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTgg                              >  1:522168/1‑70 (MQ=255)
                               ggAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTgg                              >  1:4064040/1‑70 (MQ=255)
                               ggAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGa                             <  1:3498783/71‑1 (MQ=255)
                                      ccAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCaa                       >  1:3856590/1‑70 (MQ=255)
                                      ccAAATATGAAATAGATGAGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGgcgc                                 <  1:1825884/60‑1 (MQ=255)
                                       cAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGcc                                                    <  1:3248480/40‑1 (MQ=255)
                                        aaaTATGAAATAGATGTGAAAGTCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCaataa                    <  1:3408622/71‑1 (MQ=255)
                                         aaTATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATc                                          >  1:584715/1‑48 (MQ=255)
                                             tGAAATAGATGTGAAAGTCGTTGCTGAAGCAGCCGCGAAACATCAGGTg                                     >  1:330255/1‑49 (MQ=255)
                                              gAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCa                        <  1:2168050/61‑1 (MQ=255)
                                                     aTGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTa       <  1:2769908/71‑1 (MQ=255)
                                                     aTGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGCCGCTGGAAATCAATAATTCCTCAttt          <  1:177377/68‑1 (MQ=255)
                                                         gAAAGTCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATc                         <  1:2349857/49‑1 (MQ=255)
                                                          aaaGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCt                                >  1:1208376/1‑41 (MQ=255)
                                                          aaaGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTg                               <  1:2267179/42‑1 (MQ=255)
                                                          aaaGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCaataat                   >  1:2728942/1‑54 (MQ=255)
                                                          aaaGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTc  <  1:1169656/71‑1 (MQ=255)
                                                             gCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAAt                          <  1:944372/44‑1 (MQ=255)
                                                              |                                                                  
CAGCGGCAATGTGCATATAATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTC  >  AP009048/1099485‑1099613

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: