Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,172,235 C→T 29.8% pseudogene (42/447 nt)
pseudogene (1636/2041 nt)
gatR ←
gatR ←
DNA‑binding transcriptional regulator, N‑ter fragment; ECK2083:JW5340:b2087
DNA‑binding transcriptional regulator; ECK2083:JW5340+JW2074:b4498

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,172,2350CT29.8% 85.3 / 44.3 57pseudogene (42/447 nt)
pseudogene (1636/2041 nt)
gatR
gatR
DNA‑binding transcriptional regulator, N‑ter fragment; ECK2083:JW5340:b2087
DNA‑binding transcriptional regulator; ECK2083:JW5340+JW2074:b4498
Reads supporting (aligned to +/- strand):  ref base C (30/10);  new base T (9/8);  total (39/18)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACG  >  W3110S.gb/2172175‑2172304
                                                            |                                                                     
cAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACTGTGATAtttt                                                             >  1:698832/1‑71 (MQ=255)
cAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACTGTGATAtttt                                                             >  1:1260147/1‑71 (MQ=255)
cAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACTGTGATAtttt                                                             >  1:1247099/1‑71 (MQ=255)
cAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATAtttt                                                             >  1:1115958/1‑71 (MQ=255)
cAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATAtttt                                                             >  1:1940194/1‑71 (MQ=255)
cAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATAtttt                                                             >  1:764232/1‑71 (MQ=255)
cAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATAtttt                                                             >  1:762041/1‑71 (MQ=255)
cAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATAtttt                                                             >  1:63461/1‑71 (MQ=255)
cAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGT‑‑‑TAATCACCGTGATATTTTTgg                                                          >  1:1211684/1‑71 (MQ=255)
     tAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACTGTGATATTTTTGGCg                                                        >  1:741993/1‑71 (MQ=255)
       aTGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACTGTGATATTTTTGGCGGt                                                      <  1:1600085/71‑1 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCATACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCa                                                    >  1:765580/1‑71 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCACACGCTGCCGGGAGACTGTTGGTAATCACCGTGata                                                                 >  1:794629/1‑58 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCa                                                    >  1:13519/1‑71 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCa                                                    >  1:1233503/1‑71 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCa                                                    >  1:451282/1‑71 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCa                                                    >  1:1072216/1‑71 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCa                                                    >  1:658751/1‑71 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCa                                                    >  1:263051/1‑71 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCa                                                    >  1:229324/1‑71 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCa                                                    >  1:888441/1‑71 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCa                                                    >  1:207532/1‑71 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCa                                                    >  1:112124/1‑71 (MQ=255)
         gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCa                                                    >  1:1926549/1‑71 (MQ=255)
          tCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACTGTGATATTTTTGGCGGTCAt                                                   >  1:124231/1‑71 (MQ=255)
          tCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACTGTGATATTTTTGGCGGTCAt                                                   >  1:1572982/1‑71 (MQ=255)
          tCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGATGGTAATCACTGTGATATTTTTGGCGGTCAt                                                   >  1:1909935/1‑71 (MQ=255)
          tCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGATGGTAATCACTGTGATATTTTTGGCGGTCAt                                                   >  1:302372/1‑71 (MQ=255)
          tCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGATGGTAATCACTGTGATATTTTTGGCGGTCAt                                                   >  1:298644/1‑71 (MQ=255)
            tttATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCAtt                                                  >  1:1153734/1‑70 (MQ=255)
            tttATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCAtt                                                  >  1:1174723/1‑70 (MQ=255)
            tttATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCAtt                                                  >  1:53638/1‑70 (MQ=255)
            tttATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCAtt                                                  >  1:1667296/1‑70 (MQ=255)
               aTTTTCGGAAAGGGCAAACGCTGCCTGGAGACTGTTGGTAATCACTGTGATATTTTTGGCGGTCATTAAt                                               <  1:121835/70‑1 (MQ=255)
                ttttCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACTGTGATATTTTTGGCGGTCATTAATcc                                             <  1:782544/71‑1 (MQ=255)
                ttttCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACTGTGATATTTTTGGCGGTCATTAATcc                                             <  1:649550/71‑1 (MQ=255)
                ttttCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACTGTGATATTTTTGGCGGTCATTAATcc                                             <  1:836233/71‑1 (MQ=255)
                            ggaaaGGCTGCCGGGAGACTGTTGGTAATCAGCGTGTTTTTTTTGGCGGTCATTAATCCTTCGGCTATAAg                                 <  1:780789/69‑1 (MQ=255)
                            gCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAg                                 <  1:825801/71‑1 (MQ=255)
                            gCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAg                                 <  1:1400924/71‑1 (MQ=255)
                            gCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAg                                 <  1:855627/71‑1 (MQ=255)
                             cAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGc                                <  1:1270655/71‑1 (MQ=255)
                                      ccGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAg                                 <  1:1676199/61‑1 (MQ=255)
                                                    gTAATCACTGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCgag         <  1:51998/71‑1 (MQ=255)
                                                    gTAATCACTGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCgag         <  1:1148011/71‑1 (MQ=255)
                                                    gTAATCACTGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCgag         <  1:1395609/71‑1 (MQ=255)
                                                        tCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGc                                <  1:405832/44‑1 (MQ=255)
                                                        tCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGc                                <  1:425269/44‑1 (MQ=255)
                                                        tCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGc                                <  1:1134980/44‑1 (MQ=255)
                                                        tCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGc                                <  1:76117/44‑1 (MQ=255)
                                                          aCCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAAc   >  1:1595809/1‑71 (MQ=255)
                                                          aCCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAAc   >  1:190185/1‑71 (MQ=255)
                                                          aCCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAAc   >  1:1833333/1‑71 (MQ=255)
                                                          aCCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAAc   >  1:1196922/1‑71 (MQ=255)
                                                          aCCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAAc   >  1:634214/1‑71 (MQ=255)
                                                          aCCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAAc   >  1:574127/1‑71 (MQ=255)
                                                          aCCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAAc   >  1:1669155/1‑71 (MQ=255)
                                                           ccGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACg  <  1:977384/71‑1 (MQ=255)
                                                           ccGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACg  <  1:642885/71‑1 (MQ=255)
                                                           ccGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACg  <  1:427103/71‑1 (MQ=255)
                                                           ccGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACg  <  1:1965852/71‑1 (MQ=255)
                                                           ccGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACg  <  1:1796911/71‑1 (MQ=255)
                                                           ccGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACg  <  1:1691354/71‑1 (MQ=255)
                                                           ccGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACg  <  1:155409/71‑1 (MQ=255)
                                                            |                                                                     
CAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACG  >  W3110S.gb/2172175‑2172304

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: