Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 75,677 A→C 23.9% F541L (TTT→TTG sgrR ← DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb75,6770AC23.9% 51.0 / 34.6 46F541L (TTT→TTGsgrRDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (20/15);  new base C (8/3);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

GGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAG  >  W3110S.gb/75609‑75746
                                                                    |                                                                     
ggCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTaaa                                                                     <  1:2987271/71‑1 (MQ=255)
ggCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTaaa                                                                     <  1:957833/71‑1 (MQ=255)
ggCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTaaa                                                                     <  1:2898920/71‑1 (MQ=255)
ggCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTaaa                                                                     <  1:1740588/71‑1 (MQ=255)
          tagtgATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTCAAATCGAACCAg                                                           >  1:1075109/1‑71 (MQ=255)
            gtgATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGcc                                                         >  1:674477/1‑71 (MQ=255)
            gtgATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGcc                                                         >  1:1368802/1‑71 (MQ=255)
              gATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGcc                                                         >  1:3476944/1‑69 (MQ=255)
              gATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGcc                                                         >  1:1564170/1‑69 (MQ=255)
              gATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGcc                                                         >  1:3381735/1‑69 (MQ=255)
                ttttGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGa                                                       <  1:3169166/69‑1 (MQ=255)
                 tttGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTaaaa                                                                    >  1:1453418/1‑55 (MQ=255)
                  ttGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAggg                                                    <  1:1724905/70‑1 (MQ=255)
                  ttGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAggg                                                    <  1:188837/70‑1 (MQ=255)
                  ttGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAggg                                                    >  1:2133658/1‑70 (MQ=255)
                  ttGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAggg                                                    <  1:228499/70‑1 (MQ=255)
                  ttGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAggg                                                    >  1:408769/1‑70 (MQ=255)
                  ttGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCAAGCTGATTTAAAATCGAACCAGCCGAggg                                                    <  1:3810708/70‑1 (MQ=255)
                     ttATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTAtt                                                <  1:124722/71‑1 (MQ=255)
                       aTCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTCAAATCGAACCAGCCGAGGGTa                                                  >  1:1302610/1‑67 (MQ=255)
                          aaCCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTaaaa                                                                    >  1:2600524/1‑46 (MQ=255)
                          aaCCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTaaaa                                                                    >  1:2311334/1‑46 (MQ=255)
                               gCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTCAAATCGAACCAGCCGAGGGTATTCa                                              >  1:1306164/1‑63 (MQ=255)
                               gCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTCAAATCGAACCAGCCGAGGGTATTCa                                              >  1:4022976/1‑63 (MQ=255)
                               gCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTCAAATCGAACCAGCCGAGGGTATTCa                                              >  1:4027799/1‑63 (MQ=255)
                                   tCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCg                                                                 <  1:509877/40‑1 (MQ=255)
                                   tCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCg                                                                 <  1:1405015/40‑1 (MQ=255)
                                    cATGGATCCGGTGGCGCAAACCACGCTGATTTaaa                                                                     <  1:764817/35‑1 (MQ=255)
                                    cATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCcgcg                                 >  1:2752212/1‑71 (MQ=255)
                                     aTGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAAc                                                              >  1:1819758/1‑41 (MQ=255)
                                     aTGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAAc                                                              >  1:3333361/1‑41 (MQ=255)
                                     aTGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAAc                                                              >  1:3412211/1‑41 (MQ=255)
                                     aTGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAAc                                                              >  1:1316399/1‑41 (MQ=255)
                                     aTGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCcgcg                                 >  1:856150/1‑70 (MQ=255)
                                      tGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCa                               <  1:2729643/71‑1 (MQ=255)
                                      tGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCa                               <  1:3660938/71‑1 (MQ=255)
                                        gATCCGGTGGCGCAAACCACGCTGATTTCAAATCGAACCAGCCGAGTGTATTCATGCGCAGGCCGCGCATa                             >  1:3420908/1‑71 (MQ=255)
                                        gATCCGGTGGCGCAAACCACGCTGATTTCAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATa                             >  1:1088855/1‑71 (MQ=255)
                                        gATCCGGTGGCGCAAACCACGCTGATTTCAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATa                             >  1:2704860/1‑71 (MQ=255)
                                            cGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGc                         >  1:1895878/1‑71 (MQ=255)
                                              gTGGCGCAAACCACGCTGATTTAAAATCGAACCAGcc                                                         <  1:3408243/37‑1 (MQ=255)
                                                ggCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTg                     >  1:3878749/1‑71 (MQ=255)
                                                ggCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTg                     >  1:710941/1‑71 (MQ=255)
                                                 gcgcAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCa                                              >  1:868845/1‑45 (MQ=255)
                                                 gcgcAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCa                                              >  1:873933/1‑45 (MQ=255)
                                                      aaCCACGCTGATTTCAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTg               <  1:315925/71‑1 (MQ=255)
                                                      aaCCACGCTGATTTCAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTg               <  1:2730286/71‑1 (MQ=255)
                                                      aaCCACGCTGATTTCAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTg               <  1:1534135/71‑1 (MQ=255)
                                                           cGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCa                               <  1:3055638/50‑1 (MQ=255)
                                                           cGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGa          <  1:236489/71‑1 (MQ=255)
                                                           cGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGa          <  1:1898223/71‑1 (MQ=255)
                                                             cTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCcgc                                  <  1:67996/45‑1 (MQ=255)
                                                                 tttAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGc     >  1:2991468/1‑70 (MQ=255)
                                                                   tAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAg  >  1:3906004/1‑71 (MQ=255)
                                                                    |                                                                     
GGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAG  >  W3110S.gb/75609‑75746

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: