Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,461,557 A→C 8.9% intergenic (‑465/+216) yfcV ← / ← sixA predicted fimbrial‑like adhesin protein/phosphohistidine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,461,5570AC 8.9% 399.0 / 38.7 146intergenic (‑465/+216)yfcV/sixApredicted fimbrial‑like adhesin protein/phosphohistidine phosphatase
Reads supporting (aligned to +/- strand):  ref base A (58/75);  new base C (10/3);  total (68/78)
Fisher's exact test for biased strand distribution p-value = 3.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

Only 100 of 170 total aligned reads displayed.

TAATCAAGTTTCAATTACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGATACGCG  >  W3110S.gb/2461491‑2461625
                                                                   |                                                                    
tAATCAAGTTTCAATTACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGa                                                                  >  1:1335253/1‑71 (MQ=255)
  aTCAAGTTTCAATTACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGAt                                                                 <  1:4874750/70‑1 (MQ=255)
     aaGTTTCAATTACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              <  1:3356251/70‑1 (MQ=255)
     aaGTTTCAATTACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              <  1:1663556/70‑1 (MQ=255)
     aaGTTTCAATTACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              <  1:2624347/70‑1 (MQ=255)
     aaGTTTCAATTACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              <  1:2801843/70‑1 (MQ=255)
      aGTTTCAATTACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATaa                                                               >  1:2267649/1‑68 (MQ=255)
      aGTTTCAATTACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              >  1:4825695/1‑69 (MQ=255)
         ttCAATTACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAgtg                                                         >  1:5307921/1‑71 (MQ=255)
         ttCAATTACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAgtg                                                         >  1:3036264/1‑71 (MQ=255)
              ttACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              >  1:3149517/1‑61 (MQ=255)
              ttACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              >  1:693737/1‑61 (MQ=255)
              ttACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              >  1:1719436/1‑61 (MQ=255)
              ttACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              >  1:4196019/1‑61 (MQ=255)
              ttACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCaaa                                                    <  1:4311519/71‑1 (MQ=255)
                     ctgcctgATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              <  1:2377616/54‑1 (MQ=255)
                     ctgcctgATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACat                                             >  1:4286315/1‑71 (MQ=255)
                       gcctgATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACatta                                           >  1:3709535/1‑71 (MQ=255)
                       gcctgATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACatta                                           >  1:4494524/1‑71 (MQ=255)
                       gcctgATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACatta                                           >  1:908650/1‑71 (MQ=255)
                       gcctgATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACatta                                           >  1:658953/1‑71 (MQ=255)
                            attattCAACTCTATTTTCTCCATTATCTGTATTAATaaa                                                                     >  1:3200850/1‑40 (MQ=255)
                            attattCAACTCTATTTTCTCCATTATCTGTATTAATAAATGAt                                                                 <  1:5375587/44‑1 (MQ=255)
                            attattCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGa                                                >  1:483973/1‑61 (MQ=255)
                            attattCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCa                                       <  1:2693216/70‑1 (MQ=255)
                            attattCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCAt                                      >  1:765758/1‑71 (MQ=255)
                             ttattCAACTCTATTTTCTCCATTATCTGTATTAATaaa                                                                     >  1:3473637/1‑39 (MQ=255)
                             ttattCAACTCTATTTTCTCCATTATCTGTATTAATaaa                                                                     >  1:4817530/1‑39 (MQ=255)
                             ttattCAACTCTATTTTCTCCATTATCTGTATTAATaaa                                                                     >  1:1493686/1‑39 (MQ=255)
                             ttattCAACTCTATTTTCTCCATTATCTGTATTAATaaa                                                                     >  1:872409/1‑39 (MQ=255)
                             ttattCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              <  1:4235927/46‑1 (MQ=255)
                             ttattCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              <  1:2857057/46‑1 (MQ=255)
                             ttattCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              <  1:325440/46‑1 (MQ=255)
                              tattCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAAt                                                              <  1:2218769/45‑1 (MQ=255)
                                ttaaaCTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTc                                   <  1:182247/67‑1 (MQ=255)
                                ttCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTc                                   <  1:488794/70‑1 (MQ=255)
                                  cAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTc                                   <  1:4019074/68‑1 (MQ=255)
                                   aaCTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAg                                                           >  1:5165177/1‑43 (MQ=255)
                                   aaCTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAg                                                           >  1:4253619/1‑43 (MQ=255)
                                   aaCTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAg                                                           >  1:5310601/1‑43 (MQ=255)
                                   aaCTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACa                               >  1:1223070/1‑71 (MQ=255)
                                   aaCTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACa                               >  1:570999/1‑71 (MQ=255)
                                   aaCTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACa                               >  1:554816/1‑71 (MQ=255)
                                      tctATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCaa                                 <  1:1049573/66‑1 (MQ=255)
                                      tctATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCaa                                 <  1:5171574/66‑1 (MQ=255)
                                      tctATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCaa                                 <  1:2520752/66‑1 (MQ=255)
                                      tctATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCaa                                 <  1:2492177/66‑1 (MQ=255)
                                      tctATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCaa                                 <  1:944629/66‑1 (MQ=255)
                                        tataTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGt                            <  1:3903582/69‑1 (MQ=255)
                                        tATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCaaa                                                    >  1:4761688/1‑45 (MQ=255)
                                        tATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCaaa                                                    >  1:3262055/1‑45 (MQ=255)
                                        tATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGt                            <  1:4808971/69‑1 (MQ=255)
                                        tATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGt                            <  1:2650901/69‑1 (MQ=255)
                                        tATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGt                            <  1:4023133/69‑1 (MQ=255)
                                        tATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGt                            <  1:2106066/69‑1 (MQ=255)
                                        tATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTaa                          >  1:4541933/1‑71 (MQ=255)
                                          ttttCTCCATTATCTGTATTAATAAATGATAATGAg                                                           >  1:5492885/1‑36 (MQ=255)
                                          ttttCTCCATTATCTGTATTAATAAATGATAATGAg                                                           >  1:2006158/1‑36 (MQ=255)
                                          ttttCTCCATTATCTGTATTAATAAATGATAATGAg                                                           >  1:3935027/1‑36 (MQ=255)
                                          ttttCTCCATTATCTGTATTAATAAATGATAATGAg                                                           >  1:5051377/1‑36 (MQ=255)
                                          ttttCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTaa                          <  1:1220853/69‑1 (MQ=255)
                                            ttCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACttt                      >  1:790459/1‑71 (MQ=255)
                                            ttCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACttt                      >  1:5470755/1‑71 (MQ=255)
                                            ttCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACttt                      >  1:3006118/1‑71 (MQ=255)
                                            ttCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACtt                       >  1:1316655/1‑70 (MQ=255)
                                             tctcCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACttt                      <  1:988045/70‑1 (MQ=255)
                                             tctcCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACttt                      <  1:2328211/70‑1 (MQ=255)
                                              ctcCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTg                     <  1:1518749/70‑1 (MQ=255)
                                              ctcCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTg                     <  1:2834562/70‑1 (MQ=255)
                                              ctcCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTg                     <  1:4940258/70‑1 (MQ=255)
                                              ctcCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTg                     <  1:3388739/70‑1 (MQ=255)
                                               ctcATTATCATTTATTAATACAGATAATGAGTGTCAAACAGACATTATTCAtt                                     <  1:105042/51‑1 (MQ=9)
                                                 cATTATCATTTATTAATACAGATAATGAGTGTCAAACAGACATTATTCATTTc                                   >  1:2488584/1‑53 (MQ=17)
                                                 cATTATCATTTATTAATACAGATAATGAGTGTCAAACAGACATTATTCATTTc                                   >  1:1095640/1‑53 (MQ=17)
                                                 cATTATCATTTATTAATACAGATAATGAGTGTCAAACAGACATTATTCATTTc                                   >  1:796289/1‑53 (MQ=17)
                                                 cATTATCATTTATTAATACAGATAATGAGTGTCAAACAGACATTATTCATTTc                                   >  1:1996156/1‑53 (MQ=17)
                                                 cATTATCATTTATTAATACAGATAATGAGTGTCAAACAGACATTATTCATTTc                                   >  1:2384276/1‑53 (MQ=17)
                                                 cATTATCATTTATTAATACAGATAATGAGTGTCAAACAGACATTATTCATTTc                                   >  1:3256095/1‑53 (MQ=17)
                                                   ttATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACt                        <  1:5203750/62‑1 (MQ=255)
                                                   ttATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACt                        <  1:1427982/62‑1 (MQ=255)
                                                   ttATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACt                        <  1:762941/62‑1 (MQ=255)
                                                     aTCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACtt                       >  1:4267791/1‑61 (MQ=255)
                                                     aTCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAAt             >  1:4060613/1‑71 (MQ=255)
                                                     aTCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAAt             >  1:4659227/1‑71 (MQ=255)
                                                           tATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGa        <  1:3981140/71‑1 (MQ=255)
                                                           tATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGa        <  1:82844/71‑1 (MQ=255)
                                                           tATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGa        <  1:1187555/71‑1 (MQ=255)
                                                            aTTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGa        <  1:3372096/70‑1 (MQ=255)
                                                            aTTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGa        <  1:4788494/70‑1 (MQ=255)
                                                               aataaATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGa        >  1:5315743/1‑67 (MQ=255)
                                                               aataaATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGATaag    >  1:3539043/1‑69 (MQ=255)
                                                               aataaATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGATAcg    >  1:5449089/1‑71 (MQ=255)
                                                               aataaATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGATAcg    >  1:2430312/1‑71 (MQ=255)
                                                               aataaATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGATAcg    >  1:2127594/1‑71 (MQ=255)
                                                               aataaATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGATAcg    >  1:4755290/1‑71 (MQ=255)
                                                               aataaATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGATAcg    >  1:1454649/1‑71 (MQ=255)
                                                               aataaATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGATAcg    >  1:2554500/1‑71 (MQ=255)
                                                               aataaATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGATAcg    >  1:3426151/1‑71 (MQ=255)
                                                               aataaATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGATAc     >  1:153946/1‑70 (MQ=255)
                                                                   aaTGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGATAcgcg  >  1:468958/1‑69 (MQ=255)
                                                                   |                                                                    
TAATCAAGTTTCAATTACGTCCTGCCTGATTATTCAACTCTATTTTCTCCATTATCTGTATTAATAAATGATAATGAGTGTCAAACAGACATTATTCATTTCAACATGTAACTTTGGGCAAAATCAAGATACGCG  >  W3110S.gb/2461491‑2461625

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: