Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,485,132 G→A 10.9% S162N (AGC→AAC)  dsdA → D‑serine ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,485,1320GA10.9% 131.1 / 6.4 46S162N (AGC→AAC) dsdAD‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base G (21/20);  new base A (2/3);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

TTGCTGACGCTTGATGATGACTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAGCGCCCGCATTGGCTTTAAG  >  W3110S.gb/2485062‑2485202
                                                                      |                                                                      
ttGCTGACGCTTGATGATGACTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAg                                                                        <  1:987611/71‑1 (MQ=255)
ttGCTGACGCTTGATGATGACTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAg                                                                        >  1:3151985/1‑71 (MQ=255)
ttGCTGACGCTTGATGATGACTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAg                                                                        >  1:2665428/1‑71 (MQ=255)
ttGCTGACGCTTGATGATGACTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAg                                                                        >  1:3952364/1‑71 (MQ=255)
ttGCTGACGCTTGATGATGACTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAg                                                                        >  1:913349/1‑71 (MQ=255)
ttGCTGACGCTTGATGATGACTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAg                                                                        >  1:1655860/1‑71 (MQ=255)
  gCTGACGCTTGATGATGACTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCa                                                                      >  1:766106/1‑71 (MQ=255)
      aCGCTTGATGATGACTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTGc                                                                  <  1:2225095/71‑1 (MQ=255)
        gCTTGATGATGACTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCt                                                                 <  1:1801639/70‑1 (MQ=255)
          ttGATGATGACTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCtgtg                                                              >  1:3114615/1‑71 (MQ=255)
           tgatgatgaCTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCa                                                                      >  1:862725/1‑62 (MQ=255)
             atgatgaCTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAg                                                                        <  1:2033114/58‑1 (MQ=255)
             atgatgaCTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAg                                                                        <  1:401862/58‑1 (MQ=255)
              tgatgaCTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAACATTGCTGTGGGCt                                                          >  1:4196449/1‑71 (MQ=255)
                atgaCTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCAt                                                                     <  1:379603/58‑1 (MQ=255)
                  gaCTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTg                                                                   >  1:614353/1‑58 (MQ=255)
                   aCTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAAc                                                      <  1:1920243/70‑1 (MQ=255)
                       cAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTg                                                                   >  1:1965635/1‑53 (MQ=255)
                          cAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATc                                              >  1:275045/1‑71 (MQ=255)
                            aaCTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTg                                            <  1:227018/71‑1 (MQ=255)
                            aaCTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTg                                            >  1:1469735/1‑71 (MQ=255)
                              cTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAg                                                                        >  1:1152068/1‑41 (MQ=255)
                              cTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTggg                                          >  1:3762314/1‑71 (MQ=255)
                                gCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGt                                        >  1:4177902/1‑71 (MQ=255)
                                   ttCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCa                                                                      <  1:942827/38‑1 (MQ=255)
                                   ttCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCa                                                                      <  1:3643112/38‑1 (MQ=255)
                                     ctcCGGAGTTTAAACAGTTCTTTAGCCAATACAGCAtt                                                                    >  1:4433835/1‑38 (MQ=255)
                                      tcCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCaa                                  >  1:3189336/1‑71 (MQ=255)
                                         ggAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATc                                              >  1:2789368/1‑56 (MQ=255)
                                         ggAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCg                               >  1:2820167/1‑71 (MQ=255)
                                          gAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAAt                                 >  1:3607524/1‑68 (MQ=255)
                                             tttAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAAcc                                                     <  1:3521514/45‑1 (MQ=255)
                                             tttAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAAcc                                                     <  1:4257098/45‑1 (MQ=255)
                                             tttAAACAGTTCTTTAGCCAATACAACATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCAt                           <  1:304573/71‑1 (MQ=255)
                                               tAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCATCCGGAAATCTGGGGTTATCAATCGGCAtt                          <  1:3533510/70‑1 (MQ=255)
                                                aaaCAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGTTTATCAATCGGCATTa                         <  1:2388262/70‑1 (MQ=255)
                                                 aaCAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTAt                        >  1:100669/1‑70 (MQ=255)
                                                 aaCAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATg                       <  1:3107404/71‑1 (MQ=255)
                                                     gTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATc                                    <  1:3112138/54‑1 (MQ=255)
                                                      ttCTTTAGCCAATACAACATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAgc                    <  1:3345663/69‑1 (MQ=255)
                                                       tCTTTAGCCAATACAACATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAgcgc                  <  1:3492752/70‑1 (MQ=255)
                                                        cTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAGCGccc                >  1:3570480/1‑71 (MQ=255)
                                                        cTTTAGCCAATACAACATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGa                      >  1:3491810/1‑65 (MQ=255)
                                                           tAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAGCGCCCg               <  1:3368862/69‑1 (MQ=255)
                                                           tAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAGCGCCCGc              <  1:3283479/70‑1 (MQ=255)
                                                            aGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAGCGCCCGCa             <  1:2360754/70‑1 (MQ=255)
                                                                aaTACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATc                                    >  1:4368998/1‑43 (MQ=255)
                                                                  ttcAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTa                         <  1:1923213/50‑1 (MQ=255)
                                                                   aCAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCa                                   <  1:2973821/41‑1 (MQ=255)
                                                                   aCAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCAtt                          >  1:2521557/1‑50 (MQ=255)
                                                                    cAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAGCGCCCGCATTGGCttt     >  1:3548403/1‑70 (MQ=255)
                                                                    cAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAGCGCCCGCATTGGCttt     >  1:2591785/1‑70 (MQ=255)
                                                                    cAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAGCGCCCGCATTGGCttt     <  1:1133098/70‑1 (MQ=255)
                                                                    cAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAGCGCCCGCATTGGCTTTa    <  1:3150159/71‑1 (MQ=255)
                                                                     aGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCaa                                  >  1:671633/1‑40 (MQ=255)
                                                                      gCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAGCGCCCGCATTGGCTTTAAg  >  1:3808677/1‑71 (MQ=255)
                                                                      |                                                                      
TTGCTGACGCTTGATGATGACTACAGCAAACTGCTTTCTCCGGAGTTTAAACAGTTCTTTAGCCAATACAGCATTGCTGTGGGCTCAACCGGAAATCTGGGGTTATCAATCGGCATTATGAGCGCCCGCATTGGCTTTAAG  >  W3110S.gb/2485062‑2485202

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: