Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,201,015 G→C 12.2% T56T (ACG→ACC sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,201,0150GC12.2% 277.9 / 23.9 106T56T (ACG→ACCsdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base G (38/54);  new base C (9/4);  total (47/59)
Fisher's exact test for biased strand distribution p-value = 7.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 108 total aligned reads displayed.

TGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTG  >  Exported/1200945‑1201081
                                                                      |                                                                  
tGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCGCTGTCGCTGACg                                                                    >  1:5159125/1‑71 (MQ=255)
tGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACg                                                                    >  1:139326/1‑71 (MQ=255)
   tCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGt                                                                 <  1:1062648/71‑1 (MQ=255)
   tCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGt                                                                 <  1:2388242/71‑1 (MQ=255)
   tCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGt                                                                 <  1:138073/71‑1 (MQ=255)
    cGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGt                                                                 >  1:7122/1‑70 (MQ=255)
    cGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGt                                                                 >  1:1280185/1‑70 (MQ=255)
     gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTaa                                                               <  1:3658153/71‑1 (MQ=255)
     gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                                <  1:3705772/70‑1 (MQ=255)
     gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                                <  1:3107155/70‑1 (MQ=255)
      aaaaaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTaaa                                                              <  1:3949212/71‑1 (MQ=255)
        aaaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAgg                                                            >  1:1916355/1‑71 (MQ=255)
        aaaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAgg                                                            >  1:4885595/1‑71 (MQ=255)
        aaaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAgg                                                            >  1:3083624/1‑71 (MQ=255)
         aaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGc                                                           <  1:3214937/71‑1 (MQ=255)
            gCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                        <  1:4463620/71‑1 (MQ=255)
            gCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                        <  1:1785162/71‑1 (MQ=255)
            gCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                        <  1:4456703/71‑1 (MQ=255)
            gCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                        <  1:4158358/71‑1 (MQ=255)
                   ggATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACAcc                                                  >  1:4572350/1‑70 (MQ=255)
                    gATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACTGGTAAAGGCCACCACAcc                                                  <  1:5236837/69‑1 (MQ=255)
                    gATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACAcc                                                  <  1:3455469/69‑1 (MQ=255)
                    gATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                                >  1:5047711/1‑71 (MQ=255)
                     aTAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                                <  1:3828588/54‑1 (MQ=255)
                     aTAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                                <  1:2541638/54‑1 (MQ=255)
                     aTAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                                <  1:5176416/54‑1 (MQ=255)
                      agacGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCacac                                                   >  1:5109075/4‑66 (MQ=255)
                            tACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCAGACGGGTAAAGGCCACCACACCGATATCGCCa                                        >  1:460753/1‑71 (MQ=255)
                             aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGatat                                             >  1:780093/1‑65 (MQ=255)
                             aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGatat                                             >  1:4445484/1‑65 (MQ=255)
                             aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                       >  1:3298186/1‑71 (MQ=255)
                             aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                       >  1:4016157/1‑71 (MQ=255)
                             aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                       >  1:3678227/1‑71 (MQ=255)
                              cTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                       <  1:2910008/70‑1 (MQ=255)
                               tCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                               <  1:1827148/61‑1 (MQ=255)
                               tCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                               <  1:2844025/61‑1 (MQ=255)
                               tCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                               <  1:5104891/61‑1 (MQ=255)
                               tCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                               <  1:5310405/61‑1 (MQ=255)
                                 gcgTTGCCGTGGACGTTTATGGTTCACTTTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattat                                    >  1:2162360/1‑70 (MQ=255)
                                 gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCg                                           >  1:5673414/1‑63 (MQ=255)
                                 gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattatt                                   <  1:4326109/71‑1 (MQ=255)
                                 gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattat                                    <  1:2738794/70‑1 (MQ=255)
                                 gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattat                                    <  1:205515/70‑1 (MQ=255)
                                 gcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCattat                                    <  1:2320579/70‑1 (MQ=255)
                                   gTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCAttatta                                  <  1:3583519/70‑1 (MQ=255)
                                   gTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATtattat                                 <  1:3305787/71‑1 (MQ=255)
                                     tGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATgg                               >  1:1163418/1‑71 (MQ=255)
                                      gCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATggg                              >  1:1241847/1‑71 (MQ=255)
                                      gCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATggg                              >  1:1900878/1‑71 (MQ=255)
                                      gCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATgg                               <  1:4498710/70‑1 (MQ=255)
                                       ccGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCacac                                                   >  1:2770013/1‑49 (MQ=255)
                                         gttgACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCt                           >  1:557156/4‑71 (MQ=255)
                                          tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCtt                          <  1:4504614/71‑1 (MQ=255)
                                           ggACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATg                                >  1:1442507/1‑64 (MQ=255)
                                           ggACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATg                                >  1:435331/1‑64 (MQ=255)
                                            gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                         >  1:5487791/1‑70 (MQ=255)
                                            gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                         <  1:5210681/70‑1 (MQ=255)
                                            gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                        <  1:3776526/71‑1 (MQ=255)
                                            gACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACcaca                                                    >  1:2326386/1‑43 (MQ=12)
                                            gACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCAttatta                                  >  1:5683594/1‑61 (MQ=255)
                                            gACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCtgg                         <  1:2125531/70‑3 (MQ=255)
                                             aCGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                        <  1:87060/38‑1 (MQ=255)
                                             aCGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                        <  1:5262900/38‑1 (MQ=255)
                                             aCGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGc                        <  1:2226196/70‑1 (MQ=255)
                                              cGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg                      <  1:1742101/71‑1 (MQ=255)
                                              cGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg                      <  1:4697721/71‑1 (MQ=255)
                                               gTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                                <  1:347360/44‑1 (MQ=255)
                                               gTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                                <  1:3319133/44‑1 (MQ=255)
                                                tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAgg                     >  1:3736191/1‑70 (MQ=255)
                                                tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAgg                     >  1:3126734/1‑70 (MQ=255)
                                                tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTACTATg                                <  1:4664236/59‑1 (MQ=255)
                                                  tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                                >  1:374146/1‑41 (MQ=255)
                                                  tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGt                             >  1:5135927/1‑60 (MQ=255)
                                                  tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTaa                  >  1:5435198/1‑71 (MQ=255)
                                                  tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTa                   <  1:3053035/70‑1 (MQ=255)
                                                  tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTa                   <  1:2150604/70‑1 (MQ=255)
                                                       ttCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc             <  1:481076/71‑1 (MQ=255)
                                                         cACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt           >  1:1876963/1‑71 (MQ=255)
                                                          aCTGTCGCTGACGGGTAAAGTCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc            <  1:3951882/69‑1 (MQ=255)
                                                          aCTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc            <  1:1358070/69‑1 (MQ=255)
                                                          aCTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt           <  1:106621/70‑1 (MQ=255)
                                                          aCTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt           <  1:3041287/70‑1 (MQ=255)
                                                          aCTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTg          >  1:4793423/1‑71 (MQ=255)
                                                            tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc         >  1:5175328/1‑70 (MQ=255)
                                                            tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc         >  1:5539730/1‑70 (MQ=255)
                                                            tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc         >  1:2059994/1‑70 (MQ=255)
                                                            tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc         >  1:852152/1‑70 (MQ=255)
                                                             gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                        <  1:2980099/54‑1 (MQ=255)
                                                             gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                        <  1:3867789/54‑1 (MQ=255)
                                                             gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa               <  1:2673851/63‑1 (MQ=255)
                                                             gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc             >  1:2793122/1‑65 (MQ=255)
                                                              tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg                      <  1:1518438/55‑1 (MQ=255)
                                                              tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc      >  1:1145933/1‑71 (MQ=255)
                                                                tcTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATg                                >  1:2541223/2‑43 (MQ=38)
                                                                tcTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATg                                >  1:3872182/2‑43 (MQ=38)
                                                                tcTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATg                                >  1:4313237/2‑43 (MQ=38)
                                                                  tGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTg  >  1:2724003/1‑71 (MQ=255)
                                                                  tGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAAc                 >  1:24917/1‑56 (MQ=255)
                                                                   gACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTg  >  1:2265245/1‑70 (MQ=255)
                                                                    aCGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCAACCGTg  <  1:2093796/69‑1 (MQ=255)
                                                                      |                                                                  
TGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTG  >  Exported/1200945‑1201081

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: