Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,915,033 T→C 31.2% intergenic (+29/‑66) ysaB → / → glyQ hypothetical protein/glycine tRNA synthetase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,915,0330TC31.2% 82.7 / 55.3 64intergenic (+29/‑66)ysaB/glyQhypothetical protein/glycine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (20/24);  new base C (8/12);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAAATAATATGCA  >  W3110S.gb/3914967‑3915103
                                                                  |                                                                      
aCGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATc                                                                     <  1:3023552/70‑1 (MQ=255)
aCGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCt                                                                    <  1:2209655/71‑1 (MQ=255)
  gcgGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtc                                                                  <  1:4036000/71‑1 (MQ=255)
            aTTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCgg                                                        >  1:3314985/1‑71 (MQ=255)
            aTTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCgg                                                        >  1:4096296/1‑71 (MQ=255)
            aTTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCgg                                                        >  1:2642391/1‑71 (MQ=255)
             ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtccat                                                               <  1:4630467/63‑1 (MQ=255)
             ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtccat                                                               <  1:1754150/63‑1 (MQ=255)
             ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtccat                                                               <  1:2372657/63‑1 (MQ=255)
              tttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtccatcca                                                            >  1:3578575/1‑65 (MQ=255)
               ttttGCATGTTTCCATGCGTTAAGCTCTGCTTTTTCCCGTTTCGTACTGTacat                                                                      >  1:1431297/1‑54 (MQ=255)
               ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTacat                                                                      >  1:2875569/1‑54 (MQ=255)
               ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                     >  1:1787487/1‑71 (MQ=255)
               ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATc                                                                     >  1:2667539/1‑55 (MQ=255)
               ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGta                                                     >  1:1864590/1‑71 (MQ=255)
               ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGta                                                     >  1:2996543/1‑71 (MQ=255)
                tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                     >  1:3941380/1‑70 (MQ=255)
                tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                     <  1:3438760/70‑1 (MQ=255)
                tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                     >  1:5038911/1‑70 (MQ=255)
                tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                     <  1:1521399/70‑1 (MQ=255)
                tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                     >  1:2311354/1‑70 (MQ=255)
                tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGta                                                     <  1:2716401/70‑1 (MQ=255)
                tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGta                                                     <  1:1794922/70‑1 (MQ=255)
                   gCATCTTTCCATGCGTTAAGCCTTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATAc                                                  >  1:718155/1‑70 (MQ=255)
                     aTCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGtata                                                   <  1:4919944/67‑1 (MQ=255)
                     aTCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGtata                                                   <  1:1029262/67‑1 (MQ=255)
                       cTTTCCATGCGTTAAGCCTTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCggg                                                       <  1:653318/61‑1 (MQ=255)
                       cTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAg                                                           <  1:1378241/57‑1 (MQ=255)
                       cTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAg                                                           <  1:2695137/57‑1 (MQ=255)
                       cTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAg                                                           <  1:2289176/57‑1 (MQ=255)
                        tttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACt                                                 >  1:2865101/1‑66 (MQ=255)
                         ttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCtt                                                                   >  1:1395177/1‑47 (MQ=255)
                         ttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCg                                                         >  1:2240750/1‑57 (MQ=255)
                          tCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTAGTGTATATCTTCCATCCAGCGGGTATACTGATccc                                           >  1:2125753/1‑70 (MQ=255)
                                 gTTAAGCCTTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTa                                   <  1:1947958/71‑1 (MQ=255)
                                    aaGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGtata                                                   >  1:3182458/1‑52 (MQ=255)
                                    aaGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAAt                                >  1:1718781/1‑71 (MQ=255)
                                         cTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                     <  1:3475198/45‑1 (MQ=255)
                                          tGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTa                                   <  1:3824334/62‑1 (MQ=255)
                                           gCTTTTTCCCGTTTCGTATTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCca                             >  1:3520034/1‑67 (MQ=255)
                                           gCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                          <  1:3583325/38‑1 (MQ=255)
                                           gCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                          <  1:391563/38‑1 (MQ=255)
                                           gCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                          <  1:639576/38‑1 (MQ=255)
                                           gCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                          <  1:701590/38‑1 (MQ=255)
                                           gCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                          <  1:1897285/38‑1 (MQ=255)
                                            cTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTaa                                  >  1:3737020/1‑61 (MQ=255)
                                            cTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACg                         >  1:3660285/1‑70 (MQ=255)
                                            cTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGtata                                                   >  1:1259356/1‑44 (MQ=255)
                                             tttttCCCGTTTCGTACTGTATATCTTCCATCCAGCgg                                                        >  1:794174/1‑38 (MQ=255)
                                             tttttCCCGTTTCGTACTGTATATCTTCCATCCAGCgg                                                        >  1:2491649/1‑38 (MQ=255)
                                               tttCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATc                     >  1:2880616/1‑71 (MQ=255)
                                                ttCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCttt                                    <  1:2841194/55‑1 (MQ=255)
                                                ttCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCttt                                    <  1:4771460/55‑1 (MQ=255)
                                                ttCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCcaca                           >  1:4013416/1‑64 (MQ=255)
                                                 tCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATc                     >  1:4357238/1‑69 (MQ=255)
                                                   ccGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAg                  <  1:195782/70‑1 (MQ=255)
                                                   ccGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAg                  <  1:4744929/70‑1 (MQ=255)
                                                   ccGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAg                  <  1:2925350/70‑1 (MQ=255)
                                                   ccGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAg                  <  1:763311/70‑1 (MQ=255)
                                                   ccGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAg                  <  1:776746/70‑1 (MQ=255)
                                                    cGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAg                  <  1:2756310/69‑1 (MQ=255)
                                                     gTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACTCGTATCCAg                  <  1:3581452/68‑1 (MQ=255)
                                                      tttCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATc                     <  1:4814218/64‑1 (MQ=255)
                                                      tttCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATc                     <  1:4541050/64‑1 (MQ=255)
                                                      tttCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCAc               >  1:3622412/1‑70 (MQ=255)
                                                              tGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCttt                                    <  1:3225827/41‑1 (MQ=255)
                                                                tataTCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTaa                                  >  1:2548691/1‑41 (MQ=255)
                                                                 atatCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAAATAATATg    >  1:4939881/1‑70 (MQ=255)
                                                                  tatCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAAATAATATGCa  <  1:32681/71‑1 (MQ=255)
                                                                  tatCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAAATAATATGCa  <  1:676293/71‑1 (MQ=255)
                                                                  |                                                                      
ACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAAATAATATGCA  >  W3110S.gb/3914967‑3915103

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: