Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,326,354 G→A 6.9% P181P (CCG→CCA yciO → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,326,3540GA 6.9% 175.2 / 2.4 58P181P (CCG→CCAyciOconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/24);  new base A (2/2);  total (32/26)
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.90e-01

TAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGC  >  W3110S.gb/1326291‑1326421
                                                               |                                                                   
tAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACgg                                                              >  1:1713332/1‑71 (MQ=255)
 aaaGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACgg                                                              <  1:3244087/70‑1 (MQ=255)
  aaGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGa                                                                    >  1:1551500/1‑63 (MQ=255)
     aTCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTAtt                                                         <  1:191377/71‑1 (MQ=255)
     aTCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTAtt                                                         <  1:3333981/71‑1 (MQ=255)
      tCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTAtt                                                         <  1:1804778/70‑1 (MQ=255)
       cGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGa                                                       >  1:1426579/1‑71 (MQ=255)
              aaaaaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACaa                                                                 >  1:3527438/1‑54 (MQ=255)
              aaaaaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTAcc                                                <  1:2018746/71‑1 (MQ=255)
              aaaaaCAGGTAGATTTGATTATTCACGGCGGATATCTCGGGCAGAAACCGACAACGGTTATTGATCTTAcc                                                <  1:779749/71‑1 (MQ=255)
               aaaaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTAcc                                                >  1:2128220/1‑70 (MQ=255)
               aaaaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCg                                               <  1:3374297/71‑1 (MQ=255)
               aaaaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCg                                               <  1:3720123/71‑1 (MQ=255)
               aaaaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCg                                               <  1:3065664/71‑1 (MQ=255)
                 aaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCgat                                             <  1:3694715/71‑1 (MQ=255)
                 aaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCgat                                             <  1:2020029/71‑1 (MQ=255)
                  aCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCgat                                             >  1:3393133/1‑70 (MQ=255)
                   cAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCgatga                                           <  1:3541087/71‑1 (MQ=255)
                     ggTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGa                                                                    >  1:2653564/1‑44 (MQ=255)
                         gATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGcc                                     >  1:3205344/1‑71 (MQ=255)
                          aTTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTg                                                        <  1:995121/51‑1 (MQ=255)
                          aTTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGcc                                     <  1:605027/70‑1 (MQ=255)
                          aTTTGATTATTCACGGCGGATATCTCGGGCAGAAACCGACAACGGTTATTg                                                        <  1:3225601/51‑1 (MQ=255)
                           tttGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTa                                                  <  1:4114224/56‑1 (MQ=255)
                           tttGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCgg                                   >  1:1406834/1‑71 (MQ=255)
                           tttGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCgg                                   >  1:695915/1‑71 (MQ=255)
                           tttGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCg                                    >  1:3710994/1‑70 (MQ=255)
                               attattCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCgt                               <  1:3189288/71‑1 (MQ=255)
                                ttattCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCgatga                                           <  1:2872410/58‑1 (MQ=255)
                                ttattCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCgtg                              >  1:2230375/1‑71 (MQ=255)
                                  attCACGGCGGCTATCTCGGGCAGAAACCGACAACg                                                               >  1:4356484/1‑36 (MQ=255)
                                  attCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTAtt                                                         >  1:948028/1‑42 (MQ=255)
                                      aCGGCGGCTATCTCGGGCAGAAACCAACAACGGTTATTGATCTCACCga                                              >  1:1057407/1‑49 (MQ=38)
                                      aCGGCGGCTATCTCGGGCAGAAACCAACAACGGTTATTGATCTCACCga                                              >  1:1059750/1‑49 (MQ=38)
                                        ggcggcTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTg                       >  1:1172427/1‑70 (MQ=255)
                                        ggcggcTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGa                      <  1:3453700/71‑1 (MQ=255)
                                         gcggcTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGaa                     <  1:1032621/71‑1 (MQ=255)
                                               aTCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCg                 >  1:3513695/1‑69 (MQ=255)
                                                tctcGGGCAGAAACCGACAACGGTTATTGATCTTACCga                                              >  1:1710476/1‑39 (MQ=255)
                                                tctcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAg              >  1:2621431/1‑71 (MQ=255)
                                                tctcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAg              >  1:1461089/1‑71 (MQ=255)
                                                 ctcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAg              >  1:462677/1‑70 (MQ=255)
                                                  tcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTa               >  1:723199/1‑68 (MQ=255)
                                                  tcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTa               >  1:718100/1‑68 (MQ=255)
                                                  tcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTa               >  1:3040184/1‑68 (MQ=255)
                                                  tcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGt            >  1:3474913/1‑71 (MQ=255)
                                                  tcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGt            >  1:2814724/1‑71 (MQ=255)
                                                  tcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGt            >  1:613959/1‑71 (MQ=255)
                                                  tcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGt            >  1:581173/1‑71 (MQ=255)
                                                    gggCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGt            <  1:101887/69‑1 (MQ=255)
                                                       cAGAAACCGACAACGGTTATTGATCTTAGCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGAtgt       >  1:2885557/1‑71 (MQ=255)
                                                        aGAAACCGACAACGGTTATTGATCTTACCGATTATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGAtgtg      <  1:528351/71‑1 (MQ=255)
                                                        aGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGAtgtg      <  1:1429119/71‑1 (MQ=255)
                                                        aGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGAtgt       >  1:175690/1‑70 (MQ=255)
                                                        aGAAACCAACAACGGTTATTGATCTCACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGAtgtg      <  1:2584153/71‑1 (MQ=255)
                                                        aGAAACCAACAACGGTTATTGATCTCACCGATGATACGCCGGTCGTGGTGCGTCAAGGCGTAGGTGAtgtg      <  1:1654380/71‑1 (MQ=39)
                                                         gAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGa                      <  1:2102813/54‑1 (MQ=255)
                                                            aCCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGc  >  1:2712368/1‑71 (MQ=255)
                                                               |                                                                   
TAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGC  >  W3110S.gb/1326291‑1326421

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: