Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,828,990 G→A 14.3% H328H (CAC→CAT astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,828,9900GA14.3% 167.0 / 17.1 63H328H (CAC→CATastBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base G (30/24);  new base A (6/3);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

AAAGACTTTTAGTTCGCTAATCGGGTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAA  >  W3110S.gb/1828921‑1829060
                                                                     |                                                                      
aaaGACTTTTAGTTCGCTAATCGGGTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCg                                                                        <  1:1234063/70‑1 (MQ=255)
        ttAGTTCGCTAATCGGGTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTccc                                                                 <  1:840653/69‑1 (MQ=255)
        ttAGTTCGCTAATCGGGTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTccc                                                                 <  1:4433694/69‑1 (MQ=255)
         tAGTTCGCTAATCGGGTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGa                                                               <  1:2094366/70‑1 (MQ=255)
          aGTTCGCTAATCGGGTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGaca                                                             >  1:969239/1‑71 (MQ=255)
              cGCTAATCGGGTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCATGTTCCCGACACTcc                                                         >  1:469637/1‑71 (MQ=255)
                cTAATCGGGTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTg                                                       <  1:397551/71‑1 (MQ=255)
                   aTCGGGTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAgg                                                    >  1:3497125/1‑71 (MQ=255)
                   aTCGGGTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAg                                                     <  1:1727469/70‑1 (MQ=255)
                     cGGGTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAgg                                                    >  1:1851271/1‑69 (MQ=255)
                      gggTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCATGTTCCCGACACTCCTGAGgcag                                                 >  1:3168421/1‑71 (MQ=255)
                        gTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCATGTTCCCGACACTCCTGAGgcagca                                               <  1:3524345/71‑1 (MQ=255)
                           gTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGgcagcagc                                             <  1:3720481/70‑3 (MQ=255)
                            tCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCt                                                        <  1:1537710/58‑1 (MQ=255)
                             cAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCaa                                           >  1:3082862/1‑70 (MQ=255)
                              aGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAAca                                         >  1:2700351/1‑71 (MQ=255)
                                cGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAAcat                                        >  1:3684028/1‑70 (MQ=255)
                                 ggCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGgcagc                                                <  1:1403466/61‑1 (MQ=255)
                                  gCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACt                                                           <  1:2920335/49‑1 (MQ=255)
                                          ttCATTGAGATAACCCCATACTCCGGCGTGTTccc                                                                 <  1:295805/35‑1 (MQ=255)
                                          ttCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGgcagca                                               >  1:4631719/1‑53 (MQ=255)
                                          ttCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCAc                                              >  1:1969704/1‑54 (MQ=255)
                                          ttCATTGAGATAACCCCATACTCCAGCGTGTTCCCGACACTCCTGAGGCAGCAc                                              >  1:4494207/1‑54 (MQ=255)
                                            cATTGAGATAACCCCATACTCCGGCGTGTTCCCGacac                                                            >  1:2952133/1‑38 (MQ=255)
                                            cATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACcatca                           >  1:2820/1‑71 (MQ=255)
                                            cATTGAGATAACCCCATACTCCGGCATGTTCCCGACACTCCTGAGGCAGCACCAAc                                          >  1:1892388/1‑56 (MQ=255)
                                            cATTGAGATAACCCCATACTCCGGCATGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACcatc                            >  1:1188505/1‑70 (MQ=255)
                                             aTTGAGATAACCCCATACTCCGGCATTTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACcatca                           <  1:3871718/70‑1 (MQ=255)
                                              ttGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGa                                                      <  1:3202744/42‑1 (MQ=255)
                                               tcagaTAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATcg                        >  1:3643298/3‑71 (MQ=255)
                                               tGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAgg                                                    >  1:2855720/1‑43 (MQ=255)
                                               tGAGATAACCCCATACTCCGGCATGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATcg                        >  1:4606293/1‑71 (MQ=255)
                                                gagaTAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGgcagca                                               >  1:28985/1‑47 (MQ=255)
                                                gagaTAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGgcagca                                               >  1:1797596/1‑47 (MQ=255)
                                                gagaTAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGgc                                                   >  1:250763/1‑43 (MQ=255)
                                                gagaTAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGa                                  >  1:3581321/1‑60 (MQ=255)
                                                gagaTAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGa                                  >  1:3284876/1‑60 (MQ=255)
                                                     aaCACCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACcatca                           >  1:4258937/1‑62 (MQ=255)
                                                      aCCCCATACTCCGGGGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCt                   >  1:2065632/1‑69 (MQ=255)
                                                      aCCCCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCt                   >  1:786832/1‑69 (MQ=255)
                                                      aCCCCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCt                   >  1:1301830/1‑69 (MQ=255)
                                                        cccATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTcagc               <  1:1218664/71‑1 (MQ=255)
                                                         ccATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTcagca              <  1:2095633/71‑1 (MQ=255)
                                                          cATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTcagcag             >  1:3209823/1‑71 (MQ=255)
                                                          cATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTcagcag             >  1:2327301/1‑71 (MQ=255)
                                                          cATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTcagcag             >  1:2569908/1‑71 (MQ=255)
                                                            tACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCAcc                                             <  1:327167/37‑1 (MQ=255)
                                                            tACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCAcc                                             <  1:3126728/37‑1 (MQ=255)
                                                            tACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAAcatcat                                     <  1:2481033/45‑1 (MQ=255)
                                                            tACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACcatc                            <  1:3340025/54‑1 (MQ=255)
                                                            tACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACcat                             <  1:2818843/53‑1 (MQ=255)
                                                            tACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACca                              <  1:1527076/52‑1 (MQ=255)
                                                              cTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCaa                                           >  1:3298354/1‑37 (MQ=255)
                                                              cTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTgg         >  1:1534355/1‑71 (MQ=255)
                                                              cTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTgg         >  1:4099445/1‑71 (MQ=255)
                                                              cTCCGGCATGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTg          <  1:3183156/70‑1 (MQ=255)
                                                               tCCGTCATGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGc        >  1:730389/1‑71 (MQ=255)
                                                                 cGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGc        <  1:2672445/69‑1 (MQ=255)
                                                                 cGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGc        <  1:3980129/69‑1 (MQ=255)
                                                                 cGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCt       >  1:556252/1‑70 (MQ=255)
                                                                  ggCGTGTTCCCGACACTCCTGAGGCAGCACCATCATCATGGAACCATCATCGCGGCTCAGCAGTTGGc        <  1:2478899/68‑1 (MQ=255)
                                                                  ggCGGGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGc        <  1:875396/68‑1 (MQ=255)
                                                                    cGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTg          >  1:1154728/1‑64 (MQ=255)
                                                                    cGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTa   <  1:4077354/71‑1 (MQ=255)
                                                                     gtgtTCCCGACACTCCTGAGGCAGCACCAACATCATgg                                   >  1:1920728/1‑38 (MQ=255)
                                                                     gtgtTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTaa  <  1:2058674/71‑1 (MQ=255)
                                                                     |                                                                      
AAAGACTTTTAGTTCGCTAATCGGGTTGTCAGCGGCAAGGAGTTCATTGAGATAACCCCATACTCCGGCGTGTTCCCGACACTCCTGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAA  >  W3110S.gb/1828921‑1829060

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