Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,202,086 C→T 15.5% R413R (CGC→CGT)  sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,202,0860CT15.5% 149.5 / 15.0 58R413R (CGC→CGT) sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base C (28/21);  new base T (3/6);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

GTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACG  >  Exported/1202020‑1202156
                                                                  |                                                                      
gTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCAc                                                                      >  1:762809/1‑69 (MQ=255)
gTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCAc                                                                      >  1:1288866/1‑69 (MQ=255)
gTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCAc                                                                      >  1:5285196/1‑69 (MQ=255)
    aTTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGc                                                                <  1:4469686/71‑1 (MQ=255)
       gAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCa                                                               >  1:1055199/1‑69 (MQ=255)
       gAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCa                                                               >  1:5573887/1‑69 (MQ=255)
       gAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCAcc                                                             >  1:580349/1‑71 (MQ=255)
           gTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACcgcgc                                                         <  1:58668/71‑1 (MQ=255)
           gTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACcgcgc                                                         <  1:5164087/71‑1 (MQ=255)
            tAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACcgcgcg                                                        >  1:880853/1‑71 (MQ=255)
                gCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtctc                                                    >  1:3333967/1‑71 (MQ=255)
                gCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtctc                                                    >  1:3739672/1‑71 (MQ=255)
                 ccATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtctc                                                    >  1:4229551/1‑70 (MQ=255)
                 ccATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCg                                                   >  1:2211402/1‑71 (MQ=255)
                   aTTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCg                                                   <  1:845554/69‑1 (MQ=255)
                     tGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTg                                                >  1:611908/1‑70 (MQ=255)
                     tGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTg                                                >  1:3480006/1‑70 (MQ=255)
                          ctGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCa                                                               >  1:4869360/1‑50 (MQ=255)
                          ctGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAg                                          >  1:4865052/1‑71 (MQ=255)
                          ctGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAg                                          >  1:4823271/1‑71 (MQ=255)
                           tgtgAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAgg                                         >  1:1298988/1‑71 (MQ=255)
                           tgtgAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAgg                                         >  1:4900608/1‑71 (MQ=255)
                           tgtgAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAgg                                         >  1:3969631/1‑71 (MQ=255)
                             tgAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTc                                       <  1:3333494/71‑1 (MQ=255)
                               aaGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCa                                      >  1:1598500/1‑70 (MQ=255)
                               aaGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCAt                                     <  1:1483214/71‑1 (MQ=255)
                                aGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATc                                    <  1:1486707/71‑1 (MQ=255)
                                   cGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaa                                 <  1:1944295/71‑1 (MQ=255)
                                   cGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaa                                 <  1:1098307/71‑1 (MQ=255)
                                      ttAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAAc                               >  1:3092243/1‑70 (MQ=255)
                                       tAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                             >  1:573100/1‑71 (MQ=255)
                                        aaCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                             <  1:2635090/70‑1 (MQ=255)
                                        aaCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                            >  1:4061821/1‑71 (MQ=255)
                                        aaCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAAcca                             <  1:2331851/70‑3 (MQ=9)
                                        aaCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAAcca                             >  1:1675178/1‑68 (MQ=9)
                                        aaCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAAcca                             <  1:2338470/70‑3 (MQ=9)
                                         aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                            <  1:1416027/70‑1 (MQ=255)
                                         aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                            <  1:5636342/70‑1 (MQ=255)
                                         aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                           >  1:1591409/1‑71 (MQ=255)
                                         aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                           <  1:214850/71‑1 (MQ=255)
                                          cgccgcGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAAc                               <  1:1318944/66‑1 (MQ=255)
                                            ccgcGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                             >  1:2896731/1‑66 (MQ=255)
                                             cgcacgTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACg                       <  1:2094673/71‑1 (MQ=12)
                                              gcgcgGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGa                      <  1:3272937/71‑1 (MQ=255)
                                                 cgtaTGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAAc                   <  1:3611663/68‑1 (MQ=12)
                                                   gATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCg                 >  1:3287383/1‑71 (MQ=255)
                                                      ggCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGt               >  1:1703545/1‑70 (MQ=255)
                                                       gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTa              <  1:374468/70‑1 (MQ=255)
                                                       gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTa              <  1:5028452/70‑1 (MQ=255)
                                                       gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTa              <  1:602788/70‑1 (MQ=255)
                                                        ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCg                                   <  1:1429092/48‑1 (MQ=255)
                                                         gctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaaag           <  1:3752269/70‑3 (MQ=12)
                                                         gctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaaag           <  1:155195/70‑3 (MQ=12)
                                                          ctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaa             >  1:4658923/1‑68 (MQ=16)
                                                          ctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaa             >  1:1338565/1‑68 (MQ=16)
                                                          ctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa          >  1:5130055/1‑71 (MQ=255)
                                                           tGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTAc                        <  1:279462/56‑1 (MQ=255)
                                                           tGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTAc                        <  1:5138410/56‑1 (MQ=255)
                                                             cgccgtACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGata        <  1:1741280/70‑3 (MQ=9)
                                                             cgccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGa                      >  1:4681059/1‑56 (MQ=255)
                                                             cgccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt       <  1:5233609/71‑1 (MQ=255)
                                                              gccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg      <  1:3386637/71‑1 (MQ=255)
                                                              gccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg      <  1:2945933/71‑1 (MQ=255)
                                                              gccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg      <  1:1053086/71‑1 (MQ=255)
                                                               ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa     <  1:485031/71‑1 (MQ=255)
                                                               ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa     <  1:2237289/71‑1 (MQ=255)
                                                               ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa     <  1:740483/71‑1 (MQ=255)
                                                               ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa     <  1:3133687/71‑1 (MQ=255)
                                                                cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaaa                                >  1:1721249/1‑43 (MQ=255)
                                                                cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaaa                                >  1:4772866/1‑43 (MQ=255)
                                                                cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg      <  1:122387/69‑1 (MQ=255)
                                                                cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGaa    >  1:2482802/1‑71 (MQ=255)
                                                                  cACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACg  >  1:2128450/1‑71 (MQ=255)
                                                                  cACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACg  >  1:3152718/1‑71 (MQ=255)
                                                                  |                                                                      
GTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACG  >  Exported/1202020‑1202156

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: ‑