Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,201,067 G→C 13.8% E74Q (GAA→CAA) ‡ sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,201,0670GC13.8% 312.4 / 32.9 123E74Q (GAA→CAA) ‡sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base G (52/54);  new base C (6/11);  total (58/65)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

Only 100 of 128 total aligned reads displayed.

GTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTG  >  Exported/1200999‑1201135
                                                                    |                                                                    
gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                                     <  1:912467/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:2059994/1‑70 (MQ=255)
      tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                               >  1:5539730/1‑70 (MQ=255)
        tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                            >  1:1145933/1‑71 (MQ=255)
            tGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTg                                                        >  1:2724003/1‑71 (MQ=255)
             gACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTg                                                        >  1:2265245/1‑70 (MQ=255)
              aCGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCAACCGTg                                                        <  1:2093796/69‑1 (MQ=255)
                 ggTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                                                   >  1:5640380/1‑71 (MQ=255)
                  gTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGatat                                                   >  1:16107/1‑70 (MQ=255)
                   tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                                      <  1:3137000/50‑1 (MQ=255)
                   tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                                      <  1:2963199/50‑1 (MQ=255)
                   tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                                      <  1:4425788/50‑1 (MQ=255)
                   tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                                      <  1:2645510/50‑1 (MQ=255)
                    aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGata                                                    >  1:2191146/1‑67 (MQ=255)
                     aaGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGa                                                >  1:797924/1‑70 (MQ=255)
                     aaGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGa                                                >  1:4609480/1‑70 (MQ=255)
                     aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                                <  1:3802516/70‑1 (MQ=255)
                     aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                                <  1:5313696/70‑1 (MQ=255)
                      aGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGatt                                              <  1:1721752/71‑3 (MQ=255)
                       ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                             <  1:3683581/71‑4 (MQ=255)
                       ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                             <  1:822517/71‑4 (MQ=255)
                       ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                             <  1:2724485/71‑4 (MQ=255)
                       ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                                          >  1:3985719/1‑58 (MQ=255)
                        gCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                             >  1:3360804/1‑67 (MQ=255)
                         ccaccaCACCGATATCGCCATTTTTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTa                                           >  1:5507829/1‑71 (MQ=255)
                           accacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                                     <  1:520764/43‑1 (MQ=255)
                            ccacACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                                                  <  1:1903575/61‑1 (MQ=255)
                            ccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc                                                                  <  1:5312026/45‑1 (MQ=255)
                            ccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTg                                                                <  1:1413433/47‑1 (MQ=255)
                             cacaCCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                                   >  1:4101675/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)
                               cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg                                     >  1:4176217/1‑71 (MQ=255)
                                 ccGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGt                                   >  1:1125878/1‑71 (MQ=255)
                                 ccGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGt                                   >  1:3155375/1‑71 (MQ=255)
                                  tgATATCGCCATTATTATGGGCCTTGCGGGTAACCTGCCGGATACCGTGGATATCGat                                               <  1:2910886/57‑2 (MQ=255)
                                    atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGAc                                               <  1:5570109/56‑1 (MQ=255)
                                    atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGttt                                 >  1:4485513/1‑70 (MQ=255)
                                     tatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTa                               <  1:5695205/71‑1 (MQ=255)
                                     tatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTa                               <  1:2652223/71‑1 (MQ=255)
                                      atCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                              >  1:2619947/1‑71 (MQ=255)
                                       tCGCCATTATTATGGGTCTTGCAGGTAAGGAACCTGCCACCGTGGatat                                                   >  1:8664/1‑49 (MQ=255)
                                       tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                                          >  1:799848/1‑42 (MQ=255)
                                       tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                                          >  1:181865/1‑42 (MQ=255)
                                       tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                                                   >  1:3658793/1‑49 (MQ=255)
                                       tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAtt                             >  1:4179570/1‑71 (MQ=255)
                                       tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                              <  1:1733972/70‑1 (MQ=255)
                                       tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                              <  1:5637959/70‑1 (MQ=255)
                                       tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                              <  1:349867/70‑1 (MQ=255)
                                       tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                              <  1:5274320/70‑1 (MQ=255)
                                       tCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTAt                              >  1:2231524/1‑70 (MQ=255)
                                        cGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTa                                           >  1:42484/1‑56 (MQ=255)
                                          ccATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGt                                                         <  1:334685/40‑1 (MQ=255)
                                          ccATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGt                                                         <  1:1404212/40‑1 (MQ=255)
                                          ccATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgc                          <  1:438116/71‑1 (MQ=255)
                                           cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                             <  1:3015882/51‑1 (MQ=255)
                                           cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                             <  1:166450/51‑1 (MQ=255)
                                           cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                             <  1:3332616/51‑1 (MQ=255)
                                           cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                             <  1:4998195/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)
                                              tattatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACg                      >  1:5009563/1‑71 (MQ=255)
                                              tattatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACg                      >  1:5206861/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:5295182/70‑1 (MQ=255)
                                                  atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag                   <  1:1330065/70‑1 (MQ=255)
                                                   tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaa                 <  1:2686352/71‑1 (MQ=255)
                                                     ggTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaaga               <  1:5659484/71‑1 (MQ=255)
                                                      gTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAagag              <  1:4348460/71‑1 (MQ=255)
                                                       tCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgc             <  1:1277218/71‑1 (MQ=255)
                                                       tCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgc             >  1:4706057/1‑71 (MQ=255)
                                                       tCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgc             <  1:5406685/71‑1 (MQ=255)
                                                           gCAGGTAACGAACCTGCGACCGTGGATATCGACAGTATTCCCGGTTTTATTc                            >  1:2844949/1‑52 (MQ=255)
                                                           gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAtt                             >  1:4062546/1‑51 (MQ=255)
                                                           gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTc                            >  1:3193066/1‑52 (MQ=255)
                                                           gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg          <  1:3789326/70‑1 (MQ=255)
                                                           gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg          <  1:3664178/70‑1 (MQ=255)
                                                               gTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGt     <  1:1867898/71‑1 (MQ=255)
                                                               gTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGt     <  1:83741/71‑1 (MQ=255)
                                                                tAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg          <  1:2144097/65‑1 (MQ=255)
                                                                 aaCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGtt                                  >  1:2542514/1‑40 (MQ=255)
                                                                 aaCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc           >  1:5675277/1‑63 (MQ=255)
                                                                  aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGttt                                 <  1:608601/40‑1 (MQ=255)
                                                                  aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAAc       >  1:2893734/1‑66 (MQ=255)
                                                                  aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAAc       >  1:1699940/1‑66 (MQ=255)
                                                                  aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctg  >  1:5413966/1‑71 (MQ=255)
                                                                  aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctg  <  1:1717414/71‑1 (MQ=255)
                                                                   cGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaa                 >  1:543348/1‑55 (MQ=255)
                                                                    |                                                                    
GTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTG  >  Exported/1200999‑1201135

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: