Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,340,293 C→G 6.9% L194L (CTG→CTC) ‡ ribA ← GTP cyclohydrolase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,340,2930CG 6.9% 260.8 / 6.7 87L194L (CTG→CTC) ‡ribAGTP cyclohydrolase II
Reads supporting (aligned to +/- strand):  ref base C (48/32);  new base G (3/3);  total (52/35)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

TATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTACAA  >  W3110S.gb/1340225‑1340360
                                                                    |                                                                   
tATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAATTGAATTACACAATGCAAGAGGGTTATTTGTTCa                                                                    >  1:3292158/1‑70 (MQ=255)
   ttgcttgcCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACCCAATGCAAGAGGGTTATTTGTTCAGCa                                                                 >  1:3167364/1‑70 (MQ=255)
     gcttgcCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAAt                                                              >  1:3210961/1‑71 (MQ=255)
      cttgcCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCa                                                                 <  1:7219618/67‑1 (MQ=255)
      cttgcCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATg                                                             <  1:6767150/71‑1 (MQ=255)
        tgcCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGc                                                           >  1:6071155/1‑71 (MQ=255)
         gcCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATg                                                             >  1:6429116/1‑68 (MQ=255)
          ccGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCa                                                                 <  1:795270/63‑1 (MQ=255)
          ccGGTCATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGccc                                                         <  1:4050754/71‑1 (MQ=255)
              ttATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCaaa                                                               >  1:863737/1‑61 (MQ=255)
              ttATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCaaa                                                               >  1:5327763/1‑61 (MQ=255)
              ttATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCAttt                                                     >  1:5527855/1‑71 (MQ=255)
               tATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCAttt                                                     >  1:4312907/1‑70 (MQ=255)
               tATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCAttt                                                     >  1:5609298/1‑70 (MQ=255)
                 ttttGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGc                                                           >  1:4536052/1‑62 (MQ=255)
                 ttttGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTGATTTGTTCAGCAAATGGCCCATTTtc                                                   <  1:187716/70‑1 (MQ=255)
                  tttGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTtc                                                   <  1:31619/69‑1 (MQ=255)
                    tGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCa                                                                    >  1:7404150/1‑50 (MQ=255)
                      cTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGc                                              <  1:6028568/70‑1 (MQ=255)
                       ttCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCa                                                                 <  1:7442397/50‑1 (MQ=255)
                         ccGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCtt                                            <  1:5031449/69‑1 (MQ=255)
                         ccGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCtt                                            <  1:68148/69‑1 (MQ=255)
                          cGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTgg                                         >  1:5322413/1‑71 (MQ=255)
                           ggCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAAt                                                              >  1:5816580/1‑49 (MQ=255)
                            gcaagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTa                                       >  1:5138689/1‑71 (MQ=255)
                            gcaagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTa                                       >  1:4222810/1‑71 (MQ=255)
                             caagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCgg                                               >  1:3297265/1‑62 (MQ=255)
                             caagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCgg                                               >  1:1208512/1‑62 (MQ=255)
                             caagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTAt                                      >  1:1889331/1‑71 (MQ=255)
                              aagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATc                                     >  1:6346793/1‑71 (MQ=255)
                               agcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTtc                                                   <  1:7173522/56‑1 (MQ=255)
                               agcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATc                                     >  1:1166856/1‑70 (MQ=255)
                                 caaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCt                                             <  1:1412664/60‑1 (MQ=255)
                                  aaaTGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTa                                       <  1:2029671/65‑1 (MQ=255)
                                  aaaTGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTAta                                     >  1:6725543/1‑66 (MQ=255)
                                  aaaTGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCgaga                                 <  1:6332372/71‑1 (MQ=255)
                                    aTGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGAt                                >  1:5917228/1‑70 (MQ=255)
                                    aTGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGAt                                >  1:603649/1‑70 (MQ=255)
                                     tGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGAt                                <  1:2659261/69‑1 (MQ=255)
                                      gAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCAt                                                       <  1:4263042/45‑1 (MQ=255)
                                        aTTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:4278258/1‑71 (MQ=255)
                                         ttACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGtt                          >  1:4354366/1‑71 (MQ=255)
                                         ttACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           <  1:4396341/70‑1 (MQ=255)
                                              caATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTa                     <  1:3663820/71‑1 (MQ=255)
                                              caATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTa                     >  1:6171478/1‑71 (MQ=255)
                                                aTGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAtt                   <  1:6781157/71‑1 (MQ=255)
                                                aTGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAtt                   >  1:7755092/1‑71 (MQ=255)
                                                aTGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAtt                   <  1:6124488/71‑1 (MQ=255)
                                                 tGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTg                  >  1:7213777/1‑71 (MQ=255)
                                                 tGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTg                  >  1:6683032/1‑71 (MQ=255)
                                                 tGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTg                  >  1:3972281/1‑71 (MQ=255)
                                                 tGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTg                  >  1:2976223/1‑71 (MQ=255)
                                                   cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAAt                             >  1:7473043/1‑58 (MQ=255)
                                                   cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAAt                             >  1:303449/1‑58 (MQ=255)
                                                   cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAAt                             >  1:3096978/1‑58 (MQ=255)
                                                   cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAAt                             >  1:4861256/1‑58 (MQ=255)
                                                   cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:2889967/1‑60 (MQ=255)
                                                   cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:7122809/1‑60 (MQ=255)
                                                   cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:2487549/1‑60 (MQ=255)
                                                   cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:2088155/1‑60 (MQ=255)
                                                   cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:6350047/1‑60 (MQ=255)
                                                     agagGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCtt                                            <  1:5126215/41‑1 (MQ=255)
                                                      gagGGTTATTTGTTAAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGCGt              >  1:1246257/1‑70 (MQ=255)
                                                         ggTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTa                                       <  1:4523446/42‑1 (MQ=255)
                                                         ggTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACg          <  1:1260752/71‑1 (MQ=255)
                                                          gTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    <  1:738768/60‑1 (MQ=255)
                                                          gTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGa         <  1:6937821/71‑1 (MQ=255)
                                                           ccaTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGtt                      <  1:5191618/55‑1 (MQ=38)
                                                           ccaTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGtt                      <  1:5200994/55‑1 (MQ=38)
                                                           ccaTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGtt                      <  1:4781571/55‑1 (MQ=38)
                                                             aTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    >  1:2060981/1‑57 (MQ=255)
                                                             aTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    >  1:3200027/1‑57 (MQ=255)
                                                             aTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    >  1:4634282/1‑57 (MQ=255)
                                                              tttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCt      <  1:5333106/70‑1 (MQ=255)
                                                              tttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCt      <  1:6351604/70‑1 (MQ=255)
                                                              tttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATT‑GGGTTACGACCTa     >  1:5921977/1‑70 (MQ=255)
                                                               ttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGt              <  1:1995570/61‑1 (MQ=255)
                                                               ttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGt              <  1:1726977/61‑1 (MQ=255)
                                                               ttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTa     >  1:32734/1‑70 (MQ=255)
                                                                tGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTAt                                      >  1:1447419/1‑36 (MQ=255)
                                                                tGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTTGTTATTGGGGTTACGACCTa     >  1:6742634/1‑69 (MQ=255)
                                                                tGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTAc    <  1:382965/70‑1 (MQ=255)
                                                                tGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTAc    <  1:4194662/70‑1 (MQ=255)
                                                                 gTTCAGCAAATGGCCCATTTTCTCGGGTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTAcaa  >  1:1038347/1‑71 (MQ=255)
                                                                 gTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    >  1:7315210/1‑53 (MQ=255)
                                                                 gTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTAcaa  >  1:4249525/1‑71 (MQ=255)
                                                                 gTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGGTACGACCTAca   <  1:5256991/70‑1 (MQ=255)
                                                                    |                                                                   
TATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTACAA  >  W3110S.gb/1340225‑1340360

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: ‑