Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,341,916 G→A 8.0% R725R (CGC→CGT 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,9160GA 8.0% 184.3 / 5.1 62R725R (CGC→CGTgyrADNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base A (2/3);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GTTGCGGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTAAAGCGCACCACTT  >  W3110S.gb/2341848‑2341986
                                                                    |                                                                      
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)
   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:4283143/1‑71 (MQ=255)
             ttGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCGCCTAAACGAATACCACGAaca                                                         >  1:4277951/1‑71 (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)
                       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)
                               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:3092603/69‑1 (MQ=255)
                                                    tATCGCCTTCGCCTAAACGAATACCACGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGAct                    <  1:3721397/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)
                                                               ccTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTAAAGCGcac       >  1:3818918/1‑71 (MQ=255)
                                                                 tAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCa                                 >  1:639064/1‑43 (MQ=255)
                                                                 tAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCAc                                >  1:3483744/1‑44 (MQ=255)
                                                                 tAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACagaaga                      >  1:2816574/1‑54 (MQ=255)
                                                                 tAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTAAAg            >  1:558506/1‑64 (MQ=255)
                                                                 tAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTAAAGCGcacc      >  1:17233/1‑70 (MQ=255)
                                                                  aaGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTAAAGCGcaccac    >  1:1262262/1‑71 (MQ=255)
                                                                  aaGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTAAAGCGcaccac    >  1:2414796/1‑71 (MQ=255)
                                                                  aaGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTAAAGCGcaccac    >  1:2462954/1‑71 (MQ=255)
                                                                  aaGCGAATACCGAGAACACCGGTGGTGTTGCAGCCCATCGCACGGACagaaga                      >  1:814238/1‑53 (MQ=255)
                                                                    gCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTAAAGCGCACCACtt  <  1:3078075/71‑1 (MQ=255)
                                                                    |                                                                      
GTTGCGGTGAGGATTGCGCCATCGCCACGAGGCACGATCAGAGAGACGACTTTATCGCCTTCACCTAAGCGAATACCGCGAACACCGGTGGTGTTGCAGCCCATCGCACGGACAGAAGACTCTTTAAAGCGCACCACTT  >  W3110S.gb/2341848‑2341986

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: