Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 257,873 A→G 24.6% K15K (AAA→AAG crl → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb257,8730AG24.6% 100.8 / 36.3 57K15K (AAA→AAGcrlDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (24/19);  new base G (6/8);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT  >  W3110S.gb/257803‑257941
                                                                      |                                                                    
cAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCaaaaaa                                                                      <  1:2152595/71‑1 (MQ=255)
  gTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAAt                                                                     >  1:2913965/1‑70 (MQ=255)
   ttGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAAttt                                                                   <  1:849580/71‑1 (MQ=255)
    tGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAAttt                                                                   >  1:2765204/1‑70 (MQ=255)
     gCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAAttt                                                                   >  1:4848915/1‑69 (MQ=255)
        tCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGc                                                              <  1:2643783/71‑1 (MQ=255)
          acaacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCAc                                                            >  1:4315255/1‑71 (MQ=255)
          acaacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCAc                                                            >  1:29592/1‑71 (MQ=255)
          acaacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCAc                                                            >  1:1347119/1‑71 (MQ=255)
          acaacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCAc                                                            >  1:4442179/1‑71 (MQ=255)
          acaacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCAc                                                            <  1:4436345/71‑1 (MQ=255)
          acaacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCAc                                                            <  1:2971650/71‑1 (MQ=255)
             acaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCAc                                                            <  1:1577827/68‑1 (MQ=255)
               aGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTa                                                          >  1:2028349/1‑68 (MQ=255)
                   gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTAc                                                                 <  1:871663/57‑1 (MQ=255)
                    aTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCg                                                    <  1:1111752/69‑1 (MQ=255)
                    aTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGta                                                  <  1:1805761/71‑1 (MQ=255)
                    aTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGta                                                  <  1:1322813/71‑1 (MQ=255)
                      aGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCACTAGGCCCGtata                                                >  1:4814511/1‑71 (MQ=255)
                      aGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCACTAGGCCCGtata                                                >  1:420628/1‑71 (MQ=255)
                      aGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCACTAGGCCCGtata                                                >  1:2583770/1‑71 (MQ=255)
                      aGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGtata                                                >  1:2024303/1‑71 (MQ=255)
                        cAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCaaaaaa                                                                      >  1:1906928/1‑47 (MQ=255)
                             aCGTTACCGAGTGGACACCCGGAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTg                                          <  1:712350/70‑1 (MQ=255)
                             aCGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCAc                                                            >  1:1433922/1‑52 (MQ=255)
                               gTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAg                                       >  1:1143737/1‑71 (MQ=255)
                                ttACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTAc                                                                 <  1:4474392/44‑1 (MQ=255)
                                ttACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGc                                                              <  1:8082/47‑1 (MQ=255)
                                  aCCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCACTAGGCCCGTATATTCGTGAAGGt                                     <  1:3508095/70‑1 (MQ=255)
                                  aCCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCa                                                             >  1:1639980/1‑46 (MQ=255)
                                  aCCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGt                                           >  1:4333999/1‑64 (MQ=255)
                                  aCCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAgg                                      >  1:1766416/1‑69 (MQ=255)
                                  aCCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGt                                     <  1:453274/70‑1 (MQ=255)
                                   ccGAGTGGACACCCGAAGAGCAGATTGATCaaaaaa                                                                      <  1:5027164/36‑1 (MQ=255)
                                   ccGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTa                                    <  1:3360725/70‑1 (MQ=255)
                                          gACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCaa                             <  1:3424255/70‑1 (MQ=255)
                                           acacCCGAAGAGCAGATTGATCAAAAAGTTTACCg                                                               <  1:3775333/35‑1 (MQ=255)
                                            cacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGtata                                                >  1:2892382/1‑49 (MQ=255)
                                            cacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGa                          <  1:1090905/71‑1 (MQ=255)
                                              cccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGa                                         <  1:1841357/54‑1 (MQ=255)
                                                cGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAg                                       >  1:4315422/1‑54 (MQ=255)
                                                  aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCg                    <  1:4394868/71‑1 (MQ=255)
                                                    gagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGAt                  >  1:855837/1‑71 (MQ=255)
                                                       cAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGt                                           >  1:1330242/1‑43 (MQ=255)
                                                       cAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGt                                           >  1:4399928/1‑43 (MQ=255)
                                                        aGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGAt                         <  1:1346279/60‑1 (MQ=255)
                                                        aGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGAt                         <  1:575862/60‑1 (MQ=255)
                                                         gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCttt             >  1:534152/1‑71 (MQ=255)
                                                         gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCttt             >  1:1681172/1‑71 (MQ=255)
                                                          aTTGATCAAAAAGTTTACCGCACTAGGCCCGTATATTCGTGAAg                                       <  1:2077916/44‑1 (MQ=255)
                                                          aTTGATCAAAAAGTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGa                          <  1:471685/57‑1 (MQ=255)
                                                          aTTGATCAAAAAGTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGa                          <  1:2965618/57‑1 (MQ=255)
                                                          aTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGt                                 <  1:1649728/50‑1 (MQ=255)
                                                            tGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTaa                                   >  1:3665616/1‑46 (MQ=255)
                                                            tGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTaa                                   >  1:2262668/1‑46 (MQ=255)
                                                            tGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGc                               >  1:1459890/1‑50 (MQ=255)
                                                            tGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTc          >  1:4051291/1‑71 (MQ=255)
                                                            tGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTc          >  1:283536/1‑71 (MQ=255)
                                                              aTCAAAAAGTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGa        <  1:1370589/71‑1 (MQ=255)
                                                                cAAAAAATTTACCGCACTAGGCCCGTATATTCGTg                                          >  1:333928/1‑35 (MQ=255)
                                                                 aaaaagTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTg     <  1:1596701/71‑1 (MQ=255)
                                                                 aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGa        >  1:1116871/1‑68 (MQ=255)
                                                                 aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGAt       >  1:5016235/1‑69 (MQ=255)
                                                                  aaaagTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTg     <  1:3940715/70‑1 (MQ=255)
                                                                  aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTc          <  1:1107880/65‑1 (MQ=255)
                                                                  aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTc          <  1:3693695/65‑1 (MQ=255)
                                                                  aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGt    <  1:355908/71‑1 (MQ=255)
                                                                    aaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCt  <  1:534452/71‑1 (MQ=255)
                                                                      |                                                                    
CAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT  >  W3110S.gb/257803‑257941

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: