Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,375,910 T→C 8.2% N68N (AAT→AAC groL → Cpn60 chaperonin GroEL, large subunit of GroESL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,375,9100TC 8.2% 163.6 / 5.0 61N68N (AAT→AACgroLCpn60 chaperonin GroEL, large subunit of GroESL
Reads supporting (aligned to +/- strand):  ref base T (43/13);  new base C (2/3);  total (45/16)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

GTGCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCAC  >  W3110S.gb/4375840‑4375978
                                                                      |                                                                    
gTGCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAat                                                                      <  1:2829045/71‑1 (MQ=255)
gTGCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAat                                                                      <  1:1392036/71‑1 (MQ=255)
gTGCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAat                                                                      <  1:3268812/71‑1 (MQ=255)
gTGCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAat                                                                      <  1:1758619/71‑1 (MQ=255)
        aCCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTgc                                                              >  1:117504/1‑71 (MQ=255)
            tCACCAAAGATGGTGTTTCCGTTTCTCGTGAAATCGAACTTGAAGACAAGTTAGAAAatat                                                                    <  1:3720998/61‑1 (MQ=255)
            tCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAatat                                                                    <  1:671510/61‑1 (MQ=255)
              aCCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGa                                                         <  1:142581/70‑1 (MQ=255)
              aCCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGAt                                                        <  1:84140/71‑1 (MQ=255)
                cAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAata                                                                     >  1:434547/1‑56 (MQ=255)
                cAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAata                                                                     >  1:1437428/1‑56 (MQ=255)
                  aaGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAata                                                                     >  1:3459367/1‑54 (MQ=255)
                  aaGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAata                                                                     >  1:2967734/1‑54 (MQ=255)
                  aaGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGTGCGCAg                                                          >  1:834698/1‑65 (MQ=255)
                  aaGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGTGCGCAg                                                          >  1:2231987/1‑65 (MQ=255)
                     aTGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGa                                                   >  1:2740555/1‑69 (MQ=255)
                      tGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAaag                                                <  1:3009832/71‑1 (MQ=255)
                      tGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAaag                                                <  1:1427927/71‑1 (MQ=255)
                        gtgtTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGa                                                         >  1:423832/1‑60 (MQ=255)
                        gtgtTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAaaga                                               <  1:2375206/70‑1 (MQ=255)
                            ttCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTg                                          >  1:3272124/1‑71 (MQ=255)
                             tCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAaagaa                                              <  1:137516/66‑1 (MQ=255)
                             tCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAaagaa                                              <  1:535827/66‑1 (MQ=255)
                                 ttGCTCTTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAaaga                                               <  1:2807944/61‑1 (MQ=255)
                                 ttGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCctc                                      >  1:3079814/1‑70 (MQ=255)
                                 ttGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCctc                                      >  1:151070/1‑70 (MQ=255)
                                  tGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAaagaag                                             >  1:2321706/1‑62 (MQ=255)
                                  tGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAaagaag                                             >  1:3409206/1‑62 (MQ=255)
                                  tGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAaagaag                                             >  1:1805706/1‑62 (MQ=255)
                                  tGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAaagaag                                             >  1:889901/1‑62 (MQ=255)
                                      cGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAg                                 >  1:1828786/1‑70 (MQ=255)
                                      cGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAg                                 >  1:3793502/1‑70 (MQ=255)
                                       gTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGc                                >  1:1784192/1‑70 (MQ=255)
                                        tGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCa                               >  1:315109/1‑70 (MQ=255)
                                        tGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCa                               >  1:3329936/1‑70 (MQ=255)
                                        tGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCa                               >  1:428645/1‑70 (MQ=255)
                                         gAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCaaa                             >  1:738253/1‑71 (MQ=255)
                                          aaaTCGAACTGGAGGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAac                            >  1:2723486/1‑71 (MQ=255)
                                          aaaTCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAac                            >  1:2805474/1‑71 (MQ=255)
                                          aaaTCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAac                            >  1:3300775/1‑71 (MQ=255)
                                          aaaTCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAac                            >  1:3350954/1‑71 (MQ=255)
                                             tCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTCCCTCTAAAGCAAacgac                         >  1:3582531/1‑71 (MQ=255)
                                              cGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAacgacg                        >  1:1773366/1‑71 (MQ=255)
                                              cGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAacgacg                        >  1:1098934/1‑71 (MQ=255)
                                                 aCTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTg                     >  1:1745502/1‑71 (MQ=255)
                                                   tGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGtt                                           <  1:3056217/47‑1 (MQ=255)
                                                   tGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGc                                         <  1:2197761/49‑1 (MQ=255)
                                                     gAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCa                   >  1:77550/1‑69 (MQ=255)
                                                     gAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAgg                 >  1:746888/1‑71 (MQ=255)
                                                     gAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAgg                 >  1:1656899/1‑71 (MQ=255)
                                                     gAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAgg                 >  1:1523315/1‑71 (MQ=255)
                                                     gAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAg                  >  1:2772585/1‑70 (MQ=255)
                                                     gAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCT‑‑AGCAAACGACGCTGc                    >  1:902457/1‑66 (MQ=255)
                                                       aGACAAGTTCGAAAACATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCGAACGACGCTGCAGGCg               <  1:3622824/71‑1 (MQ=255)
                                                       aGACAAGTTCGAAAACATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCGAACGACGCTGCAGGCg               <  1:1706684/71‑1 (MQ=255)
                                                       aGACAAGTTCGAAAACATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCGAACGACGCTGCAGGCg               <  1:643757/71‑1 (MQ=255)
                                                        gACAAGTTCGAAAATATGGGTGCGCAGATGGTGAaaga                                               <  1:3206116/38‑1 (MQ=255)
                                                          cAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGt                                            >  1:3712648/1‑39 (MQ=255)
                                                          cAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGt                                            >  1:1348253/1‑39 (MQ=255)
                                                          cAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGt                                            >  1:157773/1‑39 (MQ=255)
                                                          cAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGt                                            >  1:2289072/1‑39 (MQ=255)
                                                          cAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGt                                            >  1:3003225/1‑39 (MQ=255)
                                                          cAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGt                                            >  1:3005558/1‑39 (MQ=255)
                                                          cAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGAc             <  1:219977/70‑1 (MQ=255)
                                                            aGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGt          <  1:1892772/71‑1 (MQ=255)
                                                              ttCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTac        >  1:1120003/1‑71 (MQ=255)
                                                                cGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAg                                 >  1:2756680/1‑44 (MQ=255)
                                                                    aaTATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGt          >  1:2551543/1‑63 (MQ=255)
                                                                    aaTATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACcaccac  >  1:2397479/1‑71 (MQ=255)
                                                                    aaTATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACcaccac  >  1:1464848/1‑71 (MQ=255)
                                                                    aaTATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACcaccac  >  1:895741/1‑71 (MQ=255)
                                                                      |                                                                    
GTGCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCAC  >  W3110S.gb/4375840‑4375978

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: