Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,117,648 A→G F144S (TTC→TCC)  SAUSA300_RS05490 ← glycerophosphoryl diester phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,117,6480AG98.3% 153.9 / ‑6.2 60F144S (TTC→TCC) SAUSA300_RS05490glycerophosphoryl diester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (2/57);  total (2/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATTTGTGGTGCGATTGAGCCTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGAAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTGATATCTTTGAAGTGATAAGCTGCGTCTAATTTCTT  >  NC_007793/1117554‑1117733
                                                                                              |                                                                                     
aaTTTGTGGTGCGATTGAGCCTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGtt                                                                                  <  1:3013036/100‑1 (MQ=255)
     gtggtgCGATTGAGCCTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGtat                                                                             >  1:4639447/1‑100 (MQ=255)
              tccgatcTTCGTAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCtgtgtg                                                                    <  1:3874107/94‑1 (MQ=255)
              ctGAGCCTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCtgtgtg                                                                    <  1:5113334/99‑1 (MQ=255)
               tGAGGCTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCGGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCtgtgtgt                                                                   <  1:2571757/100‑1 (MQ=255)
                 ggCCTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGc                                                                 <  1:1029499/99‑1 (MQ=255)
                 gatcTTTATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGc                                                                 <  1:1247434/97‑1 (MQ=255)
                 gatcTTCATGGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGc                                                                 <  1:4018824/97‑1 (MQ=255)
                 ctccTTCGTAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGc                                                                 <  1:1734592/98‑1 (MQ=255)
                 aGCCTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGc                                                                 <  1:195368/100‑1 (MQ=255)
                 aGCCTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGAAGGTTAGTATAGCTGTGTGTGc                                                                 <  1:5731058/100‑1 (MQ=255)
                  gCCTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCa                                                                <  1:1672522/100‑1 (MQ=255)
                   gcctCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCag                                                               <  1:1031343/97‑2 (MQ=255)
                   ccTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCtgtgtgtggaa                                                               <  1:2639293/100‑4 (MQ=255)
                   ccTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCAt                                                               <  1:3073526/100‑1 (MQ=255)
                    ttccATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATg                                                              <  1:1984763/97‑1 (MQ=255)
                    cTTCATAGGATTCAGGAGCATCTTTTAGGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATg                                                              <  1:5696273/100‑1 (MQ=255)
                    cTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCaga                                                              <  1:217874/100‑3 (MQ=255)
                    cTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATg                                                              <  1:5602925/100‑1 (MQ=255)
                    cTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATAGCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATg                                                              <  1:2484347/100‑1 (MQ=255)
                    cTTCAGAGGAAGAAGGAGCATCTGTTAAGTCAACAGGTATAGACATATCAGGATATAGCATAAGCAACTCATCGGAGGTTAGTATAGCTGTGTGGGCATg                                                              <  1:611590/100‑1 (MQ=255)
                     ttCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCAtac                                                             <  1:2482627/100‑3 (MQ=255)
                     ttCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATgt                                                             <  1:302761/100‑2 (MQ=255)
                     ttCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATgc                                                             <  1:3318600/100‑2 (MQ=255)
                       tgtAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGAcc                                                           >  1:2440545/3‑100 (MQ=255)
                       cATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGAcc                                                           <  1:5734300/100‑1 (MQ=255)
                        aTAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCa                                                          <  1:4470922/100‑1 (MQ=255)
                          aGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACg                                                        <  1:2351390/100‑1 (MQ=255)
                            gATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCaggga                                                      <  1:3927142/100‑5 (MQ=255)
                              ttCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGatat                                                    <  1:1045852/100‑1 (MQ=255)
                                  ggAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGgtg                                                <  1:3693936/100‑1 (MQ=255)
                                   gAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATggagc                                               <  1:4374569/100‑4 (MQ=255)
                                    tgCATCGGTTAAGTCTATGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGtt                                              <  1:3103653/99‑1 (MQ=255)
                                    aGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGtt                                              <  1:1786091/100‑1 (MQ=255)
                                    aGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGtt                                              <  1:963343/100‑1 (MQ=255)
                                    aGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGtt                                              <  1:5515864/100‑1 (MQ=255)
                                    aGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGtt                                              <  1:1115301/100‑1 (MQ=255)
                                    aGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGtt                                              <  1:3710942/100‑1 (MQ=255)
                                    aGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGtt                                              <  1:4365484/100‑1 (MQ=255)
                                    aGCATCTTTTAAATCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGtt                                              <  1:4829696/100‑1 (MQ=255)
                                     gCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGgtgtcg                                             <  1:3159638/100‑3 (MQ=255)
                                     gCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGttt                                             <  1:5415050/100‑2 (MQ=255)
                                     gCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTa                                             <  1:4445499/100‑1 (MQ=255)
                                         cTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTcc                                         <  1:2867367/100‑1 (MQ=255)
                                            ttAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCAtt                                      <  1:2652601/100‑1 (MQ=255)
                                            ttAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCAtg                                      <  1:3930614/100‑2 (MQ=255)
                                             tAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCGCGATATGGTGTTAGTCCATTg                                     <  1:5768480/100‑1 (MQ=255)
                                             tAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTg                                     <  1:5200882/100‑1 (MQ=255)
                                             tAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTg                                     <  1:3745964/100‑1 (MQ=255)
                                             ccgatCTACGTCTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGCTAGTCCATTg                                     <  1:4753749/96‑1 (MQ=255)
                                                 tCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTGatag                                 <  1:762585/100‑2 (MQ=255)
                                                       tttATATACATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTGATATCTTTgg                           <  1:5775090/100‑2 (MQ=255)
                                                             taCATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTGATATCTTTGAAGTGAt                     <  1:554632/100‑1 (MQ=255)
                                                              aCATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTGATATCTTTGAAGTGAtg                    <  1:1234926/100‑2 (MQ=255)
                                                              aCATATCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTGATATCTTTGAAGTGATa                    <  1:3074123/100‑1 (MQ=255)
                                                                atatCAGGATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTGATATCTTTGAAGTGATGAg                  <  1:4614292/100‑1 (MQ=255)
                                                                      ggATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTGATATCTTTGAAGTGATAAGCTGCgg            <  1:5137560/100‑2 (MQ=255)
                                                                      ggATATTGCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTGATATCTTTGAAGTGATAAGCTGCGt            <  1:4715222/100‑1 (MQ=255)
                                                                         tctggCTTCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCGTTGATATCTTTGAAGTGATAAGCTGCGTCTa         <  1:4473868/96‑1 (MQ=255)
                                                                                tCAGCAACTCATCGGAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTGATATCTTTGAAGTGATAAGCTGCGTCTAATTTCtt  <  1:4314565/100‑1 (MQ=255)
                                                                                              |                                                                                     
AATTTGTGGTGCGATTGAGCCTTCATAGGATTCAGGAGCATCTTTTAAGTCTACGTTTATATACATATCAGGATATTGCTTCAGCAACTCATCGAAGGTTAGTATAGCTGTGTGTGCATGACCACGATATGGTGTTAGTCCATTGATATCTTTGAAGTGATAAGCTGCGTCTAATTTCTT  >  NC_007793/1117554‑1117733

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑