Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,013,681 G→A 59.6% V11V (GTC→GTT yqeC ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,013,6810GA59.6% 29.4 / 61.5 52V11V (GTC→GTTyqeCconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (12/9);  new base A (12/19);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTA  >  W3110S.gb/3013632‑3013743
                                                 |                                                              
cACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAAc                                            >  1:988823/1‑70 (MQ=255)
cACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAAc                                            >  1:697898/1‑70 (MQ=255)
cACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAAc                                            >  1:618303/1‑70 (MQ=255)
cACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAAc                                            >  1:363060/1‑70 (MQ=255)
cACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAAc                                            >  1:332210/1‑70 (MQ=255)
aaCCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAAc                                            >  1:1963334/2‑70 (MQ=255)
  ccAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACta                                          >  1:317055/1‑70 (MQ=255)
  ccAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACta                                          >  1:330685/1‑70 (MQ=255)
        acaGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACtttt                                   >  1:54557/1‑71 (MQ=255)
        acaGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACttt                                    >  1:1064891/1‑70 (MQ=255)
        acaGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACtttt                                   <  1:882180/71‑1 (MQ=255)
        acaGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACtttt                                   <  1:794196/71‑1 (MQ=255)
        acaGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACtttt                                   <  1:354171/71‑1 (MQ=255)
              aTCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGa                              >  1:1958915/1‑70 (MQ=255)
              aTCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAg                             >  1:1625049/1‑71 (MQ=255)
              aTCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAg                             >  1:1882367/1‑71 (MQ=255)
              aTCACAGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAg                             >  1:1039698/1‑71 (MQ=255)
                  caGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCtt                         <  1:333803/71‑1 (MQ=255)
                  caGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCtt                         <  1:1651474/71‑1 (MQ=255)
                  caGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCtt                         <  1:1771640/71‑1 (MQ=255)
                   aGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCttt                        <  1:869489/71‑1 (MQ=255)
                   aGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCttt                        <  1:294108/71‑1 (MQ=255)
                   aGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCttt                        <  1:997427/71‑1 (MQ=255)
                   aGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCttt                        <  1:1188756/71‑1 (MQ=255)
                   aGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCttt                        <  1:1764295/71‑1 (MQ=255)
                   aGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCttt                        <  1:1586092/71‑1 (MQ=255)
                   aGTGGGACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCtt                         <  1:666200/70‑1 (MQ=255)
                       ggACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACtttt                                   <  1:103523/56‑1 (MQ=255)
                        gACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGAc                   <  1:1114283/71‑1 (MQ=255)
                        gACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGAc                   <  1:986504/71‑1 (MQ=255)
                        gACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGAc                   <  1:1640933/71‑1 (MQ=255)
                        gACGTTTCTGCGCACCTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGAc                   <  1:7376/71‑1 (MQ=255)
                         aCGTTTCTGCGCACCTAAATCAATGACTAATGCAGATGGGTCAACTATACTTTTCACGAGTCttt                        <  1:1243444/65‑1 (MQ=255)
                            tttCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACTAGTCTTTATGACctct               >  1:352217/1‑71 (MQ=255)
                            tttCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACctct               >  1:602087/1‑71 (MQ=255)
                            tttCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACctct               >  1:810298/1‑71 (MQ=255)
                            tttCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACctct               >  1:1570791/1‑71 (MQ=255)
                            tttCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACctct               >  1:1593717/1‑71 (MQ=255)
                            tttCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACctct               >  1:313543/1‑71 (MQ=255)
                            tttCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACctc                >  1:1226413/1‑70 (MQ=255)
                                       cTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAAt    <  1:335244/71‑1 (MQ=255)
                                       cTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAAt    <  1:1818313/71‑1 (MQ=255)
                                       cTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGAc                   <  1:1919503/56‑1 (MQ=255)
                                       cTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGAc                   <  1:504405/56‑1 (MQ=255)
                                       cTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGAc                   <  1:82325/56‑1 (MQ=255)
                                       cTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGAc                   <  1:1175705/56‑1 (MQ=255)
                                       cTAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGAc                   <  1:93137/56‑1 (MQ=255)
                                        tAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAAt    >  1:1626177/1‑70 (MQ=255)
                                        tAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAAt    >  1:1789902/1‑70 (MQ=255)
                                        tAAATCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGAc                   <  1:1538313/55‑1 (MQ=255)
                                         aaaTCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTa  >  1:219641/1‑71 (MQ=255)
                                           aTCAATAACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAAtt   <  1:790295/68‑1 (MQ=255)
                                                 |                                                              
CACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTA  >  W3110S.gb/3013632‑3013743

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: