Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,678,974 C→A 17.8% R38R (CGC→CGA ilvY → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,678,9740CA17.8% 77.5 / 4.8 28R38R (CGC→CGAilvYDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (11/12);  new base A (2/3);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CATTTTGGCCGCAGCGCGCGGGCGATGCACGTTAGCCCATCCACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCACGGTGACGCTGACTGAAGCGGGCG  >  W3110S.gb/3678909‑3679044
                                                                 |                                                                      
cATTTTGGCCGCAGCGAGCGGGCGATGCACGTTAGCCCATCCACGCTCTCACGGCAGATTCAGCGCCTgg                                                                    >  1:2345973/1‑70 (MQ=255)
            aGCGCGCGGGCGATGCACGTTAGCCCATCCACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCgg                                                       <  1:4331501/71‑1 (MQ=255)
                       gATGCACGTTAGCCCATCCACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGt                                            <  1:1270794/71‑1 (MQ=255)
                         tggaCGTTAGCCCATCCACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGc                                               <  1:1847670/63‑1 (MQ=255)
                         tGCACGTTAGCCCATCCACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGttt                                          >  1:614499/1‑71 (MQ=255)
                             cGTTAGCCCATCCACGCTCTCACGGCAGATTCAGCGACTgg                                                                    <  1:3055685/41‑1 (MQ=255)
                               ttAGCCCATCCACGCTCTCACGGCAGATTCAGCGACTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTgcgc                                    >  1:1260772/1‑71 (MQ=255)
                                   cccATCCACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGAt                                 <  1:2833974/70‑1 (MQ=255)
                                     cATCCACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTgcgc                                    <  1:41940/65‑1 (MQ=255)
                                      aTCCACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAg                                                   <  1:218662/49‑1 (MQ=255)
                                       tCCACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGAt                                                            <  1:2889497/39‑1 (MQ=255)
                                       tCCACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCt                                              <  1:585217/53‑1 (MQ=255)
                                         cACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGAt                                                            >  1:1937372/1‑37 (MQ=255)
                                         cACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTtgtg                                       >  1:1085217/1‑58 (MQ=255)
                                         cACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCa                          >  1:1145701/1‑71 (MQ=255)
                                         cACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCa                          >  1:3624607/1‑71 (MQ=255)
                                         cACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCa                          >  1:1913951/1‑71 (MQ=255)
                                         cACGCTCTCACGGCAGATTCAGCGACTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGc                           <  1:4202721/70‑1 (MQ=255)
                                         cACGCTCTCACGGCAGATTCAGCGACTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCa                          >  1:3216770/1‑71 (MQ=255)
                                          aCGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCAc                         <  1:230704/71‑1 (MQ=255)
                                             ctctCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCACGGt                      <  1:846818/71‑1 (MQ=255)
                                             ctctCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCACGGt                      <  1:4437508/71‑1 (MQ=255)
                                              tctcACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCACGGt                      <  1:1396759/70‑1 (MQ=255)
                                               ctcACGGCAGATTCAGCGACTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCACGGt                      <  1:2949550/69‑1 (MQ=255)
                                                       aGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCACGGTGACGctgact            >  1:3684178/1‑71 (MQ=255)
                                                           tCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCACGGTGa                    >  1:1488669/1‑59 (MQ=255)
                                                           tCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGAGCGATAACCGCACGGTGACGCTGACTGAAg        >  1:3096985/1‑71 (MQ=255)
                                                            cAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCACGGTGACGCTGACTGAAGc       >  1:766752/1‑71 (MQ=255)
                                                              gcgcCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATa                                <  1:4011548/44‑1 (MQ=255)
                                                               cgcCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCACGGTGACGCTGACTGAAGCggg    >  1:3920268/1‑71 (MQ=255)
                                                                gcCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATa                                <  1:82027/42‑1 (MQ=255)
                                                                gcCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATa                                <  1:87430/42‑1 (MQ=255)
                                                                 ccTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCACGGTGACGCTGACTGAAGCGGGCg  <  1:3616308/71‑1 (MQ=255)
                                                                 |                                                                      
CATTTTGGCCGCAGCGCGCGGGCGATGCACGTTAGCCCATCCACGCTCTCACGGCAGATTCAGCGCCTGGAAGAAGATCTCGGTCAGCCGCTGTTTGTGCGCGATAACCGCACGGTGACGCTGACTGAAGCGGGCG  >  W3110S.gb/3678909‑3679044

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: