Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,201,054 T→C 17.2% G69G (GGT→GGC sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,201,0540TC17.2% 279.3 / 46.2 129G69G (GGT→GGCsdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base T (48/58);  new base C (8/14);  total (57/72)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 133 total aligned reads displayed.

TGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGC  >  Exported/1200987‑1201124
                                                                   |                                                                      
tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACCCCTATATCGCCATTATTATGGGTCtt                                                                     >  1:4107399/1‑71 (MQ=255)
tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCtt                                                                     <  1:4504614/71‑1 (MQ=255)
  gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                                                                    >  1:5487791/1‑70 (MQ=255)
  gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                                                                    <  1:5210681/70‑1 (MQ=255)
  gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                                                                   <  1:3776526/71‑1 (MQ=255)
  gACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCtgg                                                                    <  1:2125531/70‑3 (MQ=255)
   aCGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGc                                                                   <  1:2226196/70‑1 (MQ=255)
    cGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg                                                                 <  1:4697721/71‑1 (MQ=255)
    cGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg                                                                 <  1:1742101/71‑1 (MQ=255)
      tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAgg                                                                >  1:3126734/1‑70 (MQ=255)
      tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAgg                                                                >  1:3736191/1‑70 (MQ=255)
        tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGt                                                                        >  1:5135927/1‑60 (MQ=255)
        tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTaa                                                             >  1:5435198/1‑71 (MQ=255)
            gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                          <  1:912467/70‑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:2059994/1‑70 (MQ=255)
                  tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                    >  1:5175328/1‑70 (MQ=255)
                  tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                    >  1:5539730/1‑70 (MQ=255)
                   gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                                                                   <  1:2980099/54‑1 (MQ=255)
                   gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                                                                   <  1:2941518/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)
                        tGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAAc                                                            >  1:24917/1‑56 (MQ=255)
                          aCGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCAACCGTg                                             <  1:2093796/69‑1 (MQ=255)
                             ggTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                                        >  1:5640380/1‑71 (MQ=255)
                              gTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGatat                                        >  1:16107/1‑70 (MQ=255)
                              gTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                                       <  1:3545152/71‑1 (MQ=255)
                               tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTaa                                                             >  1:1116551/1‑48 (MQ=255)
                               tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                           <  1:568434/50‑1 (MQ=255)
                               tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                           <  1:3137000/50‑1 (MQ=255)
                               tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                           <  1:4425788/50‑1 (MQ=255)
                                aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGata                                         >  1:2191146/1‑67 (MQ=255)
                                 aaGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGa                                     >  1:797924/1‑70 (MQ=255)
                                 aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                     <  1:5313696/70‑1 (MQ=255)
                                 aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                     <  1:3802516/70‑1 (MQ=255)
                                   ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                  <  1:2724485/71‑4 (MQ=255)
                                   ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                  <  1:822517/71‑4 (MQ=255)
                                   ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                               >  1:3985719/1‑58 (MQ=255)
                                    gCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                  >  1:3360804/1‑67 (MQ=255)
                                    gCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGa                                           >  1:1748441/1‑61 (MQ=255)
                                       accacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                          <  1:2920708/43‑1 (MQ=255)
                                        ccacACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                                       <  1:1903575/61‑1 (MQ=255)
                                        ccacACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                  <  1:5044630/66‑4 (MQ=255)
                                        ccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc                                                       <  1:5312026/45‑1 (MQ=255)
                                        ccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc                                                       <  1:4402897/45‑1 (MQ=255)
                                         cacaCCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                        >  1:1883555/1‑43 (MQ=255)
                                          acacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGa                                           >  1:1579627/1‑55 (MQ=255)
                                          acacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTccc                           >  1:2749234/1‑71 (MQ=255)
                                           cacCGATATCGCCATTATTATGGTTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg                          >  1:688421/1‑71 (MQ=255)
                                           cacCGATATCGCCATTATTATGGGTCTTGCAGGTaa                                                             >  1:523184/1‑36 (MQ=255)
                                           cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg                          >  1:1683737/1‑71 (MQ=255)
                                           cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg                          >  1:4176217/1‑71 (MQ=255)
                                            acCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCgg                         >  1:5515268/1‑71 (MQ=255)
                                            acCGATATCGCCATTATTATGGGACTTCCACGTaa                                                             >  1:4293903/1‑35 (MQ=21)
                                             ccGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGt                        >  1:1125878/1‑71 (MQ=255)
                                                atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGAc                                    <  1:5570109/56‑1 (MQ=255)
                                                atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGAc                                    <  1:2376697/56‑1 (MQ=255)
                                                 tatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTa                    <  1:2652223/71‑1 (MQ=255)
                                                 tatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTa                    <  1:5695205/71‑1 (MQ=255)
                                                 tatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTa                    <  1:14276/71‑1 (MQ=255)
                                                  atCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                 >  1:5268274/1‑41 (MQ=255)
                                                  atCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                   >  1:2619947/1‑71 (MQ=255)
                                                   tCGCCATTATTATGGGTCTTGCAGGTAAGGAACCTGCCACCGTGGatat                                        >  1:8664/1‑49 (MQ=255)
                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                               >  1:799848/1‑42 (MQ=255)
                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                                        >  1:3658793/1‑49 (MQ=255)
                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAtt                  >  1:4179570/1‑71 (MQ=255)
                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                   <  1:349867/70‑1 (MQ=255)
                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                   <  1:1733972/70‑1 (MQ=255)
                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                   <  1:5637959/70‑1 (MQ=255)
                                                   tCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTAt                   >  1:2231524/1‑70 (MQ=255)
                                                    cGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTa                                >  1:42484/1‑56 (MQ=255)
                                                      ccATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGt                                              <  1:1404212/40‑1 (MQ=255)
                                                      ccATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGt                                              <  1:334685/40‑1 (MQ=255)
                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:4998195/51‑1 (MQ=255)
                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:3015882/51‑1 (MQ=255)
                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:166450/51‑1 (MQ=255)
                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:3332616/51‑1 (MQ=255)
                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgcg              <  1:4439293/71‑1 (MQ=255)
                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgc               >  1:458109/1‑70 (MQ=255)
                                                        attattATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:5573431/50‑1 (MQ=255)
                                                        attattATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:4008774/50‑1 (MQ=255)
                                                        attattATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGa             >  1:3070030/1‑71 (MQ=255)
                                                          tattatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACg           >  1:5206861/1‑71 (MQ=255)
                                                          tattatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACg           >  1:5009563/1‑71 (MQ=255)
                                                          tattatGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTAt                   <  1:918092/63‑1 (MQ=255)
                                                           attatGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTc                 >  1:2963194/1‑64 (MQ=255)
                                                            ttatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTa         <  1:4607497/71‑1 (MQ=255)
                                                              atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag        <  1:1330065/70‑1 (MQ=255)
                                                              atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag        <  1:5295182/70‑1 (MQ=255)
                                                               tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaa      <  1:2686352/71‑1 (MQ=255)
                                                                  gTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAagag   <  1:4348460/71‑1 (MQ=255)
                                                                   tCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgc  <  1:5406685/71‑1 (MQ=255)
                                                                   tCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgc  >  1:4706057/1‑71 (MQ=255)
                                                                   |                                                                      
TGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGC  >  Exported/1200987‑1201124

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: