Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,326,372 T→C 8.5% L187L (CTT→CTC yciO → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,326,3720TC 8.5% 152.6 / 5.7 59L187L (CTT→CTCyciOconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/26);  new base C (2/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

AAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGC  >  W3110S.gb/1326305‑1326440
                                                                   |                                                                    
aaaaaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTAcc                                                                   <  1:2018746/71‑1 (MQ=255)
aaaaaCAGGTAGATTTGATTATTCACGGCGGATATCTCGGGCAGAAACCGACAACGGTTATTGATCTTAcc                                                                   <  1:779749/71‑1 (MQ=255)
 aaaaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTAcc                                                                   >  1:2128220/1‑70 (MQ=255)
 aaaaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCg                                                                  <  1:3065664/71‑1 (MQ=255)
 aaaaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCg                                                                  <  1:3374297/71‑1 (MQ=255)
 aaaaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCg                                                                  <  1:3720123/71‑1 (MQ=255)
   aaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCgat                                                                <  1:2020029/71‑1 (MQ=255)
   aaCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCgat                                                                <  1:3694715/71‑1 (MQ=255)
    aCAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCgat                                                                >  1:3393133/1‑70 (MQ=255)
     cAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCgatga                                                              <  1:3541087/71‑1 (MQ=255)
           gATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGcc                                                        >  1:3205344/1‑71 (MQ=255)
            aTTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGcc                                                        <  1:605027/70‑1 (MQ=255)
             tttGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTa                                                                     <  1:4114224/56‑1 (MQ=255)
             tttGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCgg                                                      >  1:695915/1‑71 (MQ=255)
             tttGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCgg                                                      >  1:1406834/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)
                        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:3040184/1‑68 (MQ=255)
                                    tcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTa                                  >  1:723199/1‑68 (MQ=255)
                                    tcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTa                                  >  1:718100/1‑68 (MQ=255)
                                    tcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGt                               >  1:2814724/1‑71 (MQ=255)
                                    tcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGt                               >  1:613959/1‑71 (MQ=255)
                                    tcGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGt                               >  1:3474913/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)
                                                  aCAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGaa                                        >  1:58613/1‑48 (MQ=255)
                                                  aCAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCttt                 >  1:3912692/1‑71 (MQ=255)
                                                  aCAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCttt                 >  1:4189108/1‑71 (MQ=255)
                                                   cAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTc                <  1:2031590/71‑1 (MQ=255)
                                                   cAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTc                <  1:2150797/71‑1 (MQ=255)
                                                   cAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTc                <  1:1476534/71‑1 (MQ=255)
                                                    aaCGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGAtgt                          <  1:646555/60‑1 (MQ=255)
                                                        gTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGa                        <  1:956686/58‑1 (MQ=255)
                                                           aTTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAAc         >  1:2716603/1‑70 (MQ=255)
                                                            ttGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCt               <  1:3219306/63‑1 (MQ=255)
                                                            ttGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACaa       <  1:1073400/71‑1 (MQ=255)
                                                            ttGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACa        <  1:971093/70‑1 (MQ=255)
                                                            ttGATCTCACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGAtgtg                         <  1:156497/53‑1 (MQ=255)
                                                                 cTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGc  >  1:2517068/1‑71 (MQ=255)
                                                                  ttACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAgg                                <  1:3048869/40‑1 (MQ=255)
                                                                  ttACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACaaa      >  1:3187120/1‑66 (MQ=255)
                                                                  ttACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACaa       <  1:2102557/65‑1 (MQ=255)
                                                                  ttACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAAt     <  1:1840511/67‑1 (MQ=255)
                                                                   |                                                                    
AAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGC  >  W3110S.gb/1326305‑1326440

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: