Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,341,910 A→G 8.6% G727G (GGT→GGC gyrA ← DNA gyrase (type II topoisomerase), subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,341,9100AG 8.6% 145.1 / 5.8 58G727G (GGT→GGCgyrADNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (24/29);  new base G (2/3);  total (26/32)
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.22e-01

CGTTTTGCGTTGCGGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTAA  >  W3110S.gb/2341840‑2341974
                                                                      |                                                                
cGTTTTGCGTTGCGGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCa                                                                  >  1:2731638/1‑71 (MQ=255)
   tttGCGTTGCGGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCAcc                                                                >  1:4054356/1‑70 (MQ=255)
    ttgcgttgcgGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCt                                                               <  1:1233964/70‑1 (MQ=255)
    ttgcgttgcgGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCt                                                               >  1:1933867/1‑70 (MQ=255)
        gttgcgGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCg                                                          <  1:754245/71‑1 (MQ=255)
         ttgcgGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCg                                                          <  1:1130317/70‑1 (MQ=255)
         ttgcgGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCg                                                          <  1:4036800/70‑1 (MQ=255)
         ttgcgGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCg                                                          <  1:1850881/70‑1 (MQ=255)
           gcgGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTa                                                              >  1:4199825/1‑64 (MQ=255)
           gcgGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGaa                                                        <  1:883350/70‑1 (MQ=255)
           gcgGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGaa                                                        <  1:878318/70‑1 (MQ=255)
                 aGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACcgcg                                                 <  1:21000/71‑1 (MQ=255)
                  ggATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACcgcg                                                 <  1:1308687/70‑1 (MQ=255)
                     ttGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCGCCTAAACGAATACCACGAaca                                             >  1:4277951/1‑71 (MQ=255)
                     ttGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCGCCTAAACGAATACCACGAaca                                             >  1:4283143/1‑71 (MQ=255)
                      tGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCa                                                                  <  1:2075543/49‑1 (MQ=255)
                        cgcCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACAcc                                           <  1:1543110/70‑1 (MQ=255)
                          ccATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCggt                                        <  1:3118046/71‑1 (MQ=255)
                           cATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCggtg                                       >  1:2531845/1‑71 (MQ=255)
                             tCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCggtggt                                     >  1:979756/1‑71 (MQ=255)
                               gCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCAcc                                                                >  1:3579130/1‑42 (MQ=255)
                               gCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGgtgt                                   <  1:2236340/71‑1 (MQ=255)
                                  aCGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGaa                                                        >  1:1124007/1‑47 (MQ=255)
                                  aCGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGaa                                                        >  1:1745606/1‑47 (MQ=255)
                                  aCGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTg                                 >  1:3211023/1‑70 (MQ=255)
                                   cGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAAt                                                       <  1:1202307/47‑1 (MQ=255)
                                    gAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGaa                                                        >  1:2145536/1‑45 (MQ=255)
                                     aGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATAc                                                     >  1:1084278/1‑47 (MQ=255)
                                     aGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATAc                                                     >  1:704355/1‑47 (MQ=255)
                                     aGGCACGATCAGAGAGACGACTTTATCGCCTTAACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGc                             >  1:1941865/1‑71 (MQ=255)
                                      ggCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCggtggt                                     >  1:3037985/1‑62 (MQ=255)
                                       tcACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCa                               <  1:3822439/66‑1 (MQ=255)
                                       gCACGATCAGAGAGACGACTTTATCGCCTTCGCCTAAACGAATACCACGAACACCGGTGGTGTTGCAGccc                           <  1:3668633/71‑1 (MQ=255)
                                       gCACGATCAGAGAGACGACTTTATCGCCTTCACCTa                                                              <  1:2428514/36‑1 (MQ=255)
                                       gCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCggt                                        <  1:3493553/58‑1 (MQ=255)
                                       gCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCggt                                        <  1:3498472/58‑1 (MQ=255)
                                         aCGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCg                                                          <  1:1208932/38‑1 (MQ=255)
                                         aCGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGaa                                                        <  1:1154533/40‑1 (MQ=255)
                                         aCGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGaa                                                        >  1:1604364/1‑40 (MQ=255)
                                         aCGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCggtg                                       >  1:4045489/1‑57 (MQ=255)
                                           gATCAGAGAGACGACTTTATCGCCTTCACCTAAGCg                                                          <  1:2408132/36‑1 (MQ=255)
                                           gATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCg                       <  1:2969344/71‑1 (MQ=255)
                                            aTCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGc                      <  1:4418235/71‑1 (MQ=255)
                                              cAGAGAGACGACTTTATCGCCTTCACCTAAGCGaa                                                        >  1:447526/1‑35 (MQ=255)
                                              cAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCAc                    >  1:248444/1‑71 (MQ=255)
                                                  gagaCGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGa                 >  1:1239629/1‑70 (MQ=255)
                                                  gagaCGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGAc                <  1:2127340/71‑1 (MQ=255)
                                                   agaCGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGgtgt                                   <  1:2973484/51‑1 (MQ=255)
                                                   agaCGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCg                       <  1:3084757/63‑1 (MQ=255)
                                                      cgacTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGc                      <  1:733138/61‑1 (MQ=255)
                                                      cgacTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACagaa            <  1:2234132/71‑1 (MQ=255)
                                                            tATCGCCTTCGCCTAAACGAATACCACGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGAct        <  1:3721397/69‑1 (MQ=255)
                                                            tATCGCCTTCGCCTAAACGAATACCACGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGAct        <  1:3092603/69‑1 (MQ=255)
                                                            tATCGCCTTCACCTAAGCGAATACCGCGAACACCggt                                        >  1:401298/1‑37 (MQ=255)
                                                            tATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCa                               >  1:2677080/1‑46 (MQ=255)
                                                            tATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCa                               >  1:616134/1‑46 (MQ=255)
                                                            tATCGCCTTCACCTAAGCGAATACAGCGAACACCGGTGGTGTTGc                                >  1:1933178/1‑45 (MQ=255)
                                                              tCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCtt     >  1:3396422/1‑70 (MQ=255)
                                                                gCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTaa  <  1:1289572/71‑1 (MQ=255)
                                                                gCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTa   <  1:4246012/70‑1 (MQ=255)
                                                                      |                                                                
CGTTTTGCGTTGCGGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTAA  >  W3110S.gb/2341840‑2341974

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: