Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,201,075 C→T 16.1% A76A (GCC→GCT) ‡ sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,201,0750CT16.1% 240.7 / 33.6 106A76A (GCC→GCT) ‡sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base C (44/43);  new base T (6/11);  total (51/55)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 112 total aligned reads displayed.

TCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGCAC  >  Exported/1201007‑1201145
                                                                    |                                                                      
tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                                      >  1:1145933/1‑71 (MQ=255)
    tGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTg                                                                  >  1:2724003/1‑71 (MQ=255)
      aCGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCAACCGTg                                                                  <  1:2093796/69‑1 (MQ=255)
         ggTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                                                             >  1:5640380/1‑71 (MQ=255)
          gTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                                                            <  1:3545152/71‑1 (MQ=255)
            aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGata                                                              >  1:2191146/1‑67 (MQ=255)
             aaGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGa                                                          >  1:4609480/1‑70 (MQ=255)
             aaGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGa                                                          >  1:797924/1‑70 (MQ=255)
             aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                                          <  1:5313696/70‑1 (MQ=255)
             aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                                          <  1:3802516/70‑1 (MQ=255)
              aGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGatt                                                        <  1:1721752/71‑3 (MQ=255)
               ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                                       <  1:2724485/71‑4 (MQ=255)
               ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                                       <  1:822517/71‑4 (MQ=255)
               ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                                       <  1:3683581/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)
                 ccaccaCACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGatt                                                        <  1:5359276/68‑3 (MQ=255)
                 ccaccaCACCGATATCGCCATTTTTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTa                                                     >  1:5507829/1‑71 (MQ=255)
                    ccacACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                                                            <  1:1903575/61‑1 (MQ=255)
                    ccacACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                                       <  1:5044630/66‑4 (MQ=255)
                     cacaCTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                                                            <  1:4232881/60‑1 (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)
                       cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg                                               >  1:1683737/1‑71 (MQ=255)
                        acCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCgg                                              >  1:5515268/1‑71 (MQ=255)
                         ccGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGt                                             >  1:3155375/1‑71 (MQ=255)
                         ccGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGt                                             >  1:1125878/1‑71 (MQ=255)
                          tgATATCGCCATTATTATGGGCCTTGCGGGTAACCTGCCGGATACCGTGGATATCGat                                                         <  1:2910886/57‑2 (MQ=255)
                            atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGAc                                                         <  1:5570109/56‑1 (MQ=255)
                            atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGAc                                                         <  1:2376697/56‑1 (MQ=255)
                            atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGttt                                           >  1:4485513/1‑70 (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)
                               tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                                                    >  1:181865/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:5637959/70‑1 (MQ=255)
                               tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                                        <  1:1733972/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:166450/51‑1 (MQ=255)
                                   cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                       <  1:4998195/51‑1 (MQ=255)
                                   cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                       <  1:3332616/51‑1 (MQ=255)
                                   cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                       <  1:3015882/51‑1 (MQ=255)
                                   cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                       <  1:712724/51‑1 (MQ=255)
                                   cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgcg                                   <  1:4439293/71‑1 (MQ=255)
                                   cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgc                                    >  1:458109/1‑70 (MQ=255)
                                    attattATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                       <  1:4008774/50‑1 (MQ=255)
                                    attattATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGa                                  >  1:3070030/1‑71 (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:1330065/70‑1 (MQ=255)
                                          atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag                             <  1:5295182/70‑1 (MQ=255)
                                           tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaa                           <  1:2686352/71‑1 (MQ=255)
                                             ggTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaaga                         <  1:5659484/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)
                                                   gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAtt                                       >  1:4062546/1‑51 (MQ=255)
                                                   gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTc                                      >  1:1314832/1‑52 (MQ=255)
                                                   gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTc                                      >  1:3193066/1‑52 (MQ=255)
                                                   gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg                    <  1:3789326/70‑1 (MQ=255)
                                                   gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg                    <  1:3664178/70‑1 (MQ=255)
                                                      ggTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACg                <  1:4591619/71‑1 (MQ=255)
                                                       gTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGt               <  1:83741/71‑1 (MQ=255)
                                                       gTAACGAACATGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGt               >  1:4708744/1‑71 (MQ=255)
                                                        tAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg                    <  1:2144097/65‑1 (MQ=255)
                                                         aaCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGtt                                            >  1:2542514/1‑40 (MQ=255)
                                                         aaCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc                     >  1:5675277/1‑63 (MQ=255)
                                                         aaCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc                     >  1:2158501/1‑63 (MQ=255)
                                                          aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGttt                                           <  1:5580462/40‑1 (MQ=255)
                                                          aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGttt                                           <  1:608601/40‑1 (MQ=255)
                                                          aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAAc                 >  1:1699940/1‑66 (MQ=255)
                                                          aCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAAc                 >  1:2893734/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)
                                                             aaCCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgctg         >  1:4854142/1‑71 (MQ=255)
                                                             aaCCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgatg         >  1:1298047/1‑71 (MQ=255)
                                                               ccTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgctgct       <  1:2169312/71‑1 (MQ=255)
                                                                cTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgc                       <  1:2971960/54‑1 (MQ=255)
                                                                   ccACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGca   <  1:2405849/71‑1 (MQ=255)
                                                                    cACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGa                                  >  1:5690424/1‑39 (MQ=255)
                                                                    cACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGa                                  >  1:1940507/1‑39 (MQ=255)
                                                                    cACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGcac  <  1:4681650/71‑1 (MQ=255)
                                                                    |                                                                      
TCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGCAC  >  Exported/1201007‑1201145

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: